Deforming Black Hole and Cosmological Solutions by Quasiperiodic and/or Pattern Forming Structures in Modified and Einstein Gravity

We elaborate on the anholonomic frame deformation method, AFDM, for constructing exact solutions with quasiperiodic structure in modified gravity theories, MGTs, and general relativity, GR. Such solutions are described by generic off-diagonal metrics, nonlinear and linear connections and (effective) matter sources with coefficients depending on all spacetime coordinates via corresponding classes of generation and integration functions and (effective) matter sources. There are studied effective free energy functionals and nonlinear evolution equations for generating off-diagonal quasiperiodic deformations of black hole and/or homogeneous cosmological metrics. The physical data for such functionals are stated by different values of constants and prescribed symmetries for defining quasiperiodic structures at cosmological scales, or astrophysical objects in nontrivial gravitational backgrounds some similar forms as in condensed matter physics. It is shown how quasiperiodic structures determined by general nonlinear, or additive, functionals for generating functions and (effective) sources may transform black hole like configurations into cosmological metrics and inversely. We speculate on possible implications of quasiperiodic solutions in dark energy and dark matter physics. Finally, it is concluded that geometric methods for constructing exact solutions consist an important alternative tool to numerical relativity for investigating nonlinear effects in astrophysics and cosmology.


Introduction
A series of recent works has been devoted to constructing black hole, BH, and cosmological solutions encoding quasiperiodic gravitational and matter fields interactions and geometric flow structures [1,2,3,4,5,6,7]. We elaborated a new geometric and analytic method for generating exact and parametric solutions of (modified) Einstein equations characterized by generic off-diagonal metrics and (generalized) connections with coefficients depending on all spacetime coordinates. It is called the anholonomic frame deformation method, AFDM; for reviews of results, see [8,9,10,11,12,13,14]. This geometric formalism provides an alternative to numerical methods and simulations in modern gravity [15,16,17,18], see a review of such results in [19]. In our works, there were considered various geometric applications in particle physics and cosmology and heterotic superstring theory, and studied various models of commutative and noncommutative and/or supersymmetric modified gravity theories, MGT, and general relativity, GR, see [20,21,22,23,24,25,26] and references therein.
The AFDM was formulated by developing a geometric techniques of decoupling and integrating systems of nonlinear partial differential equations, PDEs. The general goal was to study nonlinear dynamical and/or evolution properties of (modified) Einstein equations, when various classes of solutions are generated by nonholonomic (i.e. non-integrable, equivalently, anholonomic) frame transforms with deformed symmetries and distortions of metric and connection structures. In this paper, we consider that readers are familiar with main concepts and methods of differential geometry and topology, functional analysis and theory of partial differential equations, PDEs, and ordinary differential equations, ODEs, with applications in mathematical relativity, astrophysics and cosmology (summarized in [27,28,29,30,31]). The most important feature of the AFDM is that it allows us to construct in direct form exact solutions of nonlinear systems when there are not involving simplified diagonal metric ansatz resulting in ordinary differential equations, ODEs. We can study singular and nonsingular metrics and generalized connections, quasiperiodic (non) stationary gravitational and matter field structures, nonlinear waves and solitons, deformed horizons, nontrivial topological structures etc. Additionally, one can be analysed the conditions when and if some classes of "general solutions" (with less clear physical properties/ implications) may have certain (non) smooth limits and/or nonholonomic constraints to some well-known and more "simple" physically important solutions.
We emphasize that new classes of exact and parametric quasiperiodic and pattern forming solutions in MGTs and GR generated following the AFDM (in this and partner works [4,6,7]) can not be constructed by applying former methods with diagonal metric ansatz, for instance, depending only on a "radial" like coordinate, for BHs, or on a time like coordinate, for homogeneous cosmological solutions. Prescribing from the very beginning some "simplified" diagonal ansatz, we can transform the gravitational and matter field equations into certain systems of nonlinear ODEs which can be solved in general forms with metrics parameterized by integration constants. This way, we lose various possibilities to find more general classes of nonlinear solutions, for instance, with solitonic hierarchies, pattern forming structure etc. which are determined by generating and integration functions. More realistic descriptions and understanding of physical properties of realistic nonlinear gravitational and matter field systems are possible in terms of generalized ansatz for coefficients of metrics, frames, and (non) linear connections, and off-diagonal solutions of PDEs depending on a maximal possible number of spacetime coordinates and various (non) commutative continuous/ discrete parameters.
The solutions for gravitational and matter fields with quasiperiodic structures are defined by generic offdiagonals metrics, nonlinear and (generalized) connections, and (effective) matter sources with coefficients determined by respective classes of generating and integration functions depending on all spacetime coordinates. 1 Prescribing for physical and geometric objects corresponding smooth / singularity/ symmetry conditions, solving corresponding Cauchy problems, and/ or satisfying necessary boundary conditions, we can model various (non) commutative and/or semi-classical and quantum stationary configurations and/or nonlinear gravitational and matter field interactions. Such models of generalized off-diagonal spacetimes are characterized by nonlinear / discrete / continuous (semi) classical symmetries; there are computed locally anisotropically polarized (physical) constants and stated values of integration parameters; in general, such constructions can be with (non) singular and/or topologically nontrivial structures. In mentioned partner works [4,6,7] (see also citations therein) and related papers on applications of the AFDM, various examples of exact and parametric solutions, for instance, with soliton distributions and nonlinear waves, partial derivative and gravitational diffusions processes, geometric flows, and quasicrystal like structures, locally anisotropic cosmological configurations etc., have been considered.
The goal of this paper is to analyze the conditions when certain data for generating and integration functions, (effective) sources, prescribed symmetries for integration functions and various type parameters, define offdiagonal deformations of physically important solutions. We study nonholonomic deformations of BH and cosmological metrics into new classes of exact solutions encoding quasiperiodic and pattern forming gravitational and effective matter fields. Our work reconsiders a series of ideas, methods, and results elaborated in condensed matter physics and develop them with applications in modern gravity, cosmology, and astrophysics. On relevant former results, we cite a series of papers for the physics of quasicrystals and Penrose-like telings, see [32,33,34,35,36,37,38,39,40,41,42,43,44,45]). It is considered that various web like, filament, quasiperiodic, aperiodic, singular and nonsigular cosmological evolution, solitonic like distributions and nonlinear wave structures are distinguished /found in modern cosmological observational data [46,47,48,49,50,51,52,53,54,55]. Such generic locally anisotropic and inhomogeneous cosmological configurations can not be explained exhaustively using only higher symmetry solutions (for instance, Friedmann-Lemaître-Roberstson-Walker, FLRW, and/or Bianchi type metrics) depending only on time like coordinates and constructed following methods of the theory of ODEs. Applying the AFDM, we are able to construct more general classes of exact solutions in MGTs and GR which allows us to model more rich spacetime structures and nonlinear field interactions as solutions of systems of nonlinear PDEs. In this work, it is shown how such nonlinear configurations are characterized by effective free energy functionals and nonlinear evolutions equations, which can be applied to mimic dark energy and dark matter effects, see also [44,45,2,3,4,5,6,7]. We study quasiperiodic solutions not involving small deformation parameters.
The content of the paper is organized as follows: In section 2, we provide a brief review of the AFDM (relevant details with necessary N-adapted coefficient formulas and proofs are presented in Appendix A). Free energy functionals for generating functions resulting in quasiperiodic and pattern forming structures for smooth, singular and BH like structures, and cosmological solutions etc. are defined and analysed in section 3. We construct a series of examples of stationary and nonstationary exact solutions describing quasiperiodic deformations of BH solutions in MGTs and GR in section 4. There are analyzed different classes of off-diagonal solutions determined by nonlinear or additive functionals for generating functions and effective sources. Various classes of locally anisotropic and/or inhomogeneous cosmological quasiperiodical solutions are constructed in general form and studied in section 5. Such cosmological evolution models are determined by nonstationary quasiperiodic generating functionals and effective sources. Finally, we discuss and provide concluding remarks in section 6.

A Brief Review of the AFDM
We outline and reformulate the AFDM as a geometric formalism for constructing exact generic off-diagonal stationary and cosmological solutions in MGTs and GR determined by generation and integration functions and effective matter sources, and related polarization functions. The method and most important formulas are summarized in Tables 1-3. Such constructions will be applied in order to extend the approach for solutions encoding quasiperiodic and/or pattern forming, solitonic distributions and nonlinear wave gravitational and matter field structures. Details, examples with small parametric decompositions and rigorous mathematical proofs are contained in [8,9,10,11,12,13,14,20,21,22,23,24,25,26], see also a summary of necessary N-adapted formulas in Appendix A.
The 2+2 splitting (2) is nonholonomic (equivalently, non-integrable, or anholonomic). This follows from the conditions that the co-basis e α = (dx i , e a ) (3) is dual to e α = (e i , e a ) = (e i = ∂/∂x i − N a i (u)∂/∂y a , e a = ∂ a = ∂/∂y a ), which satisfies nonholonomy conditions e [α e β] := e α e β − e β e α = C γ αβ (u)e γ , with anholonomy coefficients C γ In above formulas, we can consider necessary type frame transforms e α = e α α ′ (u)du α ′ , when g α ′ β ′ (u) = g αβ e α α ′ e β β ′ . A set of coefficients N = {N a i } from (2) and/or (4) states a N-adapted decomposition of the tangent Lorentz bundle T V into conventional horizontal, h, and vertical, v, subspaces, when the Whitney sum define a nonlinear connection (N-connection) structure. A Lorentz manifold V is called nonholonomic if it is enabled with a nonholonomic distribution. In a particular case, we shall write V = (V, g, N) for a nonholonomic spacetime manifold enabled with N-connection structure N (6). This concept was introduced and studied in details in generalized Finsler geometry and geometric and physical models on (co) vector/ tangent bundles and their supersymmetric/noncommutative generalizations, see [12] and references therein. For (pseudo) Riemannian manifolds, a conventional h-/ v-splitting can be considered, for instance, in order to state a fibered structure with y a coordinates and N-adapted frames (4). Such nonholonomic structures and deformations of geometric/ physical objects (like metrics, linear connections, (effective) sources etc.) can be introduced in such forms which allows decoupling of gravitational and matter field equations in MGTs and GR, [1,2,3,4,5,6,7]. 3 It is possible to model a nonholonomic deformation with η-polarization functions,g → g, of a 'prime' metric, g, into a 'target' metric g = g (2), if where the N-elongated basis (4) We shall subject a g to the condition that it defines a solution of (modified) Einstein equations. A general prime metric in coordinate parametrization of type (1), can be also represented in N-adapted form It can be, or not, a solution of some gravitational field equations in a MGT or GR. In this work, we shall be interested in two physically important cases wheng (8) defines a BH solution (for instance, a vacuum Kerr, or Schwarzschild, metric), or a Friedman-Lemaître-Robertson-Walker (FLRW) type metric. For such diagonalizable metrics (the off-diagonal structure of the Kerr metric is determined by rotation frames and coordinates), we can always find a coordinate system whenN b i = 0. In order to avoid singular noholonomic deformations, it is convenient to construct exact solutions of necessary type gravitational equations with nontrivial conventional "polarization" functions η α = (η i , η a ), η a i , and nonzero coefficientsN b i (u). This can be achieved by considering necessary type frame/ coordinate transforms. Having constructed an explicit form 2 If the anholonomy coefficients are nontrivial (i.e. the N-adapted frame is not integrable and can not be transformed into a coordinate one), the respective equivalent metric parameterizations (1) and (2) are generic off-diagonal. 3 We shall follow the Einstein rule on summation on "up-low" cross indices if there is not a contrary statement for an explicit formula. One shall not be considered summation on repeating "low-low", or "up-up" indices. In our works , we elaborated a system of "N-adapted notations" with boldface symbols for manifolds and fiber bundles enables with a nontrivial N-connection structure N. Details on nonholonomic differential geometry and N-connections are explained in [11,12,13,14,20,21,22,23,24,25,26] and references therein. If the anholonomy coefficients C γ αβ in (5) are nontrivial, a metric g αβ (1) can not be diagonalized in a local finite, or infinite, spacetime region with respect to coordinate frames. Such metrics are called generally off-diagonal and characterized by six independent nontrivial coefficients from a set g = {g αβ (u)}. A nonholonomic frame is holonomic if all corresponding anholonomy coefficients are zero (for instance, the coordinate frames). 4 in this paragraph, we do not consider summation on repeating indices if they are not written as contraction of "up-low" ones a d-metric (7), we can study the existence and geometric/ physical properties of solutions, for instance, when η α → 1 and N a i →N a i , or if η α = 1 and/orN a i = 0 are imposed as some noholonomic constraints. 5 We denote certain nonholonomic deformations of a prime d-metrics into a target one asg → g = [g α = η αgα , η a iN a i ]. Here we emphasize that one constructs, in general, different classes of solutions for nonlinear systems of PDEs, if such approximations are considered in (modified) Einstein equations before finding solutions, or at the end (after a class of solutions has been constructed in explicit form). This is an important property of nonlinear dynamical / evolution systems which can be subjected to additional nonholonomic constraints.
The standard formulation of (pseudo) Riemannian geometry is in terms of the Levi-Civita, LC, connection ∇, which (by definition) is metric compatible and with zero torsion. Nevertheless, any Lorentzian manifold V can be enabled additionally (following certain geometric/ physical principles) with other types of linear connection structures considering any D = ∇ + Z, characterized by a respective distortion tensor Z. A general linear connection D can be metric noncompatible and/or with nontrivial torsion. On a nonholonomic V, we can consider a N-adapted variant of linear connection structure, called a distinguished connection, d-connection, D = (hD,vD), which is defined as a metric-affine (linear) connection preserving under parallel transports the N -connection splitting into h-and v -subspaces. We can define and compute for any d-connection D (in standard form) the torsion tensor, T = {T α βγ }, and the curvature tensor, R = {R α βγδ }, where the coefficients can be written in N-adapted form with respect to necessary tensor products of bases (4) and their duals.
On a nonholonomic V, we can work in equivalent form with two different linear connections: where D = (h D, v D) is completely defined by g for any prescribed N-connection structure N. In these formulas, we denote (respectively, for D and ∇) the torsions, T and ∇ T = 0, and curvatures, R = { R α βγδ } and ∇ R = {R α βγδ }, which can be defined and computed in coordinate free and/or coefficient forms. 6 There is a canonical distortion relation The distortion distinguished tensor, d-tensor, is an algebraic combination of the coefficients of the corresponding torsion d-tensor T = { T α βγ } of D. Readers may consult Appendix A for a summary of relevant formulas and some details on the geometry Lorentz manifolds with N-adapted 2+2 variables.
We can define two different Ricci tensors, Ric = { R βγ := R γ αβγ } and Ric = {R βγ := R γ αβγ }, when Ric is characterized by h -v N-adapted coefficients, Respectively, there are two different scalar curvatures, R := g αβ R αβ and R := g αβ R αβ = g ij R ij + g ab R ab . Following the two connection approach (9), we conclude that the (pseudo) Riemannian geometry can be equivalently described by two different geometric data (g, ∇) and (g, N, D). Using the canonical distortion relation (10), we can compute respective distortions of curvature and Ricci tensors, 5 It is possible to keep a physical interpretation of a target metric g (7) with generic off-diagonal terms, as an "almost" BH, or FLRW cosmological, like metric by constructing parametric solutions with small nonholonomic deformations on some constant . Such 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 general data (ηα, η a i ). For certain subclasses of solutions, we can consider that εi, εa, ε a i ∼ ε, when one small parameter is considered for all coefficients of nonholonomic deformations. It is possible to works with mixed types of solutions and model only small diagonal deformations εi, εa, ∼ ε of metrics, for some general η a i , or to work with some nontrivial ηα but ε a i ∼ ε. In particular, one can be generated various classes of nonholonomic small deformations of solutions like in Refs. [1,2,3,4,5,6,7]. In this work, our goal is to construct and study quasiperiodic solutions with generation functions and sources without small parameters. 6 It should be noted that the well known LC-connection ∇ is a linear one but not a d-connection because it does not preserve, under general frame/coordinate transforms, a h-v-splitting. T is a nonholonomically induced torsion determined by nontrivial values (C γ αβ , ∂ β N a i , g αβ ) which is different from the Einstein-Cartan, or string theory, when additional field equations and sources for torsion fields are considered.
with corresponding distortion tensors ∇ Z and Zic. Such formulas motivate application in GR and MGTs of nonholonomic geometric methods with multiple metric and connection structures, and adapted frames. For certain well defined nonholonomic configurations, various types of gravitational and matter field equations rewritten in nonholonomic variables (g, N, D) can be decoupled and integrated in some general forms. 7

MGTs in N-adapted variables and decoupling of (modified) Einstein equations
We show how the gravitational field equations can be written in nonholonomic variables and consider the key steps for generating exact solutions in Tables 1, 2, and 3.

Gravitational field equations for the canonical d-connection
Various models of modified gravity and acceleration cosmology (see [1,2,3,4,5,6,7] and references therein) have bee elaborated for the so-called Starobinsky type R 2 gravity [47]. For well defined conditions on a class of conformal transforms, the gravitational field equations for the so-called quadratic gravity are is equivalent to the Einstein gravity with scalar field sources. In nonholonomic variables (g, N, D) and for interactions with scalar fields defined by Lagrange density m L(g, N, φ), the action for such theories can be written in the form In this formula the Planck mass M P is determined by the gravitational constant; the Lagrange density m L[φ] and the action m S = d 4 u |g| m L are postulated in such forms which, for simplicity, allow to find explicit solutions of (modified) Einstein equations. This is possible if there are considered m L[φ] depending only on the coefficients of a metric field and not on their derivatives. Applying a N-adapted variational calculus, the energy- Following N-adapted variations of S (13), one derives such gravitational field equations where := D 2 = g µν D µ D ν and certain constraints/ conditions with D | T →0 = ∇ can be imposed as in GR.
The equations (14) for MGTs can be considered as modified Einstein equations when ∇ is changed into D and the standard energy momentum tensors are nonholonomically deformed into N-adapted sources (15).
The key steps for constructing exact solutions using two different methods [1) with reduction to ODEs or 2) with integration of PDEs] of constructing exact solutions of gravitational and matter field equations in MGTs and GR are briefly reviewed in Appendix and summarized in the Tables 1, 2, and 3. In references [27,28,29,30,31], there are provided details and main examples for constructing physically important solutions by using diagonal ansatz, for instance, with spherical/cylinder symmetries reducing the (generalized) Einstein equations to certain systems of nonlinear ODEs. Proofs of results and a number of examples how the AFDM should be applied in order to generate exact solutions of gravitational, matter field and evolution nonlinear systems of PDEs are considered in [8,9,10,11,12,13,14,20,21,22,23,24,25,26].

(Off-) diagonal metric ansatz and ODEs and PDEs in (non) holonomic variables
We explain how using holonomic 3+1 and 2+2 nonholonomic variables and corresponding ansatz it is possible to transform gravitational field equations in MGTs and GR into respective systems of nonlinear ODEs and PDEs is summarized in Table 1. N-adapted formulas are provided in Appendix A.

Decoupling and integration of gravitational filed equations and stationary solutions
The key steps of AFDM for generating stationary off-diagonal exact solutions of (modified) Einstein equations described in appendices A.2.1, A.3 are summarized in Table 2. Following a nonholonomic deformation procedure, for instance, for a generating function h 4 (r, θ, ϕ), we construct a class of off-diagonal stationary solutions with Killing symmetry on ∂ t determined by sources ( h Υ, Υ) and effective cosmological constant Λ, Such solutions are, in general, with nontrivial nonholonomically induced torsion (A.2). They can be re-defined equivalently in terms of generating functions Ψ(r, θ, ϕ) or Φ(r, θ, ϕ), see (A. 19). LC-configurations in GR can be extracted for additional zero torsion constraints with a more special class of "integrable" generating functions (ȟ 4 , andΨ(r, θ, ϕ) and/orΦ(r, θ, ϕ)) for respective sourcesΥ and Λ (A.20), In terms of η-polarization functions, such d-metrics and N-connections, can be parameterized to describe nonholonomic deformations of a primary (for instance, BH) d-metricg into target generic off diagonal stationary solutions g, see (7), asg → g = [g α = η αgα , η a iN a i ].
Solutions for polarization funct.

Decoupling and integration of gravitational PDEs generating cosmological solutions
In Table 3, we state the key steps of the AFDM for generating off-diagonal locally anisotropic solutions of (modified) Einstein equations described in appendices A.2.2 and A.3.

N-connections, zero torsion
Example of a prime cosmological metricg any frame transform of a FLRW or a Bianchi metrics Solutions for polarization funct.
Polariz. funct. with zero torsion Applying the nonholonomic deformation procedure (for simplicity, we consider metrics determined by a generating function h 4 (x k , t), we construct a class of generic off-diagonal cosmological solutions with Killing symmetry on ∂ ϕ determined by sources, h Υ and Υ, and an effective nontrivial cosmological constant Λ, Such locally anisotropic and inhomogeneous cosmological solutions are, in general, with nontrivial nonholonomically induced torsion (A.2). This class of solutions can be re-defined equivalently in terms of generating functions Ψ(x k , t) and/or Φ(x k , t), see (A.34).
Cosmological configurations in GR can be extracted by imposing additional constraints for zero torsion by restricting the class of "integrable" generating functions (h 4 =ȟ 4 , andΨ(x k , t) and/orΦ(x k , t)), for respective types of sourcesΥ and Λ, as in (A.35), In terms of η-polarization functions, the coefficients of cosmological d-metrics and N-connections can be parameterized to describe nonholonomic deformations of a primary (for instance, a FLRW) d-metricg into target generic off diagonal cosmological solutions g(x i , t) → g(t), see (7), asg → g=[g α = η αgα , η a iN a i ].

Quasiperiodic/ Pattern Forming Generating Functions and Sources
The goal of this section is to analyse a series of examples of quasiperiodic and/or pattern forming distributions and related free energy functionals and systems of PDEs evolution equations on a nonholonomic Lorentzian manifold V, which have certain analogy in condensed matter physics (see, for instance, [44,45,48,49]). For small parametric dependencies, such gravitational and matter field (non) stationary distributions and nonlinear field interactions were studied in our partner works [4,6,7]. In this paper, the approach is elaborated for generating and integration polarization functions and (effective) sources stated in certain general forms without additional assumptions on decompositions on small parameters.

Nonholonomic 3+1 distributions with quasiperiodic/ pattern forming structures
Let us consider necessary smooth classes of functions q = q(x i , y 3 ), for space like distributions, and q = q(x i , y 4 = t), for locally anisotropic cosmological configurations, defined respectively in N-adapted coordinates on open regions of U ⊂ V and U ⊂ V. Such values will be used as generating functions and/or (effective) sources for different models of quasiperiodic and/or pattern forming spacetime structures. Additionally to a nonholonomic 2+2 splitting which allows us to decouple systems of nonlinear PDEs and construct exact solutions of (modified) gravitational equations, we shall consider a 3+1 splitting with local coordinates parameterized in the form u α = (uì, t), where the space like coordinates are uì = (x i , y 3 ), withì,j,k, ... = 1, 2, 3. In GR, models with 3+1 spacetime decomposition were elaborated with the aim to introduce values similar to the energy and momentum and thermodynamical like characteristics for gravitational and scalar fields, see details in Refs. [27,29]). In our recent works [4,6,7,50,51], double 2+2 and 3+1 spacetime decompositions (and various respective extra dimension generalizations) were considered for constructing generic off-diagonal solutions in MGTs and GR encoding quasi-periodic and aperiodic structures.
Any metric and/or equivalent d-metric structures parameterized in a form (2) or (7) can be re-written in certain forms with nonholonomic 3+1 splitting, For such configurations, a 4-d metric g can be considered as an extension of a 3-d metric bìj = diag(bì) on a family of 3-d hypersurfaces Ξ t parameterized by coordinate t considered as a parameter, and when b 3 = h 3 andN 2 (u) = −h 4 is defined by a lapse functionN (u). We can impose additional conditions in order to transform stationary d-metrics (A.19), or locally anisotropic cosmological d-metrics (A.34), into respective 3+1 versions (16 ). Correspondingly, we shall work with a Killing symmetry on ∂ 4 and lapse functions of typeN (x k , y 3 ), or with a Killing symmetry on ∂ 3 and lapse functions of typeN (x k , t). Such a decomposition results in a representation D = ( 3 D, t D), where 3 D defines the action of the canonical d-connection covariant derivative on space like coefficients and t D of time like coefficients. For LC-configurations, the covariant derivative operator splits as ∇ = ( 3 ∇, t ∇), when the action on a scalar field q(u) can be parameterized via frame/ coordinate transforms as ∇q = ( 3 ∇q, t ∇q = ∂ t q = q * ).

Many pattern-forming nonlinear gravitational and matter fields systems
In condensed matter physics, models with tree-waves interactions, 3WIs, for many pattern-forming systems were elaborated for explaining experimental observations of certain microscopic and quasi-classical quantum systems, see [44] and references therein. Modern cosmological data show a very complex web like quasiperiodic and/or aperiodic like structure formation, geometric anisotropic evolution and nonlinear gravitational and matter field interactions, including dark energy and and dark matter configurations [46,52]. We can apply very similar mathematical methods for geometric modeling of quasi-crystal matter or (super) galactic clusters and 3-d distributions of dark energy and dark matter and generating solutions of systems on nonliner PDEs describing such physical systems.
A prime pattern-forming field can be taken in the form v l (t)e ic l ·u + higher order terms (17) defining in flat spaces 3WIs involving two comparable wavelengths. We consider systems with two wave numbers k = 1 and k = c, 8 for 0 < c < 1, when 3WIs are modelled in two forms: 1. two waves (with wave number 1, on the outer circle) interact nonlinearly with a wave on the inner circle, with wave number c (for instance, wave vectors configurations k 1 and k 2 interact with c 1 = k 1 + k 2 ); 2. for 1/2 < c < 1, two waves on the inner circle interact with another wave on the outer circle (for instance, wave vectors configurations c 2 and c 3 interact with k 1 = c 2 + c 3 ).
For (17), we can parameterize the coefficients and respective values in such forms that there are modelled two types of 3WIs. There are involved (in the case 1 above) a triad of wave vectors k 1 , k 2 and c 1 = k 1 + k 2 , when amplitudes are subjected to equations z * 1 = µz 1 +Q zv z 2 v 1 +cubic terms, z * 2 = µz 2 +Q zv z 1 v 1 +cubic terms, and v * 1 = ξz 1 +Q zz z 1 z 2 +cubic terms. (18) Other configuration can be defined (in the case 2) for wave vectors c 2 , c 3 and In these formulas, the coefficients µ and ξ determine, respectively, the growth rates of amplitudes corresponding to wave numbers 1 and c; and, for instance, Q zv and Q zz are quadratic elements. For simplicity, we omit the cubic terms even they play also an important role in the nonlinear dynamics of such waves and result in observable effects both in condensed matter physics and at cosmological distances. Similar equations can be replicated for all possible combinations of modes describing 3WIs. Different stationary and dynamical models with nonlinear wave interactions are generated for certain ranges of values of coefficients µ, ξ, Q zv , Q zz , ... and the signs of quadratic coefficients for products of type Q zv Q zz etc. were analyzed in [44] (on experimental data in condensed matter physics, see references therein).
In a more general context, we can consider nonlinear deformations of prime waves (17),q(x i , t) → P q(x i , t), when the target field P q = q describes pattern forming configurations as solutions of nonlinear PDE, with summation on up-low indexì. This equation can be parameterized with a linear part L acting on a mode e ikx with an eigenvalue σ(k) specified by σ(1) = µ and σ(c) = ξ. Such a specification controls growth rates of the modes of interest; with dσ/dx = 0, for k = 1 and k = c, and σ(0) = σ 0 < 0, one controls the depths of minimum between k = 1 and k = c. For even functions σ of k, we can consider a 4th order polynomial on k 2 , The linear operator L in the linear part of (20) is defined by replacing k 2 by − 3 Dì 3 Dì. The nonlinear terms in that PDE can be boosted as quadratic and cubic combinations of q and its derivatives, which can re-parameterized as certain nonlinear deformations ofq (17). A standard weakly nonlinear theory expresses the values 1 Q, 2 Q, and 3 Q as certain products Q zv Q zz etc. when the solutions of (18) and (19) are nonlinearly mixed which results in different signs and constant coefficients.
Using P q = q as a solution of (20), we can generate various type pattern forming configurations in a nonholonomic spacetime. For instance, there are possible 3 important results (for simplicity, we take values of constants reproducing the results explained in details in [44] and references therein; there are cited also some papers analyzing similar structures in modern gravity, cosmology and astrophysics):

A bifurcation pattern can be seen for
There are: z hexagons, for k = 1; v hexagons for k = c; certain spatiotemporal chaos, STC, with mixed patterns defined by v stripes with patches of z rectangles; two super lattice patterns (the first one is with 6 modes on the outer circle and with 12 modes on the inner; the second one is with 6 modes on the inner circle and with 12 ones on the outer). Here it should be noted that the scales stated by µ and ξ are not uniform and that we can model additionally various diffusion and additional nonlinear wave interactions [9, 10].
2. One can be found a STC pattern with µ = ξ = 0.000707, when the correlation length is about 1-2 wave-lengths. We can ad new terms with fractional chaos and diffusion [21,22].
3. Another interesting configuration with two critical circles can be reproduced in the power section [53].
Solutions of (20) constructed as nonholonomic deformations of prime generating functions (17) with nonlinear superpositions of (18) and (19) reproduce, for explicit parameterizations of constants, various type of nonlinear interactions and pattern forming configurations. In condensed matter physics, such interactions between two waves of one wavelength with a third wave of the other wavelength are known both experimentally and theoretically. Both in condensed matter physics and at cosmological scales, that they play a key role in producing a rich variety of interesting phenomena such as web structures [52]. In this and partner works [4,6,7], we show that one can be reproduced also quasipatterns, superlattice patterns, and STC. The geometric methods and mechanism we elaborate in our works can be applied to any systems in which such 3WIs can occur at microscopic or cosmological scales. In condensed matter physics, such effects are confirmed via Faraday wave experiments, for coupled Turing systems, and some optical systems, see details in [44] and references therein. At cosmological scales, similar configurations are observed and modeled theoretically [46,52,53].

Quasicrystal like configurations in MGTs and GR
We can consider other types of quasi-periodic, or aperiodic, generating functions (not) related to 3WIs.
Formation of quasicrystalline structures and analogous dynamic phase field crystal models: Quasicrystal, QC, structures can be modeled by generating functions QC q = q(x i , t) defined as solutions of an evolution equation with conserved dynamics, In this formula, the canonically nonholonomically deformed hypersurface Laplace operator b ∆ : Such operators can be computed for any family of hypersurfaces Ξ t using formulas (10). The functional F in (21) defines an effective free energy where b = det |bìj|, δy 3 = e 3 and the operator Θ will be defined below. Such nonlinear interactions can be stabilized by the cubic term when the second order resonant interactions are varied by setting the value of Q. The average value < q > of the generating function q is conserved for any fixed t. In result, we can consider q as an effective parameter of the system and that we can choose < q > |t=t 0 = 0 since other values can be redefined and accommodated by altering Θ and Q. Using the functional (22), we can elaborate on models of dark energy, DE, for certain locally anisotropic and inhomogeneous cosmological configurations [4,6,7]. Similarly to QCs in condensed matter physics [45], we concluded that the effective free energy F [q] characterizes a 3-d phase gravitational field crystal model, when modulations are generated with two length scales for off-diagonal cosmological structures. This model consists of a nonlinear PDE with conserved nonholonomic dynamics resulting in evolution equation. It describes a time evolution of q over diffusive time scales.

3-d phase field like quasicrystal structures and evolution:
In cosmological theories, there are studied scalar fields potentials V (ϕ) modified by effective quasicrystal structures, ϕ → ϕ = ϕ 0 +ψ, where ψ(x i , y 3 , t) with (quasi) crystal like phases described by periodic or quasiperiodic modulations. Such modifications can be modelled in dynamical phase field crystal, PFC, like form [54,45]. Applying such mathematical methods in modern cosmology [4,6,7], we can elaborate models of 3-d nonrelativistic dynamics which determined by Laplace like operators 3 △ = ( 3 ∇) 2 , or b ∆ (the left label 3 emphasizes that such an operator is for a 3-d hypersurface). We write ψ instead of q in order to distinguish such QC structures (which can be generated both by gravitational and matter field with two length scales) from the class of models considered above.
In N-adapted frames with 3+1 splitting the equations for a local minimum conserving dynamics, with two length scales l i = 2π/k i , for i = 1, 2. We can elaborate on local diffusion processes determined by a free energy functional where | 3 g | is the determinant of the 3-d space metric and ǫ is a constant. For simplicity, we can restrict our constructions to only non-relativistic diffusion processes, see [22] for relativistic and N-adapted generalizations.
The functional b F [ψ] is defined by a nonholonomic deformation of the Laplace operator, 3 △ → b ∆, resulting in a nonholonomic distortion of F [ψ].

Solitonic space like distributions and nonlinear waves
Off-diagonal interactions determined by generating and integration functions and nontrivial effective sources in MGTs and GR heterotic string gravity may result in various effects with solitonic like distributions, cosmological and geometric evolution models [21,2,3,4,6,7].

Stationary solitonic distributions:
We shall use distributions sd q = q(r, ϑ, ϕ) as solutions of a respective class of solitonic 3-d equations for ǫ = ±1. The left label sd states that such a function is defined as a "solitonic distribution" when in Nadapted frames a function sd q does not depend on time coordinate. The equations (25) and their solutions can be redefined via frame/coordinate transforms for stationary generating functions parameterized in non-spherical coordinates, sd q = q(x i , y 3 ).
Generating nonlinear solitonic waves: 3-d solitonic waves with explicit dependence on time coordinate t are solutions of such nonlinear PDEs: Applying general frame/coordinate transforms, solitonic waves of type sw q = q(x i , t), = q(x 1 , y 3 , t), or = q(x 2 , y 3 , t), can be used as generating functions for certain classes of nonholonomic deformations of stationary, or cosmological metrics, and as generating sources.

Effective sources with effective quasiperiodic free energy
We can prescribe respective generating functions Φ(r, θ, ϕ), or Ψ(r, θ, ϕ), (for stationary configurations with nonlinear symmetry (A.17)), and Φ(x i , t), or Ψ(x i , t), (for cosmological solutions with nonlinear symmetry (A.32)) for quasiperiodic and/or aperiodic effective sources (15) in (14). Such configurations can be determined by additive source functionals and effective cosmological constants, or by nonlinear functionals.

Additive effective sources and cosmological constants
We shall be able to integrate in explicit form gravitational and matter fields systems of nonlinear PDEs for parameterizations of N-adapted sources (A.6) as functionals of type In these formulas, the left labels emphasize what types of effective v-sources are considered. For simplicity, the h-sources can be taken any general ones h Υ(r, θ), or h Υ(x i ). The functional dependence [...] is parameterized for such classes of functions: P q(x i , t) (20) with stationary configurations P q 0 := P q(x i , t = t 0 ) for a fixed value t 0 ; QC q = q(x i , t) (21) with stationary configurations QC q 0 := QC q(x i , t 0 ); ψ(x i , y 3 , t) (23) with stationary ψ 0 := ψ(x i , y 3 , t 0 ); sd q(r, ϑ, ϕ) (25) when cosmological configurations are generated by any source sd q 0 := sd q(r, ϑ, ϕ 0 ) for any fixed value ϕ = ϕ 0 ; and sw q = q(x i , t) (25) with stationary sw q 0 := sw q(x i , t 0 ).
For additive sources (27) and (28), respective stationary and cosmological configurations posses nonlinear symmetries: Such formulas are of type (A.17) and (A.32) and can be stated separately for all sources. The generating functions are chosen in a general form Ψ, or Φ, and, correspondingly, Ψ, or Φ. The QC like components of such quasiperiodic/ aperiodic structures are characterized by free energy functionals (22), for QC q, and (24), for ψ.

Nonlinear functionals for effective sources and cosmological constants
The modified Einstein equations can be integrated in explicit form for general nonlinear functionals The respective nonlinear symmetries (A.17) and/or (A.32) are parameterized resulting in functional dependencies of effective cosmological constants. The formulas (31) and (32) transform respectively in (27) and (28) for additional effective sources and cosmological constants. We emphasize that nonlinear effects are very important in structure formation and for multi-wave nonlinear (solitonic or other types) matter fields and gravitational interactions. Nonlinear dependencies and running of physical constants can be considered also in quantum models. Such stationary and/or cosmological solutions (without contributions of quasiperiodic fields) were studied in a series of our works [1,13,24,26], see also recent results for quasiperiodic structures [4,5,6,7].
For nonlinear effective sources (31) and respective nonlinear symmetries (33), we can define additionally free energy functionals (22), for QC q, and (24), for ψ. Such values and constants have to be determined in explicit form for astrophysical and/or cosmological configurations in order to describe observational data for dark matter and dark energy distributions with respective scales and/or quasiperiodic/ aperiodic structures.

Quasiperiodic generating functions
We can generate various types of gravitational field stationary and/or cosmological configurations using respective classes of generating functions. Such quasiperiodic/aperiodic configurations can be defined by additive quadratic functionals, or in some general nonlinear forms.

Additive quasiperiodic quadratic generating functions
We can prescribe a nontrivial cosmological constant Λ, or Λ, and consider generating functions of type where the left label "a" emphasizes that we certain additions of functionals. Nonlinear symmetries of type type (A.17) and/or (A.32) allow to compute respectively corresponding data ( a Ψ, Υ), or ( a Ψ, Υ), for certain fixed effective sources Υ = Υ(r, θ, ϕ), or Υ = Υ(x i , t). The respective formulas are The QC like components of such quasiperiodic/ aperiodic gravitational structures are also characterized by respective effective free energy functionals (22), for QC q, and (24), for ψ encoding nontrivial vacuum structures with effective cosmological constant.

Nonlinear functionals for quasiperiodic quadratic generating functions
We shall be able to generate in explicit form solutions of modified Einstein equations for nonlinear functionals Respective nonlinear symmetries (A.17) and/or (A.32) involve correspondingly an effective cosmological constant, Λ , or Λ, and nonholonomic constraints for v-sources, Υ(r, θ, ϕ), or Υ(x i , t), resulting in functional dependencies of effective cosmological constants. The formulas (39)- (42) can be re-parameterized respectively as (35)- (38) when the values of effective cosmological constant and matter sources are prescribed to be compatible with experimental data. For such configurations, we can consider structures described additionally by a free energy functional (22), for QC q, modeling QC like gravitational nonholonomic deformations.

(Non)stationary Black Hole Deformations and Quasiperiodic Structures
Generic off-diagonal nonholonomic deformations of BH like solutions were constructed in [8,9,11,12,13,23,25,2,3] by applying the AFDM in various theories of (non) commutative, generalized Finsler, supergravity and supersting MGTs and GR. For small parametric decompositions, such solutions define black ellipsoid stationary configurations, deformations of BH horizons and locally anisotropic polarizations of physical constants, deformations of vacuum solutions into nonvacuum ones. Different classes of solutions were generated with nontrivial (non)commutative backgrounds, containing solitonic distributions and/or describing propagation of black holes in extra dimensions, geometric flows of black holes etc. The approach was generalized for various quasiperiodic, quasicrystal and other type aperiodic solutions in (super) string gravity and for nonholonomic Ricci soliton configurations [4,5].
The goal of this section is to construct and study physical implications of (non) stationary generic offdiagonal solutions describing deformations of some prime BH solutions by quasiperiodic / aperiodic structures. The necessary geometric formalism is summarized in Table 2 and appendices A.2.1 and A.3.1. In this work, the solutions are generated in general forms (not depending on small parameters, see footnote 5).

Nonlinear PDEs for quasiperiodic / aperiodic stationary configurations
There are two possibilities to transform the (modified) Einstein equations (14) into systems of nonlinear PDEs (A.7)-(A.10) with quasiperiodic solutions. In the first case, one considers quasiperiodic sources determined by some additive or general nonlinear functionals. In the second case, additive/ general nonlinear functionals on quasiperiodic solutions are prescribed for generating functions. It is possible also to construct certain classes of solutions involving nonlinear functionals both for generating functions and (effective) sources.

Gravitational eqs and nonsingular solutions for stationary quasiperiodic sources
Stationary solutions with additive sources: Considering a source of type (27) (the left label as is used for "additive stationary"), when the equation (A.14) transforms into ̟ ⋄ h ⋄ 4 = 2h 3 h 4 as Υ. This equation can be integrated on y 3 = ϕ. Exact solutions for stationary configurations of the systems of nonlinear PDEs (A.13)-(A.16) can be constructed following the procedure summarized in Table 2. We can generate such off-diagonal metrics and generalized connections in general form for a generating function h 4 (r, θ, ϕ) with Killing symmetry on ∂ t determined by sources ( h Υ, as Υ) and effective cosmological constant as Λ := P 0 Λ + QC 0 Λ + ψ 0 Λ + sd Λ + sw 0 Λ related to as Υ via nonlinear symmetry transforms (29). The corresponding class of quadratic elements defining stationary solutions can be written in the form . Such solutions are, in general, with nontrivial nonholonomically induced torsion (A.2). They can be re-defined equivalently in terms of generating functions Ψ(r, θ, ϕ) or Φ(r, θ, ϕ), see (A. 19).
LC-configurations in GR determined by quaisperiodic sources can be extracted for additional zero torsion constraints resulting in a more special class of "integrable" generating functions (ȟ 4 , andΨ(r, θ, ϕ) and/oř Φ(r, θ, ϕ)) for respective sources asΥ and as Λ (A.20), Above classes of solutions define stationary off-diagonal gravitational solutions generated by quasiperiodic/aperiodic additive sources as Υ. The term ψ 0 Υ can be used for standard and/or dark matter fields but other ones ( P 0 Υ, QC 0 Υ, sd Υ, sw 0 Υ) may model dark matter stationary distributions with respective quasiperiodic/ aperiodic / solitonic configurations. Considering smooth classes of generating / integration functions and sources, we can construct various classes of nonsingular exact solutions. Applying similar methods, we can generate, for instance, generalizations of stationary models with nonlinear diffusion, fractional, self-organizing and other type processes Refs. [21,23].

Stationary solutions with nonlinear functional sources:
We can work with general nonlinear quasiperiodic / aperiodic / soliton functionals for effective source of type qp Υ(r, θ, ϕ) = qp Υ[ P q 0 , QC q 0 , ψ 0 , sd q, sw q 0 ] (31) with nonlinear symmetries (33). Stationary solutions of the nonlinear system PDEs (A.13)-(A.16) can be written as in Table 2, This formula is similar to (44) but with another type of nonlinear generation functions for (effective) sources for dark and/or usual matter sources, when as Υ → qp Υ and as Λ → qp Λ. Similar re-definitions of additive sources and cosmological constants in (45) into nonlinear functionals generate nonlinear LC-configurations We note that formulas (46) and (47) provide respective generalizations of some classes of solutions (44) and (45) considering general "functionals of functionals" with quasiperiodic (effective) real and dark matter structures which can be organized by corresponding parameters in certain forms with cosmic webs, filaments, quasiperiodic/ aperiodic and or/ solitonic distributions etc.

Nonsingular solutions for stationary quasiperiodic generating functions
Stationary solutions with additive generating functions: For this class of solutions, the quasiperiodic / aperiodic / solitonic structure is stated via generating functions (35) for nonlinear gravitational field interactions without similar prescriptions for (effective) sources as in the previous subsection. We write in brief for additive generating functions subjected to nonlinear symmetries (37). The equation (A.14) transforms into a functional equation For LC-configurations, we obtain (see the third parametrization in (A.20)) We emphasize that the (effective) source Υ in formulas (49) and (50) is not obligatory quasiperiodic. Such a source is considered for general matter fields (including both types of standard and dark matter) and stationary distributions described by an effective cosmological constant Λ. In another turn, the gravitational field distributions are with quasiperiodic / aperiodic / solitonic structure. For such classes of solutions, the gravitational fields encode certain dark energy nonlinear distributions with rich nonholonomic geometric structure and generalized nonlinear symmetries. The coefficients of this class of d-metrics can be chosen to be of necessary smooth class (for instance, nonsingular ones).

Stationary solutions with nonlinear quasiperiodic functionals for generating functions:
Above classes of generic off-diagonal solutions can be generalized for nonlinear quasiperiodic generating functionals qp Φ 2 [ P q 0 , QC q 0 , ψ 0 , sd q, sw q 0 ] (39) characterized by nonlinear symmetries of type (41 We can impose additional zero torsion constraints and extract LC-configurations as in (A.20), Considering additional assumptions and approximations for additive functionals, the formulas (51) and (52) transform respectively into (49) and (50). For such classes of solutions, the gravitational fields encode certain dark energy nonlinear distributions with a more rich nonholonomic geometric structure and generalized nonlinear symmetries when quasiperiodicity is induced from "quasiperiodicity" of matter fields. The coefficients of this class of d-metrics can be chosen to be of necessary smooth class (for instance, nonsingular ones) but can involve certain stochastic sources and fractional derivative processes.

Stationary solutions from nonlinear functionals for quasiperiodic coefficients & sources
In a more general context, we can generate nonsingural stationary off-diagonal generalized quasiperiodic solutions of the (modified) Einstein equations determined both by nonlinear functionals for generating functions, qp Φ (39) and nonlinear functionals for (effective) sources, qp Υ (31). The quasiperiodic data (for instance, scales, interaction constants and associated free energies) for the generating functions are different from the quasiperiodic data for sources. Nevertheless, such data can not be arbitrary independent ones but subjected to nonlinear symmetries generalizing (41) and (33), qp Λ qp Ψ 2 = qp Φ 2 | qp Υ| − dϕ qp Φ 2 | qp Υ| ⋄ . For additive functionals both in the gravitational and (effective sources), such a nonlinear symmetry transforms into a Λ a Ψ 2 = a Φ 2 | a Υ| − dϕ a Φ 2 | a Υ| ⋄ , which is a generalization of (37) and (29).
Following again the procedure summarized in Table 2 but for the data ( qp Ψ, qp Υ), and/or, equivalently, ( qp Φ, qp Λ), the general multi-functional nonlinear generalization of stationary solutions (46) and (51) are constructed in the form For LC-configurations, we obtain multi-functional nonlinear generalizations of (47) and (52), for stationary solutions in GR, The class of solutions (53) describes off-diagonal stationary configurations determined by multi-functional nonlinear quasiperiodic structures both for the dark energy (nonlinear gravitational distributions) and for the dark (and standard) matter fields. In explicit form, such data can be stated to be compatible with observations in modern astrophysics and cosmology.

BHs in (off-) diagonal quasiperiodic media
Various classes of generic off-diagonal stationary quasiperiodic solutions can be described in terms of ηpolarization functions as in appendix A.3.1 and following the geometric method summarized in Tables 1 and  2. In this section, we consider a primary BH d-metricg (8) defined by data [g i (r, θ, ϕ),g a =h a (r, θ, ϕ);N 3 k = w k (r, θ, ϕ),N 4 k =n k (r, θ, ϕ)] which can be diagonalized (for simplicity, we consider the Schwarzschild metric) by frame/ coordinate transforms. The stationary quasiperiodic solutions will be determined by target metrics g generated by nonholonomic deformationsg The quadratic elements corresponding to by g are parameterized in some forms similar to (7), with summation on repeating contracted low-up indices.

Singular solutions generated by stationary quasiperiodic sources
We consider qusiperiodic sources of type qp Υ[ P q 0 , QC q 0 , ψ 0 , sd q, sw q 0 ] as in (46) and compute the coefficients of (55) following formulas.
In these formulas, η 4 (r, θ, ϕ) is taken as a (non) singular generating function. Other types of generating functions are determined with nonlinear symmetries (A.17) and functionals of η 4 (r, θ, ϕ) and data for the prime d-metric, We can constrain the coefficients (56) to a subclass of data generating target stationary off-diagonal metrics of type (A.20) with zero torsion, In above formulas, the nonlinear functionals for the quasiperiodic v-source and (effective) cosmological constant can be changed into additive functionals qp Υ → as Υ and qp Λ → as Λ. The singular behaviour of such solutions is generated by some prime BH data. For certain classes of generating functions and sources, the same type of singularity is preserved. Such examples have been studied in some general forms for small parametric deformations. Similar stationary configurations can be computed for general quasiperiodic structures. The constructions depend on the type of explicit model we construct (for instance, with certain web / filament / solitonic stationary distributions). We can treat such generic off-diagonal stationary solutions as certain conventional nonholonomically deformed BH configurations imbedded into certain (non) singular media determined by stationary dark and usual matter quasiperiodic distributions.

BH solutions deformed by stationary quasiperiodic generating functions
Let us consider another class of solutions (55) when the coefficients of the d-metrics are determined by nonlinear generating functionals qp Φ 2 [ P q 0 , QC q 0 , ψ 0 , sd q, sw q 0 ] (39). Similarly, we can consider additive functionals a Φ 2 (48) including terms with possible itnegration functions h [0] 4 for h 4 and defined by some prescribed data Υ(r, θ, ϕ) and Λ. The formulas for nonlinear symmetries (41) allow us to find (recurrently) corresponding nonlinear functionals, qp η 4 (r, θ, ϕ), or additive functionals, a η 4 (r, θ, ϕ), and related polarization functions, The coefficients of (55) are computed Target stationary off-diagonal metrics (A.20) with zero torsion can be generated by polarization functions subjected to additional integrability conditions, The solutions determined in this subsection describe certain nonholonomically deformed BH configurations self-consistently imbedded into a quasiperiodic gravitational (dark energy) media.

Stationary BH deformations by quasiperiodic sources and generating functions
The most general class of nonholonomic stationary quasiperiodic deformations of BHs is determined by corresponding nonlinear quaisperiodic functionals both for the generating functions and (effective) sources. Nonlinear superpositions of solutions of type (55) and (57) are determined by coefficients of (55) computed ; η 4 = qp η 4 (r, θ, ϕ) as a generating function; (58) In such formulas, there are considered nonlinear generating functionals qp Φ 2 [ P q 0 , QC q 0 , ψ 0 , sd q, sw q 0 ] (39). Similarly, we can consider additive functionals a Φ 2 (48) for some prescribed nonlinear functionals qp Υ(r, θ, ϕ) and qp Λ (in particular, additive nonlinear functionals, as Υ (43) and as Λ can be taken). All such data are related via nonlinear symmetries generalizing (41) which allows to find (recurrently) corresponding nonlinear functionals, qp η 3 (r, θ, ϕ), , or additive functionals, a η 3 (r, θ, ϕ), for the polarization function, Imposing additional conditions for zero torsion, target stationary off-diagonal metrics (A.20) are generated Finally, we conclude that the class of stationary solutions of type (58) describes nonholonomic deformations of a BH self-consistently imbedded into quasiperiodic gravitational (dark energy) backgrounds and quasiperiodic dark/ standard matter. Such distributions can be with a very rich geometric structure (in general, with singular, or nonsingular configurations) and described by respective nonlinear symmetries.

Nonstationary deformations of BH metrics into quasiperiodic cosmological solutions
Prime BH metrics can be deformed nonholonomically into certain classes of exact quasiperiodic solutions in MGTs and GR depending in explicit form on a time like coordinate t. For such configurations, nonlinear quasiperiodic interactions can transform black hole spacetimes into locally anisotropic cosmological ones. In order to construct such exact solutions, we consider a primary BH d-metricg (8)

which via coordinate transforms is parameterized by data
The non-stationary quasiperiodic solutions will be determined by target metrics g generated by nonholonomic deformations The quadratic elements generated by such g are parameterized, with summation on low-up indices.
The quasiperiodic sources are given by nonlinear functionals for effective sources, and nonlinear functionals on cosmological constants qp Λ = Λ[ P Λ, QC Λ, ψ Λ, sd Λ, sw Λ]. In general, we can consider some nonlinear functionals for generating , when the nonlinear symmetries are stated by formulas The equation (A.14) transforms into ̟ * h * For additional zero torsion conditions, we generate target non-stationary off-diagonal metrics (A.35) defined by coefficients Non-stationary solutions generated from BH prime metrics (with zero or non-zero nonholonomically induced torsions) can be constructed in a similar manner for additive functionals for (effective) sources and generating functions which are similar to those for locally anisotropic spacetime. The integration functions and primary data are taken for BH stationary configurations.
Off-diagonal solutions of type (60) describe nonholonomic deformations of a BH solution by imbedding it self-consistently in a dynamical quasiperiodic gravitational (dark energy) background and quasiperiodic dark/ standard matter interacting fields. Such locally anisotropic cosmological like spacetimes contain BH like structures are described by relativistic evolution of geometric structure, singular, or nonsingular, and described by respective nonlinear symmetries. For certain special parameterizations, we can model "dissipation" of singular BH structure into certain types (locally anisotropic) cosmological spacetimes.

Off-diagonal Quasiperiodic Cosmological Spacetimes
Locally anisotropic and inhomogeneous cosmological solutions and accelerating universe scenarios (in MGTs, GR, and geometric flow theories) were studied in a series of works [12,24,26,1,50,51], see also references therein. Various classes of generic off-diagonal cosmological metrics were constructed by applying the AFDM as a geometric alternative to numeric methods [19]. Recently, the approach was developed by constructing quasiperiodic cosmological solutions with small parametric deformations [6,7].
The goal of this section is to study physical implications of (non) stationary generic off-diagonal solutions describing deformations of some prime cosmological spacetimes by quasiperiodic / aperiodic / solitonic and/or pattern forming structures. The necessary geometric formalism is summarized in Table 3 and appendices A.2.2 and A.3.2. We emphasize that in this work the cosmological solutions are constructed for general nonlinear or additive nonlinear functionals for generating functions and (effective) sources without additional assumptions on modelling small parameter configurations (see footnote 5).

Nonlinear PDEs for quasiperiodic/ aperiodic cosmological configurations
There are two possibilities to transform the (modified) Einstein equations (14) into systems of nonlinear PDEs (A.21)-(A.24) with quasiperiodic solutions depending in explicit form on a time like variable. In the first case, one considers quasiperiodic sources determined by some additive or general nonlinear functionals. In the second case, respective nonlinear functionals determining quasiperiodic solutions are prescribed for generating functions. It is also possible to construct certain classes of locally anisotropic and inhomogeneous cosmological solutions by considering nonlinear / additive functionals both for generating functions and (effective) sources.

Cosmological solutions for quasiperiodic sources
Cosmological solutions generated by additive functionals for sources: Let us consider an additive functional for a quasiperiodic source of type Υ(x i , t) (28),  Table 3. Such generic off-diagonal solutions are parameterized in the form For local pseudo-Riemannian configurations, we have to fix respective sign of the coefficient h 4 (x k , t) which can be considered as a generating function with Killing symmetry on ∂ 3 determined by sources ( h Υ, as Υ). Such solutions are of type (A.34) and can be re-written equivalently with coefficients stated as functionals of as Φ and as Ψ define by nonlinear symmetry formulas determined by additive (effective) sources and cosmological constants. We can extract from (62) locally anisotropic generic off-diagonal cosmological LC-configurations in GR determined by quaisperiodic sources by imposing additional zero torsion constraints. This restricts the respective classes of "integrable" generating functions (ȟ 4 ,Ψ(x i , t) and/orΦ(x i , t)) for sources asΥ and as Λ (A.35), Above linear quadratic elements define off-diagonal cosmological solutions generated by quasiperiodic additive sources as Υ. The term ψ Υ encodes contributions of standard and/or dark matter fields. The sources ( P Υ, QC 0 Υ, sd 0 Υ, sw Υ) can be chosen to model dark matter stationary distributions with respective quasiperiodic/ aperiodic / solitonic nonlinear wave interactions. Considering smooth classes of generating / integration functions and sources, we can generate nonsingular cosmological exact solutions. Applying similar methods, we can study, for instance, generalizations of cosmological models to effects determined by nonlinear diffusion, fractional, self-organizing and other type processes as we proved in Refs. [21,23].
Cosmological solutions generated by nonlinear functionals for quasiperiodic sources: Exact solutions can be generated by nonlinear quasiperiodic functionals for effective sources, qp Υ(x i , t) = qp Υ[ P Υ, QC 0 Υ, ψ Υ, sd 0 Υ, sw Υ] (32) with nonlinear symmetries (34). Using the AFDM summarized in Table  3, we construct cosmological solutions of the nonlinear system PDEs (A.28)-(A.31), This formula is similar to (44) but with another type of nonlinear generation functions for (effective) sources for dark and/or usual matter sources, when as Υ → qp Υ and as Λ → qp Λ. Similar re-definitions of additive sources and cosmological constants in (45) into nonlinear functionals generate nonlinear LC-configurations For additive functionals for cosmological sources, the formulas (64) and (65) transforms respectively into quadratic linear elements (62) and (63). Considering small parametric deformations for 4-d cosmological solutions, we can reproduce the results from [6,7].

Cosmological solutions for nonstationary quasiperiodic generating functions
In this subsection, the sources are with arbitrary data Υ µ ν = [ h Υ(x k ), Υ(x k , t)] but the generating functions are considered for some additive or general nonlinear functionals with quasiperiodic structure.

Locally anisotropic and inhomogeneous cosmological metrics with additive generating functions:
For this class of solutions, the quasiperiodic / aperiodic / solitonic structure is stated via generating functions (36) for nonlinear quasiperiodic gravitational field interactions buth without explicit prescriptions on any quasiperiodic structure for (effective) sources. In brief, such additive functionals are written with nonlinear symmetries (38), where sd 0 Φ 2 [ sd q] is taken for any functions sd q(x k ) but other components are considered as functionals on respective functions P q, QC q, ψ, and sw q on (x i , t). In result, the equation For LC-configurations, we obtain (see the third parametrization in (A.35)) (67) We can consider small parametric decompositions and frame/coordinate transforms in order to related such solutions to some well known (off) diagonal cosmological metrics.
Cosmological solutions with nonlinear quasiperiodic functionals for generating functions: The formulas (66) and (67) can be generalized in order to generate solutions from nonlinear functionals for generating functions when nonlinear quasiperiodic generating functionals  Table  3. We obtain such solutions: For zero torsion constraints in order to extract LC-configurations, Here we emphasize that there is certain duality between formulas when y 4 = t ↔ y 3 and, respectively, "overlined" values are changed into "not overlined" ones, and inversely. For such nonholonomic dual transforms, the formulas (51) and (52) transform into corresponding (68) and (69) [inverse maps can be also considered]. The coefficients of this class of d-metrics can be chosen to be of necessary smooth class and involve certain stochastic sources and fractional derivative processes. Such nonholonomic deformation and generalized transform may change the topological spacetime structure and encode dark energy and dark matter effects.

Cosmology from nonstationary functionals for quasiperiodic coefficients & sources
Conventionally, all classes of considered above cosmological solutions can be formulated in terms of generalized quasiperiodic nonlinear functionals both for generating functions, qp Φ (40) and nonlinear functionals for (effective) sources, qp Υ (32). Such data are subjected to conditions of nonlinear symmetries generalizing (42) and (34), Similar nonlinear symmetries can be considered for additive functionals both for the gravitational fields and (effective) sources, when (38) and (30) The procedure summarized in Table 3 and generalized for the data ( qp Ψ, qp Υ), and/or, equivalently, ( qp Φ, qp Λ), allows us to construct multi-functional nonlinear quasiperiodic cosmological configurations, For LC-configurations, we obtain multi-functional nonlinear generalizations of (68) and (69) which result in more general classes of locally anisotropic and inhomogeneous solutions in GR, The classes (70), of cosmological solutions in MGTs, and (71), for cosmological solutions in GR, describe offdiagonal non-stationary configurations determined by multi-functional nonlinear quasiperiodic structures. Such rich geometric nonholonomically dynamical structures can be described both for the dark energy (nonlinear gravitational distributions) and for the dark (and standard) matter fields. In general, such cosmological solutions may not have smooth limits to well-known cosmological metric (for instance, FLRW, or any Bianchi type; we shall study such configurations in next subsections). In explicit form, the geometric data for such generic off-diagonal cosmological solutions can be stated to explain various observations in modern cosmology.

Cosmological metrics evolving in (off-) diagonal quasiperiodic media
In this subsection, generic off-diagonal quasiperiodic cosmological solutions are constructed in terms of ηpolarization functions as in appendix A.3.2 and following the geometric method summarized in Tables 1 and 3. We consider a primary cosmological d-metricg (8) which can be diagonalized for a FLRW cosmological metric (in general, we can consider off-diagonal Bianchi anisotropic metrics) by frame/ coordinate transforms. The cosmological quasiperiodic solutions will be determined by target metrics g generated by nonholonomic deformationsg The quadratic elements corresponding to cosmological metrics g are similar to (2) but with a corresponding parametrization in terms of polarization functions, The target d-metrics g = g(x k , t) are characterized by N-adapted coefficients

Cosmological evolutions generated by nonstationary quasiperiodic sources
We consider qusiperiodic cosmological sources of type qp Υ(x i , t) = qp Υ[ P Υ, QC 0 Υ, ψ Υ, sd 0 Υ, sw Υ] (32) as in (64) and (65) and compute the η-polarization functions following formulas There are other types of f generating functions determined determined by η 3 (x i , t), data for the prime d-metric, subjected to nonlinear symmetries (A.32), Restricting the classes of generating functionals and sources to subclasses of solutions of (A.25) (see Table  3 and appendix A.3.2) for target off-diagonal cosmological metrics (A.35) with zero torsion, we obtain η i = e ψ /g i ; η 3 =η 3 (x i , t) as a generating function; In (72) and next formula for LC-configurations, the nonlinear functionals for the quasiperiodic v-source and (effective) cosmological constant can be changed into additive functionals qp Υ → as Υ and qp Λ → as Λ. We generate solutions of type (62) and, respectively, (63) but with a prime cosmological structure modified by certain classes of quasiperiodic generating functions and sources. For cosmological scenarios, we can generate solutions without singularities. Such examples have been studied in details for small parametric deformations in [12,24,26,1,50,51] and, for quasiperiodic configurations, in [6,7]. In this work, cosmological configurations are studied for general quasiperiodic structures. The constructions depend on the type of model we construct (for instance, different classes of generating functionals and sources have to be prescribed for certain web / filament / solitonic stationary distributions). We can treat such generic off-diagonal locally anisotropic solutions as nonholonomically deformed prime cosmological configurations imbedded into certain (non) singular media determined by evolving dark and usual matter with quasiperiodic interactions.

Cosmology from nonstationary quasiperiodic generating functions
We can construct other classes of 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 qp Φ 2 [ P q, QC q, ψ, sd 0 q, sw q] (40). In a similar manner, we can generate similar cosmological metrics by additive functionals a Φ 2 (35) for a prescribed effective source Υ(x i , t) and cosmological constant Λ. The formulas for nonlinear symmetries (42) allow us to find (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 (A.34) are recurrently computed, Target off-diagonal cosmological metrics (A.35) with zero torsion can be generated by polarization functions see appendix A.3.2 for explanations on conventions and nonholonomic constraints for functions with "inverse hats" and transforming general quasiperiodic functionals into additive ones. The solutions generated in this subsection describe certain nonholonomically deformed prime cosmological configurations (for instance, a FLRW, or Bianchi, type metric) self-consistently imbedded into a quasiperiodic gravitational (dark energy) media.

Cosmological solutions for nonstationary quasiperiodic sources & generating functions
We can construct very general classes of nonholonomic deformations of prime cosmological metrics generated by nonlinear quaisperiodic functionals both for the generating functions and (effective) sources. Nonlinear superpositions of cosmological solutions (72) and (73) are determined by coefficients of (64) computed, In such formulas, there are considered nonlinear generating functionals characterized by nonlinear symmetries of type (42) for some prescribed nonlinear functionals qp Υ(x i , t) and qp Λ. Instead of ( qp Φ, qp Υ, qp Λ), we can consider additive data ( a Φ, a Υ, a Λ). Using formulas for nonlinear symmetries, we can define general nonlinear, qp η 3 (x i , t), or additive functionals, a η 3 (x i , t), for the polarization function, LC-configurations with zero torsion for target off-diagonal cosmological metrics (A.35) are generated We note that the class of cosmological solutions (74) is "dual" on y 3 and y 4 coordinates to stationary solutions (58). For certain classes of parameterizations, such classes describe nonholonomic deformations of cosmological spacetimes, or BHs, self-consistently imbedded into quasiperiodic cosmological (dark energy) backgrounds and quasiperiodic dark/ standard matter.

Discussion and Concluding Remarks
Cosmology theories rely on the cosmological principle stating that our universe is sufficiently homogeneous and isotropic on large scales. Geometric models and computation algorithms in numerical relativity encode observational data and theoretical assumptions that expansion is that for the FLWR model governed by (modified) Friedmann equations and employing a Newtonian approximation of gravity. The transition to cosmic homogeneity begins on scales ∼ 80h −1 Mpc when the Universe is inhomogeneous and anisotropic on smaller scales, see Planck2015 data [56,57,58,59]. Modern telescopes reached a precision which shows that nonlinear general relativistic effects from inhomogeneities could be important.
There are more extreme hypotheses that inhomogeneities may provide an alternative explanation for the accelerating expansion of the Universe (for instance, it is considered the back reaction or replaced the role assigned to dark energy in the standard ΛCDM model). Such alternative approaches are based on a number of cosmological observations during the last 20 years emphasizing important phenomena of the accelerating Universe and dark energy and dark matter. The dark energy physics yields a late time acceleration of the spacetime. In its turn, the dark matter physics is "hidden" as an invisible matter (there are various models of dust, cold and hot matter, gravitational polarizations etc.) which favour the process of gravitational clustering. It is not clear how such intriguing physical effects (substances and (non) linear interactions) could emerge in the general relativity theory, GR, or should complement such a theory. The search for alternative modified gravity theories, MGTs, has become an active area for theoretical and experimental investigations. A very important task is to elaborate on new methods of constructing exact and parametric solutions of motion and evolution field equations in MGT and GR describing nonlinear gravitational and matter field interactions (as we motivated in the Introduction section). A large number of MGTs is available and studied in various details in modern literature [47,60,61,62,63,64,65,66,67,68,69,70,24,1,2,3,25,26].
The predominantly accepted ΛCDM paradigm predicts for the cosmological structure that dark matter is organized as a cosmic web structure of walls, filaments and halos [46,47,48,49,50,51,52,53,54,55]. Numerous images of the visually striking cosmic web have been created and studied theoretically in a framework of research of the large-scale structure of the universe. Recently, an advanced observational and numerical simulation technique was elaborated for pseudo-3d visualisations aimed to minimal loss of information and accurate representation of cosmic web fundamental shapes and components. In the cosmic web, one observes large scale filaments with lengths that can reach tens of Mps, see [55] and references therein. Such filaments are observed indirectly through the galaxies positions (using large galaxy surveys), or through absorption features in the spectra of hight red shift sources (with direct detecting of intergalactic medium filaments through their emission on the HI 21cm line). One estimates that gas in filaments of length l ≥ 15h −1 Mpc with relatively small inclinations to the line of sight (≤ 10 • ) can be observed during 40-100 hours with modern (radio) telescopes. It is revealed from observations of the local Universe that galaxies reside in a complex network of filamentary structures (with cosmic web). Such structures can be explained in the ΛCDM framework as resulting from nonlinear gravitational evolution. The dark matter halos (within which galaxies reside) are connected to each other through a patchwork of filaments and sheets and quasiperiodic/ aperiodic structures that constitute the structure of the intergalactic medium, IGM. In modern astronomy it is explored the possibility of using the HI 21 cm line to directly observe ICM filaments.
Various alternative approaches to structure formation in modern cosmology were elaborated using numerical relativity [15,16,17,18]. This involves modelling by numerical methods which began with evolutions of planar and spherically symmetric spacetime following the Arnowitt-Deser-Misner, ADM, formalism [27]. Latter, there were considered generalizations with Kasner and matter-fields, and for the propagation and collision and gravitational wave perturbations and linearised perturbations to homogeneous spacetimes. In order to simplify the numerical calculations there are included linear and nonlinear symmetries. Simulations with (non) linear symmetries have been performed in order to study the evolution of small perturbations to an FLRW spacetime, to explore observational implications showed differential expansion in an inhomogeneous universe. Such works indicate that the effects of nonlinear inhomogeneities may be significant.
Modern cosmological observational data and related phenomenological models emphasize a crucial importance of nonlinear physics and related mathematical methods. Distinguishing general relativistic effects determined by nonlinear structures requires more general classes of solutions of Einstein's equations in GR and application of advanced geometric and numeric methods elaborated recently for theoretical studies in MGT, non standard particle physics, astrophysics and cosmology. Post-Newtonian, small parametric (perturbative and non-perturbative) numeric approximations consist a worthwhile approach extended with methods including density perturbations in a highly nonlinear form. However the validity of nonlinear effects must be checked against more precise exact and parametric solutions constructed in analytic form.
In this work, we performed a brief review and feasibility study of geometric methods of constructing exact offdiagonal solutions to the (modified) Einstein equations for stationary, locally anisotropic BH and cosmological (in)homogeneous cosmology by geometric modeling and comparing with numeric simulating the growths of quasiperiodic / aperiodic structures and comparing to known analytic solutions. We also presented a study on the evolution of nonlinear stationary and/ or locally anisotropic configurations and analysed the resulting new classes of exact solutions in MGTs and GR. Such a research was performed by developing the anholonomic frame method, AFDM, for constructing new classes of stationary and nonstationary (cosmological) solutions of (modified) Einstein equations. Such solutions are described by generic off-diagonal metrics and generalized connections and depend, in principle, on all possible 4-d and extra dimensional space coordinates; see details, examples and various applications in [8,9,10,11,12,13,14,20,21,22,23,24,25,26].
A crucial difference from former approaches to constructing exact solutions [27,28,29,30,31] is that the AFDM allows us to work with generating and integration functions for coefficients of generic off-diagonal metrics, generalized connections and (effective) sources transforming motion and geometric evolution equations into nonlinear systems of partial differential equations, PDEs, with decoupling properties. Following this geometric method, we perform such nonholonomic deformations of some prime stationary and nonstationary solutions (for instance, BH and/or (an) isotropic cosmological solutions ) when the (generalized) Einstein equations can be decoupled in general forms and integrated for various classes of metrics g αβ (x i , y 3 , t). In particular, we can reproduce former results for diagonalizable ansatz with dependence on radial and warping coordinates as solutions of ordinary differential equations, ODEs. Nevertheless, the AFDM, is more than a constructive interference of geometric and analytic methods for constructing exact solutions for certain classes of important nonlinear systems of PDEs, in mathematical relativity and cosmology. It reflects new and formerly unknown properties and nonlinear symmetries of the (modified) Einstein equations when generic off-diagonal interactions and mixed continuous and discrete structures (quasiperiodic / aperiodic / pattern forming / solitonic ones) are considered for vacuum and non-vacuum gravitational interactions. The AFDM is appealing in some sense it is "economical and very efficient"; allowing us to proceed in the same manner but with fractional / random / noncommutative sources and their respective interaction parameters.
Finally, we emphasize that there are a number of directions in (modified) gravity, cosmology and astrophysics which can be pursued using as starting points the methods and solutions elaborated in this works. 9 This includes noncommutative and nonassociative generalizations defined by possible modified dispersion relations and/ or extra dimension contributions to dark energy/matter physics and/or quantum models with quasiperiodic and pattern forming structures.

A.1 N-adapted coefficients for curvatures and torsions
Using standard formulas, we can define and compute both in abstract and coordinate forms the torsion, T , the nonmetricity, Q, and the curvature, R, tensors for any d-connection D = (hD, vD), In literature [8,9,10,11,12,13,14], there are used terms like distinguished tensor, d-tensor, and distinguished (geometric) object, d-object, (also d-metric, d-connection) etc. for geometric and physical values defined in N-adapted form, i.e. when all values are defined in some coefficient forms preserving under parallelism the N-connection splitting (6)). Using N-adapted coefficients of a d-connection, D = {Γ γ αβ = (L i jk , L a bk , C i jc , C a bc )}, we can compute with respect to N-adapted frames (4) and (3) corresponding N-adapted coefficients We omit such cumbersome formulas which can be found in above mentioned references. The coefficients of the canonical d-connection D = { Γ γ αβ = ( L i jk , L a bk , C i jc , C a bc )} in (9) are L i jk = 1 2 g ir (e k g jr + e j g kr − e r g jk ) , C a bc = 1 2 g ad (e c g bd + e b g cd − e d g bc ) , (A.1) The coefficients of the distortion d-tensor, Z γ αβ = Γ γ αβ − Γ γ αβ (10) can be written in N-adapted form using (A.1) and LC-connection ∇ = {Γ γ αβ }, all computed with respect to (4) and (3). Using such values, we find the nontrivial d-torsion coefficients T γ αβ of D, We note that the d-torsion coefficients (A.2) vanish if in N-adapted form there are satisfied the conditions Following similar formulas, we can compute (see details in above references) the nontrivial coefficients of the Riemann d-tensor, R α βγδ , the Ricci d-tensor, R αβ (11), and the Einstein d-tensor E αβ := R αβ − 1 2 g αβ R.
A.2 Decoupling property of (modified) Einstein equations For general assumptions and using frame/ coordinate transforms, any d-metric g (2) can be parameterized g i = e ψ(r,θ) , g a = ω(r, θ, y b )h a (r, θ, ϕ), N 3 i = w i (r, θ, ϕ), N 4 i = n i (r, θ, ϕ), for ω = 1, stationary conf.; (A.4) In order to write certain formulas in compact forms, we shall use also brief notations of partial derivatives ∂ α q = ∂q/∂u α (for instance, for a function q(x k , y a )) The sources (15) for (effective) matter field configurations can be parameterized via frame transforms in respective N-adapted forms In these formulas, there are considered necessary type vielbein transforms e µ µ ′ (u γ ) and their duals e ν ′ ν (u γ ), , are considered as generating functions for (effective) matter sources imposing nonholonomic frame constraints on stationary distributions or cosmological dynamics of (effective) matter fields. For simplicity, we shall generate in explicit form certain classes of generic off-diagonal solutions with ω = 1 (i.e. with at least one Killing symmetry on ∂ t or ∂ ϕ ) when the frame/coordinate systems and transforms are compatible with the conditions ∂ b h a = 0 and [ h Υ(r, θ), Υ(r, θ, ϕ)] = 0, or [ h Υ(x i ), Υ(x i , t)] = 0. In next subsections, we shall prove that above introduced parameterizations of d-metrics and (effective) sources allows us to integrate explicitly the gravitational field equations (14). 10

A.2.1 Off-diagonal stationary configurations
In this subsection, we outline key steps for proofs of general decoupling and integrability of (modified) Einstein equations for general assumptions on coefficients d-metrics and N-connections which do not depend on y 4 with respect to a class of N-adapted frames.

Nontrivial components of the Ricci d-tensor and (modified) Einstein equations:
Introducing d-metric data (A.4) with ω = 1 into formulas (A.1) and (11) (with respective sources (A.6), and considering nontrivial N-adaped coefficients of the Ricci d-tensor), we transform the modified Einstein equations (14) into such a system of nonlinear PDEs In N-adapted frames, this system posses a decoupling property: The equations (A.7) allow us to find g 1 (or, inversely, g 2 ) for any prescribed h-source h Υ(r, θ) and given coefficient g 2 (or, inversely, g 1 ). Integrating on variable y 3 in (A.8), we can define h 3 (r, θ, ϕ) as a solution of first order PDE for any prescribed v-source Υ(r, θ, ϕ) and given coefficient h 4 (r, θ, ϕ) [we can define h 4 (r, θ, ϕ) if, inversely, h 3 (r, θ, ϕ) is given but in such cases we have to solve a second order PDE]. For well-defined values of h 3 and h 4 , the equations (A.9) transform into an algebraic linear equation for w k (r, θ, ϕ). We have to integrate two times on y 3 in order to compute n k (r, θ, ϕ) for any well-defined h 3 and h 4 . So, the decoupling property of the system (A.7)-(A.10) reflects the possibility to integrate such PDEs step by step by defining the h-coefficients, g i , and v-coefficients, h a , of a d-metric [g i = g i (x k ), g a = g a (x k , y 3 )] (2) (with data (A.4)), and, finally, the N-connection coefficients,

Extracting torsionless configurations:
We note that the LC-conditions (A.3) for stationary configurations transform into a system of 1st order PDEs, imposing additional constraints on off-diagonal coefficients of metrics of type (1). Such conditions can be imposed on d-metric and N-connection coefficients after a class of off-diagonal solutions for (14) has been constructed in an explicit form. For certain well-defined N-adapted parameterizations, the equations (A.11) can be solved in explicit forms.
Nonlinear gravitational PDEs with explicit decoupling: Let us show how the system of nonlinear PDEs (A.7)-(A.10) can be integrated in explicit form. We introduce the coefficients (A.12) Using such nontrivial and nonsingular values for ∂ 3 h a = 0 and ∂ t ̟ = 0, 11 we obtain This system can be integrated in explicit form (see below) for any generating function Ψ(r, θ, ϕ) := e ̟ and generating sources h Υ(r, θ) and Υ(r, θ, ϕ).
The value of Λ can be chosen from certain physical considerations. It can be positive or negative. Solutions with Λ = 0 have to be studied by special methods, see detils and examples in [10,11,12,13,14,20,21,22,23,24,25,26]. We can describe nonlinear systems of PDEs by two equivalent generating data (Ψ, Υ) or (Φ, Λ) [for different classes of solutions, one can be convenient to work with different types of such data]. Modules in such formulas are taken in certain forms which allows to work with real functions.

A.2.2 Off-diagonal cosmological solutions
In this paper, for simplicity, there are considered solutions g αβ (x i , t) with Killing symmetry on ∂ 3 , i.e. with ω = 1 in (A.5), which allows us to generate exact solutions in explicit form. Solutions depending on all spacetime coordinates, g αβ (x i , y 3 , t), can be constructed for nontrivial vertical conformal factors ω(x i , y a ), see details and examples in [10,11,12,13,14,20,21,22,23,24,25,26]. It should be noted that if certain classes of off-diagonal solutions for such nonholonomic cosmological configurations have been constructed in explicit form, we can impose additional nonholonomic constraints, or limits with necessary smooth classes of functions, when g αβ (x i , t) ≈ g αβ (t) are related to Bianchi type, or FLRW, like cosmological metrics.
Nontrivial components of the Ricci d-tensor for nonholonomic cosmological configurations: Let us consider d-metric data (A.5) with ω = 1 in order to compute the N-adapted and nontrivial coefficients D = { Γ γ αβ } (A.1) and R αβ (11). For nontrivial sources [ h Υ(x i )δ i j , Υ(x i , t)δ a b ] (A.6), we obtain from the modified Einstein equations (14) a system of nonlinear PDEs 12 For ∂ t h a = 0 and ∂ t ̟ = 0, 13 we rewrite the equations in the form We can integrate this system for any generating function Ψ(x i , t) := e ̟ and sources h Υ(x i ) and Υ(x k , t).
Nonlinear symmetries for generating functions and sources with effective cosmological constant: The system of two equations (A.27) and (A.29) relates four functions (h 3 , h 4 , Υ, and Ψ) which emphasizes an important nonlinear symmetry for locally anisotropic cosmological solutions and respective generating functions, (Ψ, Υ) ⇐⇒ (Φ, Λ), when This allows to introduce a new generating function Φ(x i , t) and an (effective) cosmological constant Λ = 0, which can be applyied for generating exact off-diagonal solutions in explicit forms.
Above formulas can be simplified by choosing the integration functions in the form h [0] 3 =h 3 and 1 n k =n k = ∂ k n. We use symbols with inverse "hat" following the conventions on integrability of generating functions which are similar to proofs of formula (A.20).