Off-diagonal cosmological solutions in emergent gravity theories and Grigory Perelman entropy for geometric flows

We develop an approach to the theory of relativistic geometric flows and emergent gravity defined by entropy functionals and related statistical thermodynamics models. Nonholonomic deformations of G. Perelman’s functionals and related entropic values used for deriving relativistic geometric evolution flow equations. For self-similar configurations, such equations describe generalized Ricci solitons defining modified Einstein equations. We analyse possible connections between relativistic models of nonholonomic Ricci flows and emergent modified gravity theories. We prove that corresponding systems of nonlinear partial differential equations, PDEs, for entropic flows and modified gravity posses certain general decoupling and integration properties. There are constructed new classes of exact and parametric solutions for nonstationary configurations and locally anisotropic cosmological metrics in modified gravity theories and general relativity. Such solutions describe scenarios of nonlinear geometric evolution and gravitational and matter field dynamics with pattern-forming and quasiperiodic structure and various space quasicrystal and deformed spacetime crystal models. We analyse new classes of generic off-diagonal solutions for entropic gravity theories and show how such solutions can be used for explaining structure formation in modern cosmology. Finally, we speculate why the approaches with Perelman–Lyapunov type functionals are more general or complementary to the constructions elaborated using the concept of Bekenstein–Hawking entropy.


Introduction
The most inspiring ideas in recent development of the gravity theory and cosmology are on the emergent thermodynamic nature of the spacetime geometry, when the Einstein equations can be derived using area-entropy formulas for horizons of black holes, BH, [1][2][3] and from supposed elastic properties of gravity [5][6][7]. A substantial progress includes the research on the microscopic origin of Bekenstein-Hawking entropy in string theory [8], the holographic principle [9,10], BH complementarity [11], and gauge/gravity correspondence [12,13]. Here we note a subsequent development of the (anti) de Sitter, (A)dS, and conformal field theories, CFTs, and AdS/CFT, correspondence [14] and formulating laws of the thermodynamics of 'apparent' horizons [15]. Later, it was proposed that gravity theories and generalized/modified/linearized Einstein equations are consequences of the quantum entanglement connecting nearby spacetime regions [16][17][18][19][20][21]. Recent theoretical activities are devoted to proofs that the entanglement first law and dual gravity are derived from the CFT and/or reveal a deep connection to ideas on emergence of spacetime and gravity from general quantum information principles [4,7,[22][23][24][25]. It was pointed out that in a model of dual gravity with entanglement is equivalent to the full (nonlinear) equations [26]. An intriguing conjecture that gravity links to an entropic force as a spacetime elasticity was proposed by Verlinde [6,7]. It was based on the idea that gravitational interactions result from information regarding the positions of material bodies. Emergent phenomena for gravity were investigated by many other authors and their efforts involve the holographic principle in particle physics and the information theory. There were studied geometric models and possible applications related to gravity and quantum computers; quantum gravity; cosmological inflation and acceleration; and dark energy and dark matter physics etc., see [27][28][29][30][31] and references therein.
In the quest to explore the connection between the models of emergent gravity and in modified gravity theories, MGTs, or general relativity, GR, one involve a strict area law for the BH, (A)dS, or entanglement entropy and further developments for holographic models. To derive gravitational field equations was considered that a small but nonzero volume law entropy would compete with, and at large distances involves, the area law. In certain models this is due to thermalization, elastic spacetime properties, quantum entanglement, holographic effects etc. It was proposed that such a phenomenon occurs in the dS space being responsible for the presence of a cosmological horizon. Nevertheless, in this series of two works, see [32] as a partner paper with complementary results, we deal with quite different issues on geometric flow modifications of gravity and spacetime thermodynamics. It is just shown that relativistic generalizations with a corresponding choice of evolution and thermodynamic functionals support the ideas on the origin of gravity as an effect of the entropic force but with a new type of geometric thermodynamics entropy. Such an approach with Perelman-Lyapunov entropy type functionals was developed in our works on entropic nonholonomic geometric flow evolution, nonlinear dynamics and thermodynamics for relativistic, noncommutative, fractional, supersymmetric, Finsler-Lagrange-Hamilton etc. generalizations of the theory of Ricci flow evolution and applications in modern gravity and cosmology [33][34][35][36][37][38][39][40][41].
The goal of this paper is to elaborate on geometric and physical theories relating relativistic generalizations of the Poincaré-Thurston conjecture, 1 emergent and modified gravity constructions, and Verlinde conjecture that gravity results from an entropic force as a spacetime elasticity which explain fundamental properties of dark matter, DM, and dark energy, DE, in modern cosmology [6,7]. On topology and geometry of Ricci flows, we refer to classical works [42][43][44][45] and reviews of mathematical results in monographs [46][47][48]). Here we note that Friedan published a series of works on nonlinear sigma models, σ -models, in 2 + ε dimensions, see [49][50][51] where geometric flow equations were introduced for the renorm group, RG, theories, see recent results in [52][53][54].
The other goal of this article is to develop the anholonomic frame deformation method, AFDM, (on early works see [55][56][57] and references therein), for constructing exact and parametric quasiperiodic solutions of geometric entropic flow and modified gravity equations. For reviews of recent results on black hole solutions in MGTs [58,59], space and time like (quasi) crystals, pattern forming and nonlinear gravitational wave structures and applications in modern cosmology see [60][61][62][63][64] and references therein. We elaborate on new classes of generic off-diagonal stationary and cosmological solutions with entropic geometric flows which for self similar Ricci soliton configurations result in equations considered in Verlinde works [6,7] and a covariant generalization due to Hossenfelder [65], see also critics Refs. [66,67].
Three lines of evidence motivate this article and the partner letter [32]. First, we use our former results and nonholonomic geometric methods [33][34][35][36][37][38][39][40][41] that generalized/modified relativistic flow equations and Einstein equations in GR and MGTs can be derived as systems of PDEs for modified Ricci solitons for respective nonholonomic modifications of G. Perelman's F-and W-entropy functionals. On modified gravity and applications in modern cosmology and astrophysics, see reviews [63,[68][69][70]. Second, a number of recent works invoke ideas on origin of gravity as an emergent effect of the entropic force, entanglement etc., see [4,6,7,[16][17][18][19][20][21][22][23][24][25][26]. We argue that this can be grounded and explained, at a deeper level, through modifications of the Poincaré-Thurston conjecture on geometric flows, when the F-and Wfunctionals are generalized for metrics and generalized connections on Lorentz manifolds and/or certain supersymmetric/noncommutative/fractional/stochastic generalizations. In our approach, the spacetime evolution and gravity are treated via geometric entropy values which allows to formulate respective statistical thermodynamics models. Third, new advanced methods for constructing exact solutions in MGTs and GR allows us to proceed directly toward definition of gravitational entropy and thermodynamic values making no use of holography, area-entropy relation, CFT duality etc. Due to the competition between area and volume law of generalized W-entropy, we can characterize thermodynamically new classes of BH and cosmological solutions with quasiperiodic structure, locally anisotropic inflation and acceler-ating scenarios, exhibiting memory effects in the form of an entropy displacement caused by matter etc.
This article is organized as follows: In Sect. 2, we provide an introduction to the geometry of double, 2 + 2 and 3 + 1 dimensional, spacetime fibrations defining elastic and quasiperiodic configurations both for gravitational and (effective) matter fields. Main concepts and most important results on nonlinear connection geometry on nonholonomic Lorentz manifolds, hypersurface geometric objects, and modified emergent/elastic gravity theories are outlined. There are studied the geometry of distributions, and respective Lagrange densities and geometric evolution or dynamical fields defining elastic and quasiperiodic structures. Section 3 is devoted to the theory of geometric flows and modified entropic gravity. We postulate canonical nonholonomic deformations of Perelman's F-and W-functionals encoding geometric flow evolution scenarios of entropic spacetimes with quasiperiodic structure. Such values are defined in relativistic 4-d form and for 3-d hypersurface projections. There are derived respective (generalized R. Hamilton) geometric flow equations for entropic quasiperiodic flows. The concept of nonholonomic Ricci solitons as self-similar configurations is elaborated and related modifications of the Einstein equations are analyzed. We speculate on connection between relativistic generalizations of the Poincaré-Thurston conjecture for Ricci flows and geometric proofs of E. Verlinde's conjecture.
In Sect. 4, we develop and apply the anholonomic frame deformation method, AFDM, [32,38,40,58,63,64] in order to prove general decoupling properties and integrability of nonholonomic geometric flow and Ricci soliton equations encoding elastic and quasiperiodic spacetime and (effective) matter fields properties. Such solutions are described by generic off-diagonal metrics, and generalized connections, depending on all spacetime coordinates and temperature like parameters via general classes of generating functions and (effective) sources of entropic gravity and matter fields.
Then, in Sect. 5, we consider the AFDM for constructing cosmological solutions for entropic quasiperiodic flow and MGTs. We emphasize that there are certain nonlinear symmetries relate possible classes of generating functions and (effective) sources all encoding entropic, quasiperiodic, pattern forming, space and time quasicrystal, solitonic and other type structures. There are studied cosmological configurations generated by entropic quasiperiodic sources, nonstationary generating functions, cosmological metrics evolving in (off-) diagonal elastic and/or quasiperiodic media.
Finally, we conclude our work and discuss certain perspectives of the theory of G. Perelman and E. Verlinde entropic geometric flow and emergent gravity theories in Sect. 6.

Spacetime + and 3 + 1 fibrations with elastic and quasiperiodic structures
In this section, we summarize necessary results on the geometry of Lorentz manifolds enabled with nonholonomic (i.e. non-integrable, equivalently, anholonomic) distributions defining double 2 + 2 and 3 + 1 fibrations. There are developed nonholonomic geometric methods which are important for elaborating theories of relativistic Ricci flows and possible applications in modern cosmology and astrophysics, see details in [38][39][40].
As explicit examples, we shall consider nonholonomic distributions modelling elastic and/or quaisperiodic space and time structures (for instance, quasicrystal or solitonic like configurations) [6,7,59,61,62,64,65]. It should be noted that the 2 + 2 nonholonomic splitting is important for proofs of general decoupling and integration properties of the relativistic geometric and entropic flow evolution, nonholonomic Ricci soliton and (entropic modified) Einstein equations, see Sect. 4. Additional 3 + 1 decompositions adapted to 2 + 2 splitting will be used for defining and computing entropic and thermodynamic like values for various classes of solutions of physically important systems of nonlinear partial differential equations, PDEs, see Sect. 5.

Nonlinear connections with 2 + 2 splitting of Lorentz manifolds
Let us consider a four dimensional, 4-d, Lorentzian manifold V, dim V = 4, with local pseudo-Euclidean signature (+ + +−) for a metric field g = (hg, vg). The conventional horizontal, h, and vertical, v, nonholonomic decomposition is defined by a nonlinear connection, N-connection, structure N. Such a geometric object can be always introduced as a Whitney sum where T V is the tangent bundle on V. The concept of nonholonomic manifold is used for a manifold enabled with a nonholonomic distribution. In this work, this refers to a Lorentz spacetime V := (V, N) enabled with a N-connection structure of type (1). In local coordinates, where N a i are N-connection coefficients. 2 Any set {N a i } defines subclasses of N-linear (co) frames which allows N- 2 We can parameterize the local coordinates in the form u μ = (x i , y a ), (in brief, u = (x, y)), where indices respectively take values i, j, . . . = 1, 2 and a, b, . . . = 3, 4, considering that u 4 = y 4 = t is the time like coordinate. The Einstein convention on summation on "up-low" repeating indices will be applied if contrary will not be stated for some special cases. We use boldface symbols for spaces and geometric objects adapted to a N-connection splitting. adapted diadic decompositions of geometric and physical objects. 3 On any nonholonomic manifold V, we can consider covariant derivatives determined by affine (linear) connections which are, or not, adapted to a N-connection structure. A distinguished connection, d-connection, is a linear connection D = (hD, vD) which preserves under parallel transport a h-v-decomposition (1). 4 For any d-connection D, we can define and compute in standard form the d-torsion, T, the nonmetricity, Q, and the d-curvature, R, tensors where X and Y are vector fields (i.e. d-vectors) on T V. 5 Any metric tensor g = (hg, vg), on a nonholonomic V can be written as a distinguished tensor, d-tensor (d-metric), with respective splitting into h-and v-indices, where the nonholonomic dual frame structure e α is chosen in a form when the matrix of metric coefficients g αβ is considered in diagonal form for g α := g αα , g i := g ii and g a := g aa . With respect to a dual local coordinate basis du α the same metric field is expressed where g αβ = g i j + N a i N b j g ab N e j g ae N e i g be g ab .
3 Such N-adapted local bases, e ν = (e i , e a ), and cobases, e μ = (e i , e a ), are defined by formulas . Holonomic (integrable) configurations are obtained if and only if W γ αβ = 0. 4 In general, a linear connection D is not adapted to a prescribed Nconnection structure, i.e. it is not a d-connection. In such a case, one should be not used a boldface symbols for respective geometric objects determined by D. 5 We can compute in N-adapted form the coefficients of any d- Using frame transforms (in general, not N-adapted), we can transform any metric into a d-metric (2) an off-diagonal form with N-coefficients. For nontrivial anholonomy coefficients, such a metric is generic off-diagonal.
For our geometric constructions, there are two important linear connections determined by the same metric structure: Here we note that if we prescribe a N-connection structure N, we can define a canonical d-connection D and compute certain nontrivial torsion coefficients hvT completely defined by certain off-diagonal coefficients containing N a i (u) in (3) and/or nontrivial anholonomy coefficients W γ αβ , see footnote 3. Of course, we can introduce an infinite number of metric compatible d-connections but not all such connections allow to decouple physically important systems of nonlinear PDEs (in our case, for nonholonomic geometric flows and modified, or Einstein, gravity). A D (4) allows us to prove the decoupling and integration properties of such equations in Sect. 4.
The LC-connection ∇ (4) can be defined uniquely by a metric g without any N-connection structure but ∇ can be distorted always to a necessary type d-connection allowing a general decoupling and integrability of certain important physically important systems of nonlinear PDEs. In our previous works [38][39][40] (see there necessary geometric details and [59,61,62,64] for applications of the AFDM), we used a "hat" symbol (like D) for the canonical d-connection in (4). In this paper, we shall work only with ∇ and D = D and omit "hats" on respective geometric objects. We note that all constructions performed for ∇ and D are related by a distortion relation, D[g, N] = ∇[g, N] + Z[g, N], where Z is the distortion tensor determined in standard algebraic form by the torsion tensor T; all values are completely defined by the metric tensor g adapted to N. 6 The Ricci tensors of D and ∇ are defined and computed in standard forms for different linear connection structures but defined by the same metric tensor by contracting respective indices. We denote them, respectively, Ric = {R αβ := R γ αβγ } and Ric = {R αβ := R γ αβγ }. Any (pseudo) Riemannian geometry can be equivalently described by both geometric data (g, ∇) and (g, N, D), when the canonical distortion relations R = ∇ R + ∇ Z and Ric = Ric + Zic, with respective distortion d-tensors ∇ Z and Zic, are computed for the canonical distortion relations D = ∇ + Z, see details 6 The values hT and vT are respective torsion components which vanish on conventional h-and v-subspaces, but there are nontrivial components hvT defined by certain anholonomy (equivalently, nonholonomic/nonintegrable) relations. Such a d-torsion is induced by nonholonomic configurations.
Using N-adapted coefficients of the canonical Ricci dtensor, we can compute the scalar of canonical d-curvature, s R := g αβ R αβ = g i j R i j + g ab R ab . This geometric object is different from the scalar curvature of the LC-connection, R := g αβ R αβ .
Using ∇, the Einstein equations in GR are written in standard form, In these formulas, m T αβ is the energy-momentum tensor of matter fields A ϕ determined by a general Lagrangian m L(g,∇, A ϕ), where is the gravitational coupling constant for GR. 7 We can define nonholonomic gravitational field equations using the Ricci d-tensor (5) for a canonical d-connection D Such equations are equivalent to (6) if there are imposed additional nonholonomic constraints, or found some smooth limits, for extracting LC-configurations, D | T =0 = ∇, for instance, of type In (7), a matter fields source ϒ μν can be constructed using a N-adapted variational calculus for m L(g, D, A ϕ), when In such formulas, we consider m T = g μν m T μν for We note that any (pseudo) Riemannian geometry and gravity theory, and various metric-affine modifications (for instance, F(R)-modified theories [63,68]), can be formulated equivalently using geometric data (g, ∇) and/or (g, D). There is an important motivation to use nonholonomic variables of type (g, D) because that they allow to decouple and integrate in general form various modified and standard Einstein equation. Such solutions can be with generic offdiagonal metrics and coefficients depending on all spacetime coordinates [55][56][57]. A recent review of the so-called anholonomic frame deformation method, AFDM, of constructing 7 We use abstract left labels A and m in order to distinguish the values from similar notations of pure geometric objects, for instance, T γ αβ . exact solutions in GR and MGTs, geometric flow theory, and applications in modern cosmology and astrophysics can be found in [63]. In this work, we shall develop the AFDM for constructing exact solutions in entropic geometric flow and gravity theories.
2.2 Nonholonomic 3 + 1 splitting adapted to 2 + 2 structures We outline some basic concepts on the geometry of 3+1 foliations of a nonholonomic Lorentzian manifold (V, g, N) of signature (+ + +−) into a family of non-intersecting space like 3-d hypersurfaces t parameterized by a "time function", t (u α ), stated as a scalar field as described as follows. Such spacetime decompositions are useful for elaborating various thermodynamic, locally anisotropic kinetic [71] and geometric evolution or hydrodynamic flow models [38] when a conventional splitting into time and space like coordinates is necessary. This allows definition of physical important values (for instance, entropy, effective energy etc.) and deriving fundamental geometric evolution equations. In our approach, we generalize the well-known geometric 3 + 1 formalism in GR (see, for instance, [72]) to the case of nonholonomic manifolds [38][39][40]. For a 3-d manifold , we consider an one-to-one image to a hypersurface = ζ( ) ⊂ V constructed as an homeomorphism with both continuous maps ζ and ζ −1 , when does not intersect itself. Left "up" or "low" labels by a vertical bar " " will be used in order to emphasize that certain geometric objects refer to 3-d manifolds/hypersurfaces. Such a 3-d space is supposed to be locally defined as a set of points for which a scalar field t on V is constant (for instance, i.e. t ( p) = 0 for any point p ∈ ). It is assumed also that t spans the real line R and that any is a connected submanifold of V with the topology of R 3 . 8 It should be noted that any 2 + 2 splitting by a nonholonomic distribution N (1) induces a N-connection structure for a hypersurface , i.e. an induced N-connections N : Using the coefficients of such an induced N-connection, any induced 3-metric tensor q can be 8 We can label local coordinates for a 3 + 1 splitting in u α = (x`ı , t), where indices α, β, . . . = 1, 2, 3, 4 andì,j, . . . = 1, 2, 3 are related to a 2 + 2 splitting as in previous subsection (in brief, we shall write u = (ȗ, t)). The continuous maps ζ can be parameterized to "carry along" curves/vectors in to curves/vectors in V, for ζ : (x`ı ) −→ (x`ı , 0). This way, it is possible to define and relate respective local bases ∂`ı := ∂/∂ x`ı ∈ T ( ) and ∂ α := ∂/∂u α ∈ T V. The coefficients of 3-vectors and 4-vectors are expressed correspondingly, a = a`ı ∂`ı and a = a α ∂ α (we shall use also capital letters, for instance, A = A`ı ∂`ı and A = A α ∂ α ). Similar formulas are considered for dual forms to vectors, 1-forms, when the dual bases dx`ı ∈ T * ( ) and du α ∈ T * V. The 1-forms will be parameterized for respective 3 and 4 dimensions, Ã = A`ı dx`ı andÃ = A α du α . We shall omit the left/up label by a tilde ∼ (writing A and A) if that will not result in ambiguities. written in N-adapted frames as a d-tensor (d-metric) in the form choosing common frame and coordinate systems for ⊂ V. We can extend naturally such a 3-d metric q to a 4-d d-metric g (2) reparameterized in a form adapted both to 2 + 2 and 3 + 1 nonholonomic splitting, (10), the lapse functionN (u) > 0 is defined as a positive scalar field which ensues that the d-vector n is a unite one. An "inverse hat" symbol is used in order to distinguish such a symbol from N is used traditionally in literature on GR [72]. Here we note that in another turn, the symbol N a i is used traditionally for the N-connection and this also motivates a new symbolN .
We note that for any quadratic line element ds 2 = g αβ du α du β of a metric tensor g there are such frame transforms to parameterizations wheng`ı`j = q`ı`j = g αβ e αì e βj is just the induced metric on t . In result, the determinants of 4-d and 3-d metrics are computed √ |g| =N |g| =N √ |q|.
Using certain coordinates (x`ı , t) being N-adapted on respective hypersurfaces, the time partial derivatives are computed £ t q = ∂ t q = q * and the spacial derivatives are computed q ,`ı := e ὰ ı q ,α . There are two types of induced linear connections completely determined by an induced 3-d hypersurface metric q, Such formulas are related to 4-d similar ones (4). Both linear connections, ∇ and D, are subjected also to a distortion relation For 3-d configurations, we can compute the N-adapted coefficient formulas for nonholonomically induced torsion structure T = { T`ı`jk}, determined by D, and for the Riemannian tensors R = { R`ı`jkl } and R = { R`ı`jkl }, determined respectively by ∇ and D. Using 3-d subsects of coefficient formulas, we can compute respective N-adapted hypersurface coefficients of the Ricci d-tensor, R`jk, and the Einstein d-tensor, E`jk. Contracting indices, we obtain the Gaussian curvature, R = q`jk R`jk, and the Gaussian canonical curvature, s R = q`jk R`jk, of ( , q, N). It should be noted that all this types of N-adapted and not N-adapted geometric objects can be defined in abstract form which do not depend on the type of embedding of a nonholonomic 3-d manifold ( , q, N) into a 4-d one (V, g, N).

Quasiperiodic space and time QC configurations
Let us consider two examples of space and time quasiperiodic structures defined in a curved spacetime following our works on quasicrystal, QC, models in modern cosmology [61,63,64] (alternative models are studied in [73]). Our approach was elaborated following Wilczek and co-authors works in condensed matter physics [74][75][76][77]. As a toy model, we consider one dimensional, 1-d, time quasicrystals, TQCs with time structure equations generalizing those introduced in [77]). Then we introduce some important formulas on three dimensional, 3-d, QC structures -in general, such configurations are called space-time quasicrystal structures, STQC, and studied in [64]. In this work we use a different system of notation for partial derivatives when, for instance, ∂q/∂

1-d relativistic time QC structures
We consider a scalar field ς(x i , y a ) on a space-time (V, g, N) and respective Lagrange densitý In this formula,V (ς ) is a nonlinear potential and e α are Nadapted partial derivatives. Corresponding N-adapted variational motion equations are The field ς defines a 1-d time QC structure, 1-TQC, if it is a solution of these motion equations. 9

3-d QC structures on curved spaces
QC structures and analogous dynamic phase field crystal models can be elaborated as flow evolution theories on real parameter τ (in next section, this parameter will be identified with a geometric flows one). Such a QC structure can be 9 For non-relativistic limits with g αβ = [1, 1, 1, −1] and ς → ς(t), , which leads to an effective energy E = defined by a generating function q = q(x i , y 3 , τ ) subjected to the condition that it is a solution of an evolution equation with conserved dynamics, Such evolution is considered on 3-d spacelike hypersurface t when the canonically nonholonomically deformed hypersurface Laplace operator indices rung valuesì,j, . . . 1, 2, 3. This operator is a distortion of := ( ∇) 2 constructed in 3-d Riemannian geometry, see previous subsection. The functional F in (13) is characterized by an effective free energy where q = det |q`ı`j |, δy 3 = e 3 and the operators and Q are defined and explained in [63,64]. Such nonlinear interactions are stabilized by the cubic term with Q and the second order resonant interactions are varied by setting observable values of such constants (they are different for cosmological models, in astrophysics or condensed matter physics). The average value b is conserved for any fixed time variable t and/or evolution parameter τ 0 . We can fix b |τ =τ 0 = 0 when other values are accommodated by redefining values and Q.

Distributions defining spacetime elastic configurations
In letter [32], we shown that models of entropic gravity can be derived from nonholonomic modifications of the Wfunctional when the (modified) Einstein equations are equivalent to certain nonholonomic Ricci soliton equations. Here, we shall study the conditions when entropic elastic scenarios can be modelled as nonholonomic Ricci solitons in Sect. 3.3. We shall consider certain examples of nonholonomic distributions and related Lagrange densities on a Lorentz manifold V which are used in entropic gravity theories [6,7,32,65]. Using such geometric constructions, we shall elaborate in next section on elastic flow evolution models and their selfsimilar nonholonomic Ricci soliton configurations. There are three important values: These geometric/physical objects are determined by a conventional displacement vector field u α , cosmological constant and some constants α, β, γ ; there are used short hand notations: u := √ |u α u α |, ε = ε β β , and n α := u α /u.
On V, there are considered nonholonomic distributions for corresponding total, effective gravitational, usual matter, interaction and kinetic terms of Lagrangians postulated in the form In these formulas, the Plank gravitational mass is denoted M P and the gravitational Lagrangian g L is taken as in modified gravity [63,[68][69][70]. We can fix z = 1 if we search for compatibility with [65], or z = 2 if we search for a limit to the standard de Sitter space solution [66,67] (as we use in [32]).
To model STQC structures in entropic gravity and related geometric flow theories we can consider that the displacement vector field u α [ς, b] is a functional of functions ς, b subjected to certain conditions of type (12) and/or (13) [in principle, we can consider functionals for pattern forming, nonlinear wave soliton structures, fractional and diffusion processes etc.]. The energy-momentum tensors considered in above formulas and/or derived from respective Lagrangians in (14) and computed using variations on g μν similarly to m T μν (9) (in N-adapted form, details of such computations are provided in [39,58,59,63]). For the full system, the effective energy-momentum tensor is computed We can model "pure" elastic spacetime modifications of the Einstein gravity if we fix F( s R) = s R and consider restrictions to the Levi-Civita connection D = ∇. For such conditions, we obtain respective formulas for int T μν and χ T μν which are similar to formulas (10)- (13) in [66]. 10 In this work, the generalized (effective) source for MGT (7) splits into four components, where is determined in standard form by the Newton gravitational constant G. We need additional terms and parame- 10 We use a system of notations which is similar (but without "hats") to [35,[38][39][40]63]; such notations are different from those used in [65][66][67].
terizations in order to describe structure formation in modern cosmology and to model dark energy and dark matter properties.

Relativistic geometric flows and modified entropic gravity
Grigory Perelman's proof of the Poincaré conjecture [45] on geometric flow evolution of 3-d Riemannian metrics [42][43][44] provided fundamental results in geometric analysis and topology. There were also studied possible applications in modern mathematical and particle physics. We cite [46][47][48] for reviews of rigorous mathematical results. For early applications, we refer to D. Friedan works [49][50][51] (he considered geometric evolution equations related to renorm group equations before the Hamilton-Poincaré theory was elaborated). Further developments and applications were performed in [52][53][54] and a series of works [33][34][35][36][37][38][39][40][41], see also references therein. In those works on theories of nonholonomic/noncommutative/supersymmetrics, fractional, diffusion etc. geometric flows, there were studied statistical and thermodynamic evolution models derived from certain Lyapunov type functionals. Such F-and W-entropy functionals are called in literature the Perelman functionals. The Wentropy has properties of "minus entropy" of statistical thermodynamics systems. In [32], we elaborate on the idea that such a W-entropy can be considered for formulating Verlinde type entropic gravity theories [6,7]. We study self-similar configurations of nonholonomic geometric flows resulting in entropic Ricci solitons (see Sect. 3.3). We note that Perelman suggested in his first preprint [45] that the geometric flow theory may have certain implications in black hole physics and string theory. Nevertheless, the original theory of Ricci flows was formulated in a nonrelativistic form. To consider further generalizations and applications in modern physics and cosmology we elaborated on relativistic models of geometric flow theories [38][39][40]. Such constructions can be re-defined for nonholonomic configurations modeling elastic and quasiperiodic spacetime structure as in Sect. 2.4 and allows rigorous geometric motivations for emergent entropic theories of type [6,7,32,65].
The goal of this section is to study generalizations of the hypersurface 3-d and relativistic 4-d F-and W-functionals and elaborate on respective geometric evolution scenarios supporting the Verlinde entropic gravity conjecture [6,7]. Such constructions can be considered in the framework of s a modified relativistic variant of the Poincaré-Thurston conjecture which was proven only for certain classes of Riemannian and Kähler metrics, see details in [45][46][47][48]. For relativistic configurations, we can only elaborate on geometric evolution of certain 3-d hypersurface configurations subjected to the conditions that such 3-metrics can be extended to cer-tain classes of 4-d metric and (non) linear connection structures derived as exact/parametric solutions of certain nonholonomic/entropic geometric flow equations. There will be considered also generalizations of the Hamilton equations for the entropic flow theory. The conditions for generating entropic modified Einstein equations as nonholonomic Ricci solitons will be also analysed. We emphasize that in the main part of this article there are studied relativistic geometric flow models with a temperature like evolution parameter.

Modified spacetime and hypersurface Perelman's functionals
Let us consider families of nonholonomic 4-d manifolds V(τ ) parameterized by a positive parameter τ, 0 ≤ τ ≤ τ 0 (it can be considered as a temperature like parameter) and enabled with a double nonholonomic 2 + 2 and 3 + 1 splitting [38][39][40]. Such manifolds are determined by respective families of metrics g(τ ) = g(τ, u) and N-connections N(τ ) = N(τ, u) (we shall write only the parametric dependence if that will not result in ambiguities) for which canonically corresponding d-connection structures can be constructed D(τ ) = D(τ, u). We also suppose that on V(τ ) there are defined corresponding families of Lagrange densities g L(τ ), for gravitational fields in a MGT or GR, and tot L(τ ), as total Lagrangians for effective and matter fields (14). For a double 2 + 2 and 3 + 1 splitting, we can consider local coordinates labeled as

Generalized Perelman functionals for entropic geometric flows and MGTs
For this class of theories, we postulate the modified Perelman's functionals in the form The condition is imposed on the normalizing function f (τ, u). For topological considerations, such a normalisation is not important. Nevertheless, it imposes certain nonholonomic constraints on geo-metric objects which do not allow to solve derived geometric flow evolution equations in explicit form. For applications to entropic gravity and associated thermodynamic models, we can consider f as an undetermined scalar function which can be related to possible conformal transforms or re-parameterizations. In result, we can prove certain general decoupling and integration properties of corresponding systems of nonlinear PDEs. Fixing a class of solutions, we can chose such integration functions and constants which reproduce/predict certain experimental and/or observational data.
Corresponding values of f depend on systems of reference and coordinates.
Let us explain and motivate the difference of (16) and (17), introduced in the first partner work [32], from the original Grisha Perelman F-and W-functionals [45] postulated for the Ricci flows of 3-d Riemannian metrics, see details in monographs [46][47][48]. In this work, we study geometric entropic flows of canonical geometric data (g(τ ), N(τ ), D(τ )) for nonholonomic Lorentz manifolds and various generalizations for MGTs following the program elaborated in [35,[38][39][40]63], where possible connections to emergent gravity were not analyzed. In formulas (16) and (17), we consider the gravitational Lagrangian g L = F( s R) as a functional of the scalar curvature for D, or g L = R[∇] for considering as particular cases models of geometric evolution of exact solutions in GR. The key difference from previous works is that in such relativistic functionals the term tot L is introduced, which is responsible for geometric evolution of configurations with elasticity and quasiperiodicity. Those functionals can be generalized on a temperature like parameter τ and used as certain alternative geometric functionals, for instance, for W-entropy. Nevertheless, only nonholonomic elastic quasiperiodic functionals of type (16) and (17) result for self-similar configurations (see next subsections) in entropic gravity equations of Verlinde type [6,7,32,65] and/or with quasiperiodic structure [61,63,64].
In this and partner [32] papers, we work with generalized geometric flow and entropy functionals determined by F( s R) + tot L and D, respectively, instead of the Riemannian values R and ∇ used in the former mathematical works. In our nonholonomic approaches, above F-and W-functionals characterize relativistic thermodynamic models with analogous nonlinear hydrodynamic flows of families of entropic values, metrics and generalized connections, encoding interactions of gravitational and matter fields as it is motivated in [38][39][40]. In general, it is possible to work with any class of normalizing functions f (τ, u) which can be redefined in order to include geometric and matter Lagrange terms and certain constant values and parameters. In many cases, such a function is chosen in a non-explicit form. This allows us to study non-normalized geometric flows but with nonholonomic constraints. For such conditions, there found various general decoupling and integration properties of respective physically important systems of nonlinear PDEs. In result, generic off-diagonal solutions can be constructed in explicit form as in [32,35,[38][39][40]63], but with entropic and quasiperiodic modifications. The existence of such solutions validates our nonholonomic geometric flow entropic approach, involving metrics with pseudo-Euclidean signature even analogs of the Poincaré-Thurston conjecture have not been formulated and proven for the Lorentzian spacetimes. Nevertheless, explicit constructions of exact solutions with elastic and quasiperiodic gravitational metrics and effective matter sources, which will be provided in Sects. 4-6, support E. Velinde conjecture on entropic gravity which in our works is proven for modified Poincaré functionals.

Nonholonomic 3-d space like hypersurface F-and W-functionals
We can redefine and compute relativistic entropies (16) and (17) for any 3 + 1 splitting with 3-d closed hypersurface fibrations t as we described above in Sect. 2.2. Let us denote by D = D | t the canonical d-connection D defined on a 3-d hypersurface t , when all values depend on a temperature like parameter τ (τ ) with possible scale re-definitions for another parameter τ etc. We define also in a family of d-metrics (10), the Perelman's functionals parameterized in N-adapted form are constructed in the form: These functionals are derived respectively from the previous 4-d elastic ones when the values F( s R) and tot L are computed as projections on a 3-d hypersurface for a redefined normalization function f. Using frame/coordinate transform and re-definition of the temperature like parameter, we can always chose a necessary type scaling function f which satisfies normalization conditions t μ |q`ı`j |dx 3 = 1 for The functionals (18) and (19) transform into standard Perelman functionals [45] for 3-d Riemannian metrics on t if D → ∇, F( s R) = s R and tot L = 0. In order to describe possible contributions on 3-d hypersurfaces of spacetime elasticity and quasiperiodic structure in entropic gravity, it is necessary to analyze physical effects of such nonholonomic deformations.

Geometric flow equations for modified gravitational and matter fields
Applying a variational procedure for a corresponding Ffunctional for geometric flows of 3-d Riemannian metrics, Perelman [45] provided a proof for R. Hamilton's equations [42][43][44]. For self-consistent configurations with a fixed flow parameter τ 0 , one obtains 3-d Ricci soliton equations which are equivalent to the vacuum Einstein equations for ∇ with an effective cosmological constant. In similar forms to rigorous mathematical proofs in [45][46][47][48] but elaborating on N-adapted variational procedures, for instance, for the functional F(τ ) (16) with a canonical D used instead of ∇ (see details in [32,[38][39][40]), we obtain a system of nonlinear PDEs generalizing the R. Hamilton equations for entropic and quasiperiodic geometric flow evolution determined by canonical data In these formulas, (τ ) = D α (τ )D α (τ ) and tot ϒ αβ (τ ) is chosen for geometric flows of (effective) sources of entropic gravity (15) (if we fix any τ = τ 0 ). We note that the dependence on a flow parameter τ for such (effective) matter sources is determined by certain evolutions of g(τ ) and D(τ ).
In such theories, we do not consider nonholonomic deformations and evolution of classical matter fields. For instance, we do not consider geometric flow evolution equations for the electromagnetic potential A α (τ ) with evolution terms of type ∂ τ A α even such theories were studied in our previous works [33][34][35][36][37][38][39][40][41], see references therein. The conditions R ia = 0 and R ai = 0 for the Ricci tensor } are necessary if we want to keep the metric g(τ ) to be symmetric under nonholonomic Ricci flow evolution determined by (20). Geometric flow evolution and nonholonomic gravity of theories with nonsymmetric metrics were studied in [56], see references therein. In principle, we can work with any type normalization function f which allows a general decoupling and integration of such systems of nonlinear PDEs. Such a normalization depends on frame and coordinate transforms and may encode (effective) cosmological constants, matter sources etc. We note that similar variational and/or geometric methods allows to derive from W(τ ) (17) certain types nonlinear evolution equations which are equivalent to (20). It is more difficult to solve explicitly such PDEs but a Wfunctional allows to elaborate directly on certain classes of thermodynamic models, see Sect. 5.

Entropic gravity and gravitational field equations as Ricci solitons
For self-similar point τ = τ 0 configurations when ∂ τ g μν = 0, with a corresponding choice of the normalizing geometric flow function f, the equations (20) transform into relativistic nonholonomic Ricci soliton equations which are equivalent to (modified) Einstein equations in (MGT) GR for corresponding definitions of tot ϒ αβ . A class of MGTs and GR can be formulated as geometric theories of entropic elastic origin which is similar to the idea of emergent gravity put forward by Verlinde [6,7], i.e. in the form (7) with (effective) entropic and quasiperiodic source tot ϒ αβ (15). We conclude that an emergent gravity model in the Verlinde sense [6,7,65], can be constructed for Lagrange distributions (14) and respective sources (15) introduced as generating data for the nonholonomic Hamilton equations (20) and respective relativistic Ricci solitons. Such geometric flow evolution theories and their spacetime elastic, quasiperiodic and thermodynamic properties are determined by the generalized W-entropy (17).

Decoupling and integrability of entropic flow equations
In this section, we prove that the system of nonlinear PDEs (20) describing spacetime elastic and quasiperiodic flows and entropic gravity theories can be formally integrated in very general forms for generic off-diagonal metrics and canonical d-connections (in particular, for LC-configurations). The coefficients of geometric objects for such solutions depend on all spacetime coordinates via generating and integration functions and (effective) matter sources. The anholonomic frame deformation method, AFDM, for constructing exact solutions in MGTs and GR is developed for generating new classes of solutions encoding entropic quasiperiodic modifications in g(τ ) (2), D(τ ) (4), and tot ϒ αβ (τ ) (15). For similar details and mathematical proofs, we refer readers to our previous works [55][56][57][58][59][60]63,64], on exact solutions in MGTs, and [33][34][35][36][37][38][39][40][41], for solutions with nonholonomic Ricci flows, and citations therein.

Geometric flows with parametric modified Einstein equations
Introducing effective sources, entropic geometric flow equations can written as modified Einstein equations with dependence on a temperature like parameter τ . We show that such systems of nonlinear PDEs can be decoupled in general forms.

Entropic quasiperiodic flow modifications of gravitational field equations
Using nonholonomic frame transforms and tetradic (vierbein) fields, we introduce effective sources which in Nadapted form are parameterized Such families of vielbein transforms e μ μ (τ ) = e μ μ (τ, u γ ) and their dual e ν ν (τ, u γ ), when e μ = e μ μ du μ can be chosen for any frame/coordinate transforms of a N-splitting structure (1). In result, the system of nonholonomic entropic R. Hamilton equations (20) can be written in the form (7) but with geometric objects depending additionally on a temperature like parameter τ and for effective source (22), We note that such geometric evolution equations are for an undetermined normalization function f (τ ) = f (τ, (τ, u γ ) which can be defined explicitly for respective classes of exact or parametric solutions. For self-similar point τ = τ 0 configurations with ∂ τ g μν (τ 0 ) = 0, this system of nonlinear PDEs transforms into the nonholonomic entropic Ricci soliton equations (21).
The AFDM results in more simple and explicit (still very general classes) of solutions if we work with nonholonomic configurations possessing at least one Killing symmetry, for instance, on ∂ 4 = ∂ t for stationary solution or on ∂ 3 = ∂ ϕ , locally anisotropic solutions. 11 We shall use brief notations of partial derivatives ∂ α q = ∂q/∂u α when a function q(x k , y a ), For respective Killing symmetries, the effective sources (τ, x, y) in (22) can be parameterized Considering as typical examples two types of a Killing space symmetry or time like Killing symmetry for effective generating sources, we shall construct and study properties of two general classes of exact solutions (the first one will be for stationary configurations which may contain BH solutions and the second one will be for cosmological type solutions).

Nontrivial Ricci d-tensors and decoupling of entropic flow equations
In this subsection, we outline the key steps for proofs of general decoupling and integrability of (modified) Einstein equations with effective sources (26).

Off-diagonal metric ansatz, (non) holonomic variables, and ODEs and PDEs
Let us summarize in Table 1 below the data on nonholonomic 3 + 1 and 2 + 2 variables and corresponding ansatz which allows to transform geometric and entropic flow equations and, a nonholonomic Ricci solitons, gravitational field equations in entropic MGTs and GR into respective systems of nonlinear ordinary differential equations, ODEs, and partial differential equations, PDEs. All formulas will be proven in next subsections. We model a nonholonomic deformation 11 In principle, we can construct for (23) certain classes of exact and parametric off-diagonal solutions generically depending on all spacetime coordinates (x k , y a ) but that would result in hundreds of pages with a cumbersome formulas for respective geometric techniques, see [55][56][57][58][59] and references therein.
with η-polarization functions,g → g(τ ), of a 'prime' metric, g, into a family 'target' d-metrics g(τ ) (2), if where the target N-elongated basis is determined by N a η a (τ, x k , y b ) and η a i (τ ) = η a i (τ, x k , y b ) are called respectively geometric/entropic flow or gravitational polarization functions, or η-polarizations. Any g(τ ) is subjected to the condition that it defines a solution of modified Einstein equations resulting in entropic quasiperiodic geometric flows and/or via nonholonomic deformations. A general prime metric in a coordinate parametrization is of typeg = g αβ (x i , y a )du α ⊗du β , which can be also represented equivalently in N-adapted form Such a d-metric can be, or not, a solution of some gravitational field equations in a MGT or GR but it nonholonomic deformations to a target metric (27) are subjected to the condition to define an exact or parametric solutions of certain entropic flow evolutions equations. In our works, we are interested usually in two physically important cases wheng (28) defines a BH solution (for instance, a vacuum Kerr, or Schwarzschild, Kerr-(anti) de Sitter metric), or a Friedman-Lemaître-Robertson-Walker (FLRW) type metric, or any Bianchi anisotropic metrics. For diagonalizable prime metrics (the off-diagonal structure of the Kerr metric is determined by rotation frames and coordinates), we can always find a coordinate system when N b i = 0. To avoid nonholonomic deformations with singular coordinates is convenient to construct exact solutions with nontrivial functions η α = (η i , η a ), η a i , and nonzero coefficientsN b i (u). We have to consider necessary type frame/coordinate transforms. For a d-metric (27), we can analyze the conditions of existence and geometric/physical properties of some target and/or prime solutions, for instance, when η α → 1 and N a i →N a i . The values η α = 1 and/or N a i = 0 can be imposed as some special nonholonomic constraints. 13 In brief, we shall denote certain nonholonomic 12 We do not consider summation on repeating indices if they are not written as contraction of "up-low" ones. 13 We can consider flow evolution of a physical important target metric g (27) with generic off-diagonal terms as an "almost" BH, or FLRW cosmological, like metric. Such parametric solutions are constructed for small nonholonomic deformations on some constant parameters η α = (η i , η a ), η a i , for 0 ≤ ε α , ε b i 1, when η i [N-adapt. fr.] τ, r, θ), g a (τ, r, θ, ϕ), or g i (τ, r, θ), g a (τ, r, θ, t), d-metrics r, θ, t),  Table 1 outlines the key steps for developing the AFDM to theories of entropic quasiperiodic geometric flows. In this work, the formulas depend on a flow temperature like parameter τ and the constructions are for effective matter sources encoding entropic nonholonomic flows and deformations.

Cosmological Ricci d-tensors, LC-conditions, and nonlinear symmetries
For locally anisotropic cosmological configurations, we can consider geometric data when coefficients of the geometric objects do not depend on a space like y 3 with respect to certain classes of N-adapted frames. Using d-metric data (25) with ω = 1 and a source h (τ, x i ), ( τ, x i , t) (26), we write the entropic flow modified Einstein equations (23) in , x k , y b ). Parametric ε-decompositions can be performed in a self-consistent form by omitting quadratic and higher terms after a class of solutions have been found for some evolution or nonholonomic deformation data (η α , η a i ). For certain subclasses of solutions, we can consider that ε i , ε a , ε a i ∼ ε, when only one small parameter is considered for all coefficients of nonholonomic deformations. We can work with mixed types of solutions and model only small diagonal deformations ε i , ε a , ∼ ε of metrics, for some general η a i . Alternatively, we consider nontrivial η α but ε a i ∼ ε. the form 14 : This system of nonlinear PDEs can be transformed, respectively, into a system of equation for stationary configurations Such a duality exists for Lorentz manifolds with a Killing symmetry when y 3 is a space like coordinate and y 4 = t is a time like coordinated. This duality simplifies various applications of the AFDM when we can redefine the procedure considered in the previous subsection for stationary nonholonomic configurations to certain time dependent ones. It allows use to prove decoupling properties and generate cosmological like solutions if there are known certain stationary configurations, or inversely to find certain stationary metrics as analogs of corresponding cosmological ones. 15 We can rewrite the nonlinear PDE (29)- (32) in an explicit decoupled form if we introduce the coefficients We can integrate such equations "step by step" for any generating function (x i , t) := e and sources h ϒ(x i ) and ϒ(x k , t), see next subsection. Nonlinear symmetries for generating functions and sources with effective cosmological constant: The system (34) with respective coefficients relates four functions (h 3 , h 4 , ϒ, ) when a very important nonlinear symmetry for locally anisotropic cosmological solutions and respective generating functions, ( (τ ), ϒ(τ )) ⇐⇒ ( (τ ), (τ )) can be found, This nonlinear symmetry allows us to introduce a new generating function (x i , t) and an (effective) cosmological constant (τ ) = 0, which can be applied both for generating exact off-diagonal solutions in explicit forms and elaborating on locally anisotropic cosmological scenarios with cosmological constants, for instance, considered in entropic gravity. 15 The LC-conditions (8) for stationary configurations transform into equations with coefficients depending on t,

Integrability of entropic quasiperiodic geometric flow equations
We generate and study geometric properties of two classes of generic off-diagonal solutions with elasticity and quasiperiodic structures of the system of nonlinear PDEs (23). The first one is for stationary configurations and the second one is considered for locally anisotropic cosmological models.

Off-diagonal cosmological solutions with elastic quasiperiodic structures
Applying the AFDM, we can construct cosmological solutions (which, in general, are locally anisotropic) of the entropic flow modified Einstein equations (23) with Nadapted sources h (τ ) = h (τ, x k ) and (τ ) = (τ, x k , t), see parameterizations for cosmological configurations in (26). Integrating "step by step" the system of the nonlinear PDEs (29)-(32) decoupled in the form (34), we obtain such d-metric coefficients for (2), The N-connection coefficients are computed, In these formulas, h [0] 3 (τ, x k ), 1 n k (τ, x i ), and 2 n k (τ, x i ) are integration functions encoding various possible sets of (non) commutative parameters and integration constants running on τ for geometric evolution flows. We can chose different generating data ( (τ, x i , t), ϒ(τ, x i , t)) or ( (τ, x i , t), (τ )) which are related by nonlinear differential/integral transforms (35), and respective integration functions. Such values should be chosen in explicit form following certain topology/symmetry/asymptotic conditions for some classes of exact/parametric cosmological solutions. The coefficients (36) and (37) define generic off-diagonal cosmological solutions if the corresponding anholonomy coefficients are not trivial. Such locally cosmological solutions are with nontrivial nonholonomically induced d-torsion and N-adapted coefficients which can be computed in explicit form. In order to generate as particular cases some well-known cosmological FLRW, or Bianchi, type metrics, we have to consider data of type ( (τ, t), ϒ(τ, t)), or ( (τ, t), (τ )), with integration functions which allow frame/ coordinate transforms to respective (off-) diagonal configurations g αβ (τ, t).
Let us analyze certain important nonholonomic evolution properties of above locally anisotropic cosmological solutions using the formulas for effective sources (22) with cosmological parameterizations (26). In N-adapted form, we obtain a system of equations with first order evolution derivatives ∂ τ when the v-part of vierbeinds depend on a time like coordinate y 4 = t, e  (x k , t), when the dependence on y 3 can be omitted because of Killing symmetry on ∂ 3 . We can consider frame transforms for generating effective sources, e f f In these formulas, we can prescribe any values for the matter sources tot ϒ μν (τ ) in a cosmological or spacetime QC model. Then, for simplicity, we can consider N-adapted diagonal configurations and integrate on τ and determine a cosmological evolution flow of g α (τ, x k , t) modelled as a nonholonomic and nonlinear geometric diffusion process. All geometric constructions are performed with respect to a new system of reference determined by e μ μ (τ, x k , t). We have to prescribe some locally anisotropic generating values [e μ μ (τ, x k , t), tot ϒ μν (τ, x k , t)] which are compatible with certain observational data, for instance, in modern cosmology and dark matter and dark energy physics.

Quadratic line elements for off-diagonal cosmological configurations with elastic flows
can be considered also as a generating function, for instance, for entropic quasiperiodic configurations. Using formulas (36), we find 2 = −4 (τ )h 3 (τ, r, θ, t) Introducing such values into the formulas for h a and in (36) and (37), we construct locally anisotropic cosmological solutions parameterized by d-metrics (2) with N-adapted coefficients (25), Such solutions posses a Killing symmetry on ∂ 3 and can be re-written in terms of η-polarization function functions for target locally anisotropic cosmological metrics g = [g α = η αgα , η a iN a i ] encoding primary cosmological data [g α ,N a i ].

Off-diagonal Levi-Civita entropic and quasiperiodic cosmological configurations
We can extract and model entropic flow evolution of cosmological spacetimes in GR. To satisfy the zero torsion conditions (8), see equations in footnote 15, let us consider a special class of generating functions and sources when, for  r, θ, t). Using these formulas, we conclude that the coefficient h 4 (τ ) =ȟ 4 (τ, x i , t) can be considered also as generating function for entropic cosmological solutions. For such LC-configurations, there are some parametric on τ functionsǍ(τ, x i , t) and n(τ, x i ) when the Nconnection coefficients are computed Summarizing above formulas, we construct new classes of locally anisotropic cosmological solutions as ub GR defined as subclasses of solutions (39) with zero torsion but with entropic quasiperiodic geometric flow evolution, or gener. funct.ˇ , sourceˇ ; Such cosmological metrics are generic off-diagonal and define new classes of solutions if the anholonomy coeffi-cients are not zero for N 3 k (τ = ∂ k n and N 4 i (τ = ∂ iǍ . They encode entropic quasiperiodic structures. We can analyze certain nonholonomic cosmological configurations determined, for instance, by data (ˇ ,ˇ , h [0] 3 ,ň k ), when ∂ k n → 0 and w i = ∂ iǍ → 0. Zero values can be fixed also by certain additional nonholonomic constraints. Choosing data (ˇ (τ, t),ˇ (τ, t), h [0] 3 = const,ň k = const), we can generate (off-) diagonal entropic metrics of Bianchi, or FLRW, types and generalizations to other type configurations g αβ (τ, t) in GR modified under geometric flow evolution.

Entropic quasiperiodic flows and cosmological solutions
The goal of this section is to consider physical implications of models with entropic and quasiperiodic flow evolution of locally anisotropic and inhomogeneous cosmological spacetimes.

The AFDM for entropic flow cosmological solutions
We outline the key steps on the AFDM for generating cosmological solutions with geometric flows and Killing symmetry on ∂ 3 . Considering a nonholonomic deformation procedure for a generating function g 3 (τ ) = h 3 (τ, x i , y 3 ) (36), cosmological constants (τ ) and sources h (τ ) = h (τ, x k ) and (τ ) = (τ, x k , t), see parameterizations (26) and nonlinear symmetries (35), we construct exact solutions of the system of nonlinear PDEs for emergent cosmology (34). Typical cosmological solutions of this class are parameterised Such quadratic line elements are time dual to the stationary ones which can be obtained via transforms (33).

Nonlinear PDEs for entropic quasiperiodic cosmology
We analyse two possibilities to transform the entropic flow modified Einstein equations (23) into systems of nonlinear PDEs (29)-(32) with generic off-diagonal or diagonal solutions depending in explicit form on a evolution parameter, a time like variable and two space like coordinates. In the first case, there are considered entropic quasiperiodic sources determined by some additive or general nonlinear functionals for effective matter fields. In the second case, respective nonlinear functionals determining quasiperiodic solutions for entropic configurations are prescribed for generating functions subjected to nonlinear symmetries (35). We also note that is possible to construct certain classes of locally anisotropic and inhomogeneous cosmological solutions using nonlinear/additive functionals both for generating functions and (effective) sources.

Cosmological solutions for entropic quasiperiodic sources
Cosmological configurations generated by additive entropic functionals for sources: For this class of cosmological solutions, we consider an additive functional for an entropic quasiperiodic source of type (τ, x i , t) (26), There is also an associated additive cosmological constant as (τ ) (40) (34) can be integrated following the AFDM explained in details in [55][56][57][58][59][60]63,64] (see sections on exact and parametric cosmological solutions in MGTs) and [33][34][35][36][37][38][39][40][41], for solutions with nonholonomic cosmological Ricci flows, and citations therein. Such generic off-diagonal cosmological solutions are parameterized in the form In such quadratic linear elements, we have to fix a sign of the coefficient h 3 (τ, x k , t) which describes relativistic flow evolution with a generating function with Killing symmetry on ∂ 3 determined by sources ( h (τ ), as (τ )). Such entropic and quasiperiodic flow solutions are of type (39) and can be re-written equivalently with coefficients stated as functionals of as (τ, x i , t) and as (τ, x i , t).
We can extract from off-diagonal d-metrics (42) certain cosmological LC-configurations determined by entropic quasiperiodic sources by imposing additional zero torsion constraints. Such anholonomy conditions restrict the respective classes of generating functions (ȟ 3 and/orˇ (τ, x i , t)) for n(τ, x i ),Ǎ(τ, x i , t)) and sources asˇ (τ ) (41) and as (τ ) (40), The d-metrics (42) and/or (43) define off-diagonal cosmological solutions generated by entropic quasiperiodic additive sources as (τ ) and/or asˇ (τ ). The terms (41) encode and model respectively contributions of standard and/or dark matter fields and effective entropic evolution sources. Such values can be can be prescribed in certain forms being compatible to observational data of cosmological (and geometric/entropic) evolution for dark matter distributions with possible quasiperiodic, aperiodic, pattern forming, solitonic nonlinear wave interactions. Cosmological solutions for nonlinear entropic quasiperiodic functionals for sources: Such classes of exact cosmological solutions can be generated by nonlinear quasiperiodic functionals for effective sources, , see functional dependencies in (41), subjected to nonlinear symmetries (35). Applying the AFDM, we construct entropic flow cosmological solutions of with nonlinear sources, For LC-configurations, we obtain For additive functionals for cosmological entropic and quasiperiodic sources, the formulas (44) and (45) transforms respectively into quadratic linear elements (42) and (43). Fixing a value τ 0 , we obtain cosmological solutions for Ricci solitons and MGTs.

Cosmological configurations with nonstationary entropic generating functions
In this section, the sources for (effective) matter fields and geometric flows are defined by arbitrary functions The quasiperiodic structure will be stated for some additive or general nonlinear functionals for the generating functions. Cosmological metrics with additive entropic generating functions: Such functionals are The values int are functionals on certain quasiperiodic (space time like QC or other type aperiodic, solitonic structures) given by functions ς and/or b subjected to conditions of type (12) and/or (13). The terms in the sum (46) correspond to the Lagrange densities (14) and energy-momentum tensors (15) nonholonomically parameterised by sources (26). Changing systems of references and coordinates, we can compute sums of functionals for sources of type (41) using nonlinear symmetries (35) considering an associated additive cosmological constant as (τ ) (40). For a (τ ) (46), the Eq. analogs of (31) and (32) in the decoupled system of nonlinear PDEs (34). Such equations and their solutions can be written equivalently in different forms with additive functionals of type a h 3,4 and/or a and respective nonlinear functionals for the coefficients in respective equations. Applying the AFDM as in the Sect. 4.3.1, we generate a class of parametric solutions of (34) parameterized in a form similarly to (39) for integration functions h [0] 3 (τ, x i ), 1 n k (τ, x i ) and 2 n k (τ, x i ). LC-configurations can be extracted from (47) imposing additional integrability conditions on coefficients resulting in zero nonholonomic torsion. Such solutions can be considered both for relativistic entropic Ricci solitons and in GR, being parameterized as in (40), We can consider small parametric decompositions and frame/coordinate transforms (in a more general context, we can elaborate a formalism of ηand ε-polarization functions as for BH solutions but with time like dependence for cosmological configurations) in order to relate new classes of solutions (47) and/or (48) to some well known (off-) diagonal cosmological metrics. Cosmological configurations with nonlinear entropic functionals for generating functions: Instead of additive generating functionals a (τ ) (46), we can work with nonlinear quasiperiodic generating functionals qp (τ ) = qp (τ, (τ )] (τ ), which can be solved together with other equations with decoupling (34). We obtain such solutions: where h [0] 3 (τ, x i ), 1 n k (τ, x i ) and 2 n k (τ, x i ) are integration functions.
For zero torsion constraints, we extract LC-configurations, The coefficients of d-metrics (49) and (50) can be chosen to be of necessary smooth class and involve certain entropic, quasiperiodic, stochastic sources and fractional derivative processes. Such nonholonomic deformation and generalized transforms can be constructing with changing the topological spacetime structure and modeling certain dark energy and dark matter effects as results of entropic quasiperiodic flows or nonholonomic deformations of certain prime cosmological solutions.

Emergent quasiperiodic cosmology from both generating functionals and sources
We can generate entropic flow cosmological solutions using generalized quasiperiodic nonlinear functionals both for generating functions, qp (τ ), and nonlinear functionals for (effective) sources, qp (τ ), see above formulas. Such data are connected via nonlinear symmetries of type (35), when Similar nonlinear symmetries exist for additive functionals both for the gravitational fields and (effective) sources, with (46) and (41) and can be written in a particular form a (τ) a 2 (τ ) = a 2 | a (τ )| − dt a 2 | a (τ )| * . Applying the AFDM (using data ( qp (τ ), qp (τ )), and/or, equivalently, ( qp (τ ), qp (τ ))), we construct multi-functional nonlinear entropic quasiperiodic cosmological configurations, where for integration functions For LC-configurations, we obtain multi-functional nonlinear generalizations of (49) and (50) modelling locally anisotropic and inhomogeneous solutions in entropic flow gravity and GR, The classes of cosmological solutions (51) and (52) describe off-diagonal entropic non-stationary configurations determined by multi-functional nonlinear quasiperiodic structures.

Cosmological metrics evolving in entropic quasiperiodic media
Generic off-diagonal entropic and quasiperiodic cosmological solutions are constructed in terms of η-polarization functions following the AFDM method for parametric τ and time like depending evolution. We consider a primary cosmological d-metricg (28) which can be diagonalized for a FLRW cosmological metric by frame/ coordinate transforms. 16 The cosmological entropic quasiperiodic y target metrics g(τ ) of type (27) are generated by nonholonomic deformations when overlined symbols are used for distinguishing cosmological d-metrics from stationary ones studied in previous sections. The quadratic line elements corresponding to target locally anisotropic and inhomogeneous cosmological metrics g(τ ) = g(τ, x k , t) can be written in terms of gravitational polarization functions, The η-coefficients will be constructed in some explicit forms determined by entropic quasiperiodics generation functions and/or effective sources as solutions of the system of nonlinear PDEs (34).

Cosmological evolutions generated by nonstationary entropic sources
We consider entropic sources of type qp (τ ) = qp (τ, (44) and (45) and compute the η-polarization functions following formulas η i (τ ) = e ψ(τ,x i ) /g i ; generating function where 1 n k (τ, x i ) and 2 n k (τ, x i ) are integration functions. Using η 3 (τ, x i , t) as a generating function, we can compute other types of generating functions of the same target cosmological d-metric subjected to nonlinear symmetries 16 In general, we can consider off-diagonal Bianchi anisotropic cosmological metrics or any cosmological solution in GR or MGTs. (35), For integrable generating functionals and sources, when the constructions (15) are subjected to target off-diagonal cosmological metrics (40) with zero torsion, we obtain In (53) and (54), the nonlinear functionals for the entropic quasiperiodic v-source and (effective) cosmological constant can be changed into additive functionals qp (τ ) → as (τ ) and qp (τ ) → as (τ ) which generates another classes of cosmological solutions.

Cosmology from nonstationary entropic generating functions
We can construct locally anisotropic and inhomogeneous cosmological solutions as nonholonomic deformations of some prime cosmological metrics when the coefficients of the d-metrics are determined nonlinear generating functionals (49). It is possible also to generate similar cosmological metrics by additive functionals a (τ ) (46) for prescribed families of effective sources (τ ) = (τ, x i , t) and cosmological constants (τ ). Using formulas for nonlinear symmetries (35), we compute (recurrently) corresponding nonlinear functionals, qp η 3 (τ, x i , t), or additive functionals, a η 3 (τ, x i , t), and related polarization functions, The coefficients of quadratic elements of type (39) are recurrently computed, In a series of our works [34][35][36][37][38][39][40][41]78] , we considered various modifications of the F-and W-entropy functionals for constructing (non) commutative and/or supersymmetric geometric evolution models; investigating nonlinear (fractional and/or locally anisotropic) stochastic and kinetic processes in curved spaces; and elaborating on cosmological scenarios encoding memory and quantum interactions of gravitational and (effective) matter fields. Here we note that the concepts of complex manifolds, supermanifolds and noncommutative geometry, are defined by geometric constructions with nonholonomic distributions on curved spaces. In such a general context, for various entropic spacetime and emergent gravity models, and in geometric flow evolution theories with nonholonomic constraints, we have to develop an unified geometric formalism for metric-affine spaces, generalized Finsler-Lagrange-Hamilton geometry, almost Kähler and noncommutative geometries etc. Such constructions were performed for Hořava-Lifshits, f (R), R 2 , and other types MGTs (see reviews [68][69][70]79,80]) with developments for models of thermodynamic/entropic/entanglement and emergent gravity [4][5][6][7]16,17,[19][20][21][22][23][24][25]31,[65][66][67]; and for locally anisotropic kinetic and thermodynamic theories on curved spaces (see [38,71,78,[81][82][83] and references therein) etc. We have found that using certain classes of nonholonomic variables the (relativistic) geometric flow evolution and Ricci soliton equations, and related motion equations in entropic and other type MGTs, can be decoupled and integrated in some very general forms. This allows us to construct various classes of exact and parametric solutions with generic off-diagonal metrics and generalized connections. The coefficients of new classes of stationary and (in general, locally anisotropic and inhomogeneous) cosmological solutions depend on all spacetime and associated phase space (kinetic and/or thermodynamic) coordinates via generating and integration functions and various types of commutative and noncommutative parameters and integration and physical constants. Such geometric and analytic techniques of constructing exact solutions in geometric flow evolution and MGTs have been developed in the framework of the so-called anholonomic frame deformation method, AFDM [55][56][57]. For details, examples of exact solutions, and various applications, we cite [33,[38][39][40][41]41,[58][59][60]63] and references therein.
It should be noted that there were not elaborated corresponding topological methods and a well-defined analytic formalism for investigating geometric evolution equations of metrics of pseudo-Euclidean signature. A number of such conceptual and fundamental issues in nonlinear functional analysis and the geometry of Lorentz manifolds have not addressed or solved by mathematicians. The standard geometric flow paradigm was proposed as the Hamilton-Perelman program for Riemannian metrics defining entropy type functionals and deriving nonlinear evolution equations.
To elaborate realistic relativistic physical models we deal with Ricci tensors which in the limit of weak gravitational/elastic flows approximate to the d'Alembert (wave) operator and not to the Laplace (diffusion) one used for Euclidean signatures. So, the original approach to the topology and geometric flows of Riemannian metrics has to be generalized in certain relativistic and nonholonomic forms which are compatible with modern experimental particle physics data and observational cosmology. The Poincaré-Thurston conjecture can be formulated and proven for nonrelativistic evolution of any 3-d space like hypersurface. In this and partner [32] works, we advocate that using E. Verlinde conjecture on elastic emergent gravity, we can elaborate on generalizations of nonholonomic Ricci flow theories as certain models of relativistic flow evolution. Such models are determined by extensions of Perelman's functionals for 3-d Riemannian metrics to certain modified 4-d Fand W-entropy nonholonomic analogs which are extended on a time like coordinate and/or a temperature like evolution parameter. For certain approaches with rich gravitational vacuum structure, we work with geometric relativistic kinetic/hydrodynamic/thermodynamic models (see details in [38]) or (for instance, in this work) with a τ -parametric theory describing geometric entropic flows determined by relativistic and elastically spacetime modified nonholonomic F-and W-functionals.
As it was mentioned above, there is not yet formulated a rigorous mathematical approach to the theory of relativistic geometric/entropic flows of metrics with Lorentz signature and generalized connections. Nevertheless, we shown that such theories are characterized by certain classes of generalized R. Hamilton equations with effective parametric sources which may encode entropic and quasiperiodic structures (these are necessary for explaining, for instance, the complex structure of dark matter and energy in modern cosmology). Using the AFDM, we proved that such geometric flow evolution equations and their Ricci soliton variants can be decoupled and integrated in very general forms. In Sect. 5, we constructed and analyzed possible physical implications of respective entropic locally anisotropic cosmological solutions. A series of our former results on astrophysical and cosmological models with quasiperiodic, pattern forming, quasicrystal time like structures, see [59][60][61][62][63][64] and references therein, were used for developing in this work similar models for the entropic geometric flow and gravity theories. Such exact and parametric solutions provide also explicit examples of entropic gravity models developed in a phenomenological manner in [6,7,65]. So, we conclude that the E. Verlinde conjecture on entropic character of gravity can be related to a relativistic extension of the Poincaré-Thurston conjecture. Even such geometric ideas have not been proven as explicit theorems in modern geometric analysis, there are rigorous exact solutions of respective systems of nonlinear PDES which support such entropic gravity and flow evolution ideas.
In Sects. 4 and 5, we shown how to construct in explicit form entropic and quasiperiodic solutions for relativistic geometric flows, nonholonomic Ricci solitons and generalized gravitational field equations. Such a techniques was elaborated similarly to MGTs with quasiperiodic structure (constructed and studied in this and our partner works [32,63,64]) and involves generic off-diagonal metric and nonholonomically deformed non-Riemannian linear and nonlinear connections. We emphasized, see also [41], that such configurations are not characterized, in general, by certain entropyarea, holographic or duality conditions. As a consequence, it is not possible to elaborate on thermodynamic models of entropic MGTs and physical properties of their exact or parametric solutions using only the concepts related to the Bekenstein-Hawking entropy. We consider that there is an alternative and more general way when stationary and cosmological solutions in geometric/entropic flow evolution theories, MGTs and GR, can be defined and characterized by nonholonomic deformations of Perelman's W-entropy. Such constructions are similar to the well-known results on relativistic locally anisotropic thermodynamics and kinetics [38,71] and can be generalized for emergent classical and quantum gravity theories.
Finally, we note that our geometric/entropic flow approach to MGTs provides new mathematical methods and applications in the theory of classical and quantum informatics, for research of quantum systems with entanglement, models of quantum and emergent gravity, and accelerating cosmology and dark energy/matter interactions etc. Our further research programs are related to developments in such directions.