Off-diagonal deformations of Kerr metrics and black ellipsoids in heterotic supergravity

Geometric methods for constructing exact solutions of equations of motion with first order α′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha ^{\prime }$$\end{document} corrections to the heterotic supergravity action implying a nontrivial Yang–Mills sector and six-dimensional, 6-d, almost-Kähler internal spaces are studied. In 10-d spacetimes, general parametrizations for generic off-diagonal metrics, nonlinear and linear connections, and matter sources, when the equations of motion decouple in very general forms are considered. This allows us to construct a variety of exact solutions when the coefficients of fundamental geometric/physical objects depend on all higher-dimensional spacetime coordinates via corresponding classes of generating and integration functions, generalized effective sources and integration constants. Such generalized solutions are determined by generic off-diagonal metrics and nonlinear and/or linear connections; in particular, as configurations which are warped/compactified to lower dimensions and for Levi-Civita connections. The corresponding metrics can have (non-) Killing and/or Lie algebra symmetries and/or describe (1+2)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(1+2)$$\end{document}-d and/or (1+3)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(1+3)$$\end{document}-d domain wall configurations, with possible warping nearly almost-Kähler manifolds, with gravitational and gauge instantons for nonlinear vacuum configurations and effective polarizations of cosmological and interaction constants encoding string gravity effects. A series of examples of exact solutions describing generic off-diagonal supergravity modifications to black hole/ellipsoid and solitonic configurations are provided and analyzed. We prove that it is possible to reproduce the Kerr and other type black solutions in general relativity (with certain types of string corrections) in the 4-d case and to generalize the solutions to non-vacuum configurations in (super-) gravity/string theories.

The supergravity/superstring and MGT gravitational field equations are formulated as sophisticated systems of nonlinear partial differential equations, PDEs. Various types of advanced analytic and numeric methods for constructing exact and approximate solutions of such equations have been explored. For GR, a number of examples of exact and physically important solutions are summarized in monographs [31,32]. In generalized (super-) gravity theories, most of the solutions with a variety of different vacua, such as in string and brane theories, comes from corresponding choices of the internal manifold. Toroidal compactifications, special geometric cases as Calabi-Yau and, more generally but with less supersymmetry, with SU (3) structure manifolds were considered. The bulk of solutions in various super/noncommutative/extra-dimension/modified grav-ity theories are described by metrics, frames, and connections with coefficients of fundamental geometric/physical objects depending on one and/or two coordinates in 4-d to 10-d spacetimes. For well-known classes of solutions, diagonalizations of metrics are possible via coordinate transforms and the linear connection structures are mostly of Levi Civita, LC, and Kähler type. Additional distortions by torsion structures are also considered. Various generalizations of well-known and physically important exact solutions for the Schwarzschild, Kerr, Friedman-Lemaître-Robertson-Walker (FLRW), wormhole spacetimes etc were constructed. These classes of diagonalizable metrics (the off-diagonal terms in the Kerr solutions are determined by rotations and respective frame/coordinate systems) are generated by a certain ansatz where equations of motion are transformed into certain systems of nonlinear second order ordinary equations (ODEs), 2-d solitonic equations etc. These systems of ODEs have Killing vector symmetries which result in additional parametric symmetries [33][34][35] and depend on integration constants.
A number of physically important solutions with black hole, wormhole, cosmological, monopole, and instanton configurations etc. were constructed in different (super-) string/gravity and MGTs for a diagonalizable ansatz depending generically on one spacetime coordinate (in certain cases, with dependencies on two coordinates and with Killing or other type symmetries). For the majority of such solutions, the corresponding equations of motion transform with the corresponding diagonal ansatz for metrics and frame transforms from nonlinear systems of PDE in nonlinear systems of ODEs. The integrals of such ODEs depend on certain integration constants which are defined from some boundary/asymptotic/initial and prescribed symmetry conditions. For many years, physicists and mathematicians concentrated their efforts on constructing further generalizations and applications of "diagonalizable" solutions in string and (super-) gravity theories because it was more "easy" to find analytic and numeric solutions of resulting systems of ODEs. In these approaches, the integration constants can be related to certain physical constants considering the Cauchy problem, or by using various assumptions on asymptotic/boundary/symmetry conditions. All versions of (supersymmetric) modified Einstein equations consist of very sophisticated off-diagonal nonlinear systems of PDEs. In general form, the main properties of nonlinear/nonholonomic/parametric interactions of such (super-) gravitational and matter field systems are described by PDEs and not by approximations to ODEs. To find general classes of solutions in analytical form, understand their geometric properties and search for possible physical implications is of great importance in modern mathematics, physics, and cosmology. Imposing only a "simple" diagonalizable ansatz of higher symmetry (as it is usually considered for generating new classes of exact solutions in various gravity theories), we "cut" the possibility to find infinite numbers of classes of generic off-diagonal solutions determined by generating and integrations functions depending on 3, 4, . . . spacetime coordinates, with various commutative and noncommutative parameters etc from the very beginning. A typical example is that of solitonic wave solutions depending on three variables, which cannot be generated for a very "simple" diagonal ansatz and factorized dependent on certain generated functions. It is possible that various problems in modern acceleration cosmology (with structure formation, dark energy and dark matter, etc.) could be solved by certain (super-) string/gravity solutions related to generic off-diagonal nonlinear configurations in GR and it may not be necessary to radically modify the standard gravity and particle physics theories.
The problem of constructing generic off-diagonal exact solutions 1 with coefficients of metrics and connections and other physically important geometric objects depending on three-, four-and extra-dimensional spacetime coordinates is much more difficult. For instance, there are a maximum of six independent components of a metric tensor from ten components in a 4-d (pseudo-) Riemannian spacetime. 2 Any such ansatz with metrics depending on three-four spacetime coordinates transforms the Einstein equations into systems of nonlinear PDEs, which cannot be decoupled and integrated in a general analytic form if the constructions are performed with respect to local coordinate frames and for the LC-connection. The condition of zero torsion imposes various types of contractions between the coefficients of the linear connection, reference frames and various tensor fields which do not allow any general decoupling of the corresponding systems of PDEs. To generate solutions with generic off-diagonal metrics and generalized connections for higherdimensional configurations (for instance, in 5d-10-d string gravity and MGTs) is a technically more difficult task than in 3d-4d theories.
In a series of publications [35][36][37][38][39][40][41][42][43][44], the so-called anholonomic frame deformation method, AFDM, of constructing exact solutions in commutative and noncommutative (super-) gravity and geometric flow theories has been explored. By straightforward analytic computations, it was proven that it is possible to decouple the gravitational field equations and generate general classes of solutions in various theories of gravity with metric and nonlinear, N-, and linear connections structures. The geometric formalism is based on spacetime fibrations determined by nonholonomic distributions with splitting of dimensions, 2 (or 3) +2 + 2 + · · · . In explicit 1 The metric fields corresponding to such solutions cannot be diagonalized in a finite or infinite spacetime region via coordinate transforms. 2 Four components from a maximum of ten can be fixed to be zero using coordinate transforms, which is related to the Bianchi identities. form, certain classes of N-elongated frames of reference, the considered formal extensions/embeddings of 4-d spacetimes into higher-dimensional spacetimes are introduced and necessary types of adapted linear connections are defined. These connections are called distinguished, d-connections, defined in forms which preserve the N-connection splitting. In Einstein gravity, a d-connection is considered auxiliary, which in certain canonical forms can be uniquely defined by the metric structure following the conditions of metric compatibility and some other geometric conditions (for instance, that certain zero values for "pure" horizontal and vertical components contain nonholonomically induced torsion fields). Surprisingly, such a canonical d-connection allows us to decouple the equations of motion in general forms and generate various classes of exact solutions in generalized/modified string and gravity theories. Having constructed a class of generalized solutions in explicit form (depending on generating and integration functions, generalized effective sources and integration constants), we can constrain the induced torsion fields to zero and "extract" solutions for LC-configurations and/or Einstein gravity. It should be emphasized that it is important to impose the zero-torsion conditions at the end, i.e. after we have found a class of generalized solutions. We cannot decouple and solve the corresponding systems of PDEs in general forms if we use the LC-connection from the very beginning. Here it should be noted that it is important to work with nontrivial torsion configurations in order to find exact solutions in string gravity and gauge gravity models.
In this paper, we apply methods in the geometry of nonholonomic and almost-K ähler manifolds in order to study heterotic supergravity derived in the low-energy limit of heterotic string theory [45][46][47]. This publication is associated with another paper [48], where an approach to heterotic string gravity is formulated in the language of nonholonomic and almost-Kähler geometry. We cite also [1,Section 4.4] for a summary of previous results and similar conventions on warped configurations and modified gravitational equations. 3 The main goal of this work is to develop a geometric method for integrating in generic off-diagonal forms, and for generalized connections, the equations of motion of heterotic supergravity, up to and including terms of order α . As a secondary goal, we shall construct explicit examples of exact solutions describing nonholonomic deformations of the Kerr metric. In general, it is possible to formulate conditions for effective sources and generic off-diagonal when string gravity may encode/mimic equivalent solutions in massive gravity and/or modified f (R, T ) gravity. For reviews and original results related to massive and other types MGTs, we cite references [20][21][22][23]49,50].
In this work, a series of exact and/or small parameter depending solutions which for small deformations mimic rotoid Kerr-de Sitter-like black holes/ellipsoids, self-consistently embedded into generic off-diagonal backgrounds of 10-dimensional spacetimes are constructed. Such backgrounds can be of solitonic/vertex/instanton type. We study string gravity modifications with respect to nonholonomic frames and via re-definition of generating and integrations functions and coefficients of sources. These modifications can be analyzed in the framework of Einstein gravity but modeled by effective polarized cosmological constants and off-diagonal terms. Using solutions for heterotic string gravity, it is possible to mimic physically important effects in modified gravity. In a series of associated papers, we studied the acceleration of the universe, certain dark energy and dark matter locally anisotropic interactions, effective renormalization of quantum gravity models [8,24,51] via nonlinear generic off-diagonal interactions on effective Einstein spaces. These constructions can developed for models elaborated in the framework of string theory.
The solutions of heterotic supergravity which are constructed in this and the associated [48] work describe (1+3)dimensional walls endowed with generic, off-diagonal metrics, warped to an almost-Kähler 6-d internal space in the presence of nonholonomically deformed gravitational and gauge instantons. The generalized instanton contributions are adapted to a nontrivial nonlinear connection structure determined by generic off-diagonal interactions which allows us to solve the Yang-Mills, YM, sector and the corresponding Bianchi identity at order α (related to the gravitational constant in 10-d). Such 10-d solutions preserve two real supercharges, which correspond to the N = 1/2 supersymmetry. The almost-Kähler internal 6-d structure can be defined for various classes of solutions in 10-d gravity if we prescribe an effective Lagrange type generating function. In this approach we can work both with real nonholonomic gravitational and YM instanton configurations to consider deformed SU (3) structures.
The paper is organized as follows: We begin Sect. 2 with a summary on formulation of the heterotic supergravity theory in nonholonomic variables which was performed in [48]. This will allow us to derive a general decoupling property of motion equations in further sections. The geometric formalism on nonlinear and distinguished connections and adapted metrics to nonholonomic 2 + 2 + 2 + · · · splitting of higherdimension spacetimes is outlined. We develop the AFDM as a geometric method for constructing exact solutions of equations of motion in heterotic string theory and related 4d-10d modified Einstein gravitational equations. The contributions of gauge-like NS 3-forms, curvature of interior almost-Kä hler configurations, effective scalar and gauge fields etc. are encoded into certain effective sources and generic off-diagonal terms of metrics, with effective N-connection structure, defining a nontrivial vacuum structure. Such generalized/modified gravitational equations are formulated in adapted variables and with a very general off-diagonal ansatz for the metrics, (non-) linear connections and effective matter fields, which allows a decoupling of corresponding nonlinear systems of PDEs in very general forms. We show that using this nonholonomic decoupling property, it is possible to construct classes of exact solutions depending on various sets of generating and integration functions and integrations constants, on all 10-d spacetime coordinates. The existence of a very important nonlinear symmetry is proven. This allows a re-definition of the generating functions and effective sources to other equivalent data with effective cosmological constants induced by off-diagonal, or warped, and effective sources interactions. It is shown how the geometric constructions can be performed for the "simplest" case of one Killing symmetry in the 4-d case and then generalized to non-Killing configurations and higher dimensions. The conditions of generating solutions with zero torsion are analyzed. A self-consistent formalism of constructing small N-adapted nonholonomic stationary deformations is also explored.
Section 3 is devoted to a rigorous geometric study of nonholonomic generic off-diagonal deformations of exact solutions in heterotic string gravity containing the Kerr solution as a "primary" configuration. We show by using the AFDM it is possible to generate the Kerr solution as a particular nonholonomically constrained case and small parametric deformations with a well-defined physical interpretation. Then we construct solutions with general off-diagonal deformations of the Kerr metrics in 4d to 10d effective gravity with heterotic string modifications. We provide examples of (non-Einstein) metrics with nonholonomically induced torsions and study small off-diagonal modifications of the Kerr metrics determined by warped and general almost-Kähler internal space structures. Separate subsections are devoted to ellipsoidal 4d, 6-d, and 10-d deformations of the Kerr metric resulting in target vacuum rotoid, or Kerr-de Sitter, configurations, all self-consistently defining exact solutions of equations of motion in heterotic string theory.
Finally (in Sect. 4), we summarize the paper, provide conclusions and speculate on physical meanings of solutions with generic off-diagonal metrics and generalized connections constructed using the AFDM for the heterotic string theory.

Heterotic supergravity in nonholonomic variables
In this section, we outline a nonholonomic geometric approach to the heterotic supergravity modeled in the lowenergy limit of heterotic string theory as a N = 1 and 10d supergravity coupled to super Yang-Mills theory; see for details [48]. The nonholonomic variables will be parameterized in forms which allow a general decoupling of equations of motion and generating exact solutions depending, in principle, on all spacetime coordinates. Such a higher-dimension spacetime is modeled as a 10-d manifold M, equipped with a Lorentzian metricǧ of signature (+ + − + + + + + ++), with a time-like third coordinate. The heterotic supergravity theory is defined by a couple (M,ǧ), an NS 3-formȞ , a dilaton fieldφ and a gauge connection A∇ , with gauge group SO (32) or E 8 × E 8 . In our approach, we elaborate a system of notations which is adapted correspondingly for applications of geometric methods of constructing exact solutions in [40][41][42][43][44]. The notations in [48] and this work are different from the "standard" system of notations in string theory (see, for instance, [1]). A so-called nonholonomic "shell by shell", or 2 + 2 + · · · , splitting should be elaborated with corresponding left shell labels and shell indices. This minimizes the procedure of separating the equations of motion in certain general forms and for constructing exact solutions. In "shell" diadic variables, certain important (non) linear symmetries are explicitly shown and the type of generating and integration functions which can be considered are emphasized. Using only standard 4-d, 6-d and/or 10-d indices, like in former superstring and supergravity work, it is not possible to understand how the AFDM can be applied for generating off-diagonal solutions in (super-) string/gravity theories.
2.1 Geometric conventions on nonholonomic 2 + 2 + · · · splitting For geometric spacetime models on a 10-d pseudo-Riemannian spacetime M with a time-like coordinate u 3 = t and other coordinates being space-like, 4 we consider conventional splitting of dimensions dim M = 4 + 2s = 10; s = 0, 1, 2, 3. The AFDM, allows us to construct exact solutions with arbitrary signatures of metricsǧ, but our goal is to consider extra-dimensional space generalizations of the Einstein theory to heterotic supergravity models. In most general forms, this is possible if we use the formalism of nonlinear connection splitting for higher-dimensional (super-) spaces and strings, which was originally considered in (super-) Lagrange-Finsler theory [53,54]. We shall not consider Finsler type gravity in this work, but we follow an approach re-defined for nonholonomic distributions on (super-) manifolds [40,43,44]. The same geometric technique can be applied for tangent (super-) bundles or for (super-) manifolds enabled with certain classes of nonholonomic distributions like nonholonomic frames, nonlinear connections etc. For vector/tangent (super-) bundles, certain x-coordinates are used on the base (super-) manifold, but y-coordinates are considered for the typical fiber (super-) vector space. In the case of fibrations, the (x, y)-coordinates are used for the definition of a fibered structure. We consider (abstract, or coordinate) indices and coordinates u α s = (x i s , y a s ) for an oriented number of twodimensional, 2-d, "shells" added to a 4-d spacetime of signature (+ + −+). For s = 0, we write u α = (x i , y a ) and then extend "shell by shell" to a local system of 10-d coordinates, The corresponding subsets of indices are labeled in the form: 3 . . . = 9, 10; and, for instance, , or we shall write only i s . For brief denotations, we shall write 0 u = ( 0 x, 0 y); 1 u = ( 0 u, 1 y) = ( 0 x, 0 y, 1 y), 2 u = ( 1 u, 2 y) = ( 0 x, 0 y, 1 y, 2 y) and 3 u = ( 2 u, 3 y) = ( 0 x, 0 y, 1 y, 2 y, 3 y). Here we note in modern gravity, the so-called ADM (Arnowitt-Deser-Misner) formalism is largely used with 3 + 1 splitting, or any n + 1 splitting; see for details [52]. For the purposes of this work, such splittings are not convenient because it is not possible to elaborate a technique of general decoupling of the gravitational field equations and generating off-diagonal solutions. In these cases, the conventional one-dimensional "fibers" result in certain degenerate systems of equations. To construct exact solutions in four-ten-dimensional theories, it is more convenient to work with a correspondingly defined non-integrable 2 + 2 + · · · splitting; see details in [39,40,43]. For the heterotic supergravity, such geometric constructions are explored in greater detail in [48]. In order to connect "shell by shell" indices and coordinates to standard ones for supergravity theories (see [1,2]), we can consider small Greek indices without sub indices, and respective coordinates x μ , when α, μ, . . . = 0, 1, . . . , 9. The identification of shell coordinates with the standard ones follows such a rule: x 0 = u 3 = t (time-like coordinate) and (for space-like coordinates): On M, we can consider local frames/bases, e α s = e α s α s ( s u)∂/∂u α s , where the partial derivatives ∂ β s := ∂/∂u β s define local coordinate bases. We shall underline indices if it is necessary to emphasize that such values are defined with respect to a coordinate frame. The (co) frames, e α s = e α s α s ( s u)du α s , can be defined as dual to respective e α s . For our purposes, it is convenient to work with nonholonomic (non-integrable) distributions defining a 2 + 2 + · · · spacetime splitting. Such a distribution can be introduced to define a nonlinear connection, N-connection, structure via a Whitney sum This formula states a conventional horizontal (h) and vertical (v) "shell by shell" splitting. We shall use boldface letters in order to define spaces and geometric objects enabled/adapted to a N-connection structure. In local form, an N-connection is defined by its coefficients N a s i s when s N = N a s i s ( s u)dx i s ⊗ ∂/∂y a s .
A nonholonomic manifold is enabled with a nonholonomic distribution of type (2), when (for instance, for the "zero" shell) V hM ⊕ vM. The term N-anholonomic manifold is also used. This definition comes form the fact that the N-connection coefficients determine a system of Nadapted local bases, with N-elongated partial derivatives, e ν s = (e i s , e a s ), and cobases with N-adapted differentials, e μ s = (e i s , e a s ), For s = 0, , e a 0 = ∂ ∂ y a 0 , e i 0 = dx i , e a = e a 0 = dy a + N a i dx i ; or/and , e a s = ∂ ∂ y a s , The anholonomy relations are computed W b s i s a s = ∂ a s N b s i s and W a s j s i s = a s i s j s , where the curvature of N-connection is defined as the Neijenhuis tensor, a s i s j s := e j s (N a s i s ) − e i s (N a s j s ). Any metric structure on M can be parametrized via nonholonomic frame transforms as a distinguished metric (dmetric, in boldface form), s g = {g α s β s }, s g = g i s j s ( s u) e i s ⊗ e j s + g a s b s ( s u)e a s ⊗ e b s = g i j (x) e i ⊗ e j + g ab (u) e a ⊗ e b + g a 1 b 1 ( 1 u) e a 1 ⊗ e b 1 + · · · + g a s b s ( s u)e a s ⊗ e b s .

(5)
Introducing dual coframes e μ s = (e i s , e a s ) decomposed with respect to N-elongated differentials du a s as in (3), we obtain a generic off-diagonal form s g with respect to a dual coordinate frame. A metric s g can be diagonalized by coordinate transforms in a finite region if and only if the anholonomy coefficients in (4) vanish, i.e. W γ s α s β s = 0 A linear connection s D can be introduced on M in standard form. A distinguished connection, d-connection, s D is a linear connection preserving the N-connection spitting under parallelism. One defines the curvature R α s β s and torsion T α s in standard forms; see details in [48]. In this paper, we shall work with two very important linear connection structures determined by the same metric structure. These linear connections are uniquely defined following the geometric conditions: Here we note that the LC-connection s ∇ = { α s β s γ s } can be introduced without any N-connection structure. It can always be canonically distorted to a necessary type of d-connection It should be noted that this torsion is different form torsions in Einstein-Cartan gauge type and string gravity theories with absolutely antisymmetric torsion. Additional sources are not necessary because a d-torsion (7) is determined by the nonholonomic structures. In generalized theories, the torsion fields which are independent from the metric and vielbein fields may posses proper sources. Considering additional assumptions, we can relate the values, (7) for instance, to a subclass of nontrivial coefficients of an absolute antisymmetric torsion in string gravity. We can always extract LCconfigurations with zero torsion if we additionally impose for (7) the conditions Such additional nonholonomic constraints may be stated in non-explicit forms and without certain limits with small parameters and smooth functions. It should be noted that, in general, W γ s α s β s (4) may not be zero even if the conditions (8) are satisfied.
Any (pseudo-) Riemannian geometry can be equivalently formulated in nonholonomic variables ( s g, s N, s D) or using the standard data ( s g, s ∇). Because both linear connections s ∇ and s D are defined by the same metric structure, there is a canonical distortion relation The N-adapted coefficients of a curvature d-tensor R α s β s = {R α s β s γ s δ s } of the canonical d-connection s D and s g can be computed, respectively, for all shells s = 0, 1, 2, 3; see for details [48]. The Ricci d-tensor Ric = { R α s β s := R τ s α s β s τ s } is a respective contracting of the coefficients of curvature tensor, Considering the inverse d-metric to s g, we define and compute the scalar curvature of s D, with respective h-and v-components of the scalar curva-

The AFDM for heterotic supergravity
We develop the "anholonomic frame deformation method", AFDM, and apply these geometric techniques for applications to 10-d gravity and heterotic supergravity formulated in nonholonomic variables. The heterotic supergravity field equations were formulated in N-adapted form in [48]. They include terms of order α , equivalent to the equations of motion of heterotic nonholonomic supergravity considered in [1]. In explicit form, where the Hodge operator * , s D = { D μ s }, the canonical nonholonomic d'Alembert wave operator := g μ s ν s D μ s D ν s , the Ricci d-tensor R μ s ν s , and scalar s R are determined by a d-metric g (5). The curvature d-tensor R μ s α s β s γ s is taken for an almost-Kähler structure on shells s = 1, 2, 3 as we described above. The gauge field A corresponds to the Nadapted operator s A D and curvature F = F( 1 ψ) via a map constructed in [1,48]).
Equation (12) can be written as effective Einstein equations for the canonical d-connection s D, The sources ϒ β s δ s can be formally defined as in GR but for extra dimensions and in N-adapted form, when The system (16) can be derived and formulated in variational N-adapted form by considering a nonholonomic gravitational Lagrange density of type g L = R and an effective Lagrange density for matter m L. For simplicity, we shall consider Lagrangians depending only on the coefficients of metric field and matter field but not on their derivatives (for such configurations, it will be possible to construct exact solutions in explicit form). The energy-momentum d-tensor is computed by the definition, for | s g| = det | s g μν |. We conclude that by following an N-adapted variational calculus with action we can elaborate a 10-d nonholonomic gravity theory with gravitational field equations (16). A changing of geometric data ( s g, s D) → ( s g, s ∇) is possible via the canonical distorting relations (9), or imposing the zero torsion condition at the end T = 0 (8) for extracting LC-configurations s D | T =0 = s ∇ stated by (17). The matter Lagrange density m L can be chosen in such a form that, via corresponding frame transforms, the source ϒ β s δ s will encode the terms with the effective matter fields φ, H α s β s μ s , interior space fields R μ s α s β s γ s , and gauge fields F μ s α s . The nonholonomic equations of motion in heterotic string gravity (12) can be written in the form (16) with an effective source (related via nonlinear transforms of generating functions to certain effective cosmological constants), The traces of the above sources, and the respective effective cosmological constants, are related via the condition (13).

Decoupling and integration of nonholonomic equations of motion
We show how the heterotic string gravitational field equations (12) written in the form (16) with sources (19) and possible constraints (17), can be formally integrated in very general forms for generic off-diagonal metrics with coefficients depending on all spacetime coordinates. Additional conditions for extracting LC-configurations with s ∇ will be analyzed at the end after certain classes of general solutions have been constructed.

Ansatz for metrics, N-connections, and gravitational polarizations
In the simplest form, the decoupling property for any shell s = 0, 1, 2, 3 can be proven for certain ansatz with at least one Killing symmetry on the corresponding shell. Using an additional conformal transform on the necessary shells, we can extend the constructions for non-Killing configurations. Such a procedure is described for arbitrary finite shell s in Refs. [40][41][42][43]. In this work, we state the formulas for 10-d gravity and heterotic string theory in explicit form when the prime solutions are stationary ones, i.e., do not depend on the time-like coordinate y 3 = t. Let us consider metrics of type (5), which via frame transforms g α s β s = e α s α s e β s β s g α s β s can be parametrized 5 where e a = dy a + N a for N 9 α 2 = 3 n α 2 (u β 2 , y 10 ), N 10 Such an ansatz has a Killing vector ∂/∂y 9 because the coordinate y 9 is not contained in the coefficients of such metrics. We propose that via coordinate transforms we can eliminate dependence on y 3 = t and can parametrize h 4 = 1 and h a s = h(y 4 )ȟ a s (x i , y b s ) for shells s = 1, 2, 3, as it was considered in [48], if the configurations are with warping on y 4 . The ansatz (24) can be considered as a target d-metric of a prime d-metric with flat domain wall considered in that associated paper. In this work, we restrict our considerations to stationary solutions in heterotic string gravity (which do not depend on t). Inhomogeneous and locally anisotropic cosmological configurations in such string models, with generic dependence on t (see examples in [26][27][28][29][30]) will be studied in our future publications.

N-adapted sources and nonholonomically modified Einstein equations
We will be able to integrate nonholonomic equations of motion in heterotic string gravity in explicit form for very general assumptions if the source ϒ β s δ s (19) is parametrized in N-adapted form. This can be determined by five independent effective sources, respectively, on h-and s v-subspaces, which will be related to certain effective cosmological constants corresponding to the formulas For certain general configurations, it will be possible to fix generating and effective sources of type h = s = for all s, when a value for the corresponding contribution of fields can be zero, or nonzero, for instance, This depends on the type of vacuum, nonvacuum, or effective vacuum model we study. It is possible to compensate for contributions into an effective source of a field with contributions of another field, for instance, to get = φ + H + F + R = 0 even when not all the values of such effective cosmological constants are zero. We note that by prescribing certain values of sources (34) we can relate via nonholonomic frame transforms (in coordinate and/or N-adapted form) (19), in heterotic string gravity, or (18), in a 10-d nonholonomic generalization of GR. For any effective m L we can solve a system of quadratic algebraic equations for e α s α s in a form compatible with transforms to N-adapted frames of the metric/d-metric components, In this section, we shall construct general solutions of the generalized nonholonomic Einstein equations (16) with Ricci d-tensors (25)-(33) and effective sources (encoding contributions from heterotic supergravity) (19) parameterized in the form (34) written in N-adapted form as Similar equations can be written recurrently for arbitrary finite extra dimensions.

The ansatz for effective matter fields in heterotic string gravity
Assumptions on nontrivial φ-configurations We shall find dmetrics of type (5) for which the contributions of the field φ can be included by nonholonomic deformationsg → s g into certain generic off-diagonal terms, i.e. into N-connection coefficients. For instance if under a prime configuration This equation is equivalent to a linear system of equations, which can be solved in explicit form if the N-connection coefficients w i are defined (see next sections for how such values can be found in explicit form). For such configurations, we state that φ ϒ μ s ν s and φ in (19) and (20) are zero and the fields φ with such configurations contribute to possible heterotic supergravity effects only via possible nontrivial off-diagonal interactions and not via effective sources. The terms with φ also vanish in all nonholonomic equations of motion (12)-(15) written with respect to N-adapted frames if the conditions (36) are satisfied. Nontrivial coefficients with φ and its partial derivatives will appear if certain physical equations are written, for instance, in coordinate frames. We have chosen a special type of configurations for φ in order to simplify the procedure of finding generic off-diagonal solutions in heterotic string gravity in explicit form.

Nonholonomic gauge configurations In general form, we can consider any
Nonholonomic deformations of string gauge fields can be parametrized in the form where, respectively, For instance, kink-like almost-Kähler configurations can be encoded as a prime configuration 0 H α s β s μ s modifying for the target configuration H α s β s μ s the solutions for off-diagonal metrics via an effective source. Such parametrizations are possible by considering the ansatz with effective cosmological constants with absolute antisymmetric -tensors (we refer the reader to similar details for the 4-d case to [42] and the references therein). For such an ansatz, we obtain effective energymomentum tensors were we considered Eqs. (34) and (19)- (23) and, for instance the numbers n F = tr[internal F] and n R = tr[internal R] depends on the representation of the Lie algebra for F and on the representation of Lie groups on the internal space.

Decoupling of nonholonomic equations of motion and effective
Considering the ansatz (24) for g i (x k ) = i e q(x k ) , i = 1, in (35) with respective sources, we obtain this nonlinear system of PDEs: In equivalent form, the same equations are obtained recurrently if we write, for instance, s w i s instead of s w τ s−1 , s n i s instead of s n τ s−1 etc. (some equations and solutions can be parametrized in more simple forms if we follow the first convention, other equations will be more "compact" if we follow the second convention). In these equations, the generating functions (we shall use this in the next sections, on convenience formulas with -and/or -values) and the α-, β-, γ -coefficients on corresponding shells are defined, respectively: when the frame/coordinate systems are chosen in such a way that nonzero conditions for the partial derivatives are satisfied: = 0 (this allows us to avoid singular nonholonomic deformations).
Equations (39)-(47) prove a very important decoupling property of the heterotic string equations (12) and 4d-10d (generalized) Einstein equations (16) with respect to correspondingly N-adapted frames. In explicit and in certain simple forms, such formulas can be obtained for metrics with at least one Killing symmetry. Nevertheless, the constructions can be generalized for non-Killing configurations in any finite extra-dimension gravity; see a corresponding technique in [40,43]. In brief the decoupling property of the AFDM can be explained for the 4-d configurations: 1. Equation (39) is just a 2-d Laplacian, which can be solved for any h-source h ϒ(x k ). (40) contains only the partial derivative ∂ 4 , equivalently * -derivative, and constraints by a system of two equations, together with the formula for the value y 4 ) and (x i , y 4 ) and source ϒ(x k , y 4 ). Prescribing any two such functions, we can define (integrating on y 4 ) two other such functions. We note that h a are coefficients of a d-metric, is a generating function, and ϒ is related to extradimensional string contributions. 3. Using h 3 and , or , from the previous point, we compute the coefficients α i and β; see (48). This allows us to define w i from the algebraic equations (42). 4. Having computed the coefficient γ (48), the N-connection coefficients n i can be defined after two integrations on y 4 in (41).

Equation
We can repeat the steps 2-4 recurrently on shells s = 1, 2, 3, adding dependencies on 2 + 2 + 2 extra-dimension coordinates, respectively. In this way, we involve new classes of generating functions and effective sources. Solving the systems (43)-(47) recurrently, we can compute all d-metric and N-connection coefficients for a 10-d ansatz (5).

Integration of nonholonomic equations of motion by generating functions and effective sources
The system of nonlinear PDEs (39)-(47) with coefficients (48) and effective sources (37) with contributions from string gravity can be integrated in general forms on any shell up to 10-d.

4-d non-vacuum configurations
The coefficients g i = e q(x k ) are defined by solutions of the corresponding 2-d Poisson equation (39) as we mentioned above (see point 1 at the end of previous subsection). We can integrate the system of nonlinear PDEs consisting of the first equation in (48) and (40) for arbitrary source ϒ(x k , y 4 ) and with generating function (x k , y 4 ) = e . The solutions will be generated for a stationary d-metric when the coefficients do not depend on time-like coordinate y 3 = t, when * = 0. The vertical effective gravitational field equations (40)- (42) This system can be solved in very general forms by prescribing ϒ and and integrating the equations "step by step". Introducing the function (the sign − is motivated by the pseudo-Euclidean signature), we express (49) and (50) Using h * 3 from the second equation (54) in the first equation, we write This value, together with the second equation of (54) and a further integration on y 4 , result in where is an integration function. Considering Eqs. (55), (53), and (56), we compute The first parts of the N-connection coefficients are found by integrating two times on y 4 in (51) written in the form for the coefficient γ defined in (48). Integrating two times on y 4 for explicit values of (57) and (56), we compute containing also a second set of integration functions 1 n k (x i ) and the redefined second integration function 2 n k (x i ).
We can solve the linear algebraic equations (52) and express Putting together all the above formulas and writing the effective source (34) in explicit form, we obtain the formulas for the coefficients of a d-metric and a N-connection determining a class of stationary solutions for the system (39)- (42), as a solution of 2-d Poisson equations (39); Using coefficients (58), we define such a class of quadratic elements for off-diagonal 4-d stationary configurations in heterotic supergravity with nonholonomically induced torsion, Such a class of metrics also contains exact solutions for the canonical d-connection D in R 2 gravity with effective scalar field encoded into a nonholonomically polarized vacuum with a special parametrization of the source ϒ; see for details [14,70].

Examples of 4-d-10-vacuum configurations
The configurations with zero source cannot be constructed as particular cases of former off-diagonal solutions with ϒ = 0 because these limits are not smooth for metrics (59). In string heterotic gravity, such conditions can be satisfied if different effective fields compensate their mutual contributions and result in an effective vacuum gravitational configuration. For the ansatz (24), we can construct solutions when the nontrivial coefficients of the Ricci d-tensor (25) In 4-d, there are three general classes of off-diagonal metrics which result in zero coefficients (26)-(28). Such constructions can be generalized for 10-d generic off-diagonals with mixing standard vacuum configurations for the 4-d spacetimes.
1. If h * 3 = 0 but h 3 = 0, h * 4 = 0 and h 4 = 0, we obtain only one nontrivial equation (27), where h 4 (x i , y 4 ) and w k (x i , y 4 ) are arbitrary generating functions. Integrating two times on y 4 , we obtain n k = 1 n k + 2 n k dy 4 / h 4 (60) with integration functions 1 n k (x i ) and 2 n k (x i ). In 4-d, this defines a quadratic line element Recurrently  , y 4 , y a 1 , y a 2 , y 10 ) dy 10 For instance, the bracket {...} type1,s=2 states that the vacuum metric is of type 1 on shell 2. The coefficients of such d-metrics do not depend on variables (t, y 9 ) and, on respective shells, the values 0 h 3 , 1 n k , 2 n k , 0 h 5 , 1 1 n k 1 , 1 2 n k 1 , 0 h 7 , 2 1 n k 2 , 2 2 n k 2 , 0 h 9 , 3 1 n k 3 , 3 2 n k 3 are integration functions and q, h 4 is any generating function. The coefficients n k should be found from (27); see (60). Such a family of 4-d vacuum generic off-diagonal metrics is described by where the integration function 2 n k (x i ) contains certain integration coefficients. We can extend such metrics on any shell s = 1, 2, 3 preserving the conditions of a zero effective source and adding respective generating functions gener. functs.
In 10-d, the vacuum solutions (62) are generalized to for integration functions In the above shell formulas, h * 1 5 = 0 and h * 1 6 = 0; h * 2 7 = 0 and h * 2 8 = 0; h * 3 9 = 0 and h * 3 10 = 0. 3. We consider that h * The coefficients n i are found from (27) and/or, equivalently, to (41) with γ = 3 2 |h 3 | * . Integrating on y 4 , this subclass of N-coefficients are computed with integration functions 1 n i (x k ) and 2 n i (x k ), or redefined 2 n i = − 2 n i /2c 2 . The quadratic line element for this class of solutions for vacuum metrics is described by On extra shells, this type of nonholonomic vacuum solutions are given by quadratic elements The generating functions are and the integration functions are , y a 1 ), 2 c 2 (x k , y 4 , y a 1 ), 2 1 n i 2 (x k , y 4 , y a 1 ), 2 2 n i 2 (x k , y 4 , y a 1 ), 0 h 8 (x k , y 4 , y a 1 ); 3 c 1 (x k , y 4 , y a 1 , y a 2 ), 3 c 2 (x k , y 4 , y a 1 , y a 2 ), 3 1 n i 2 (x k , y 4 , y a 1 , y a 2 ), 3 2 n i 2 (x k , y 4 , y a 1 , y a 2 ), 0 h 10 (x k , y 4 , y a 1 , y a 2 ). In certain cases, the type may be the same on two or three shells, but the parametrizations of mixed types of nonholonomic vacuum solutions is given by different dependencies on shell coordinates of the generating and integration functions.

Extra-dimensional non-vacuum solutions
The solutions for extra dimensions can be constructed in certain forms which are similar to the 4-d ones with (58) using new classes of generating and integration functions with dependencies on extra dimensional coordinates. For instance, we can generate solutions of the system (43)- (45) with coefficients (48) on the shell s = 1 following a formal analogy when ∂ 4 → ∂ 6 , i.e. * → * 1 , (x k , y 4 ) → 1 (x k 1 , y 6 ), with 1 (x k 1 , y 4 ) = e 1 ; ϒ(x k , y 4 ) → 1 ϒ(x k 1 , y 6 ) · · · and applying the same procedure as for the 4-d case but extending the number of shell coordinates, respectively, for generating and integration functions.
The solutions on s = 1 are also generated as stationary d-metrics when the coefficients do not depend on timelike coordinate y 3 = t and do not depend on the coordinate y 5 . For nontrivial 1 ϒ, we have to choose 1 * 1 = 0. The effective gravitational field equations (43)-(45) and respective equation for 1 from (48) transform (respectively) into Prescribing 1 ϒ and 1 , such equations for 6-d gravity can be integrated as in the case of s = 0 with two vertical coordinates in 4-d gravity, i.e. integrating equations "step by step". Introducing the function ( 1 ρ) 2 := h 5 h 6 , (for extra shells the sign is different because we work with space-like coordinates) we express the first two equations above as Using h * 1 5 from the second equation in the first equations, we get 1 ρ = 1 2 1 * 1 1 ϒ . Substituting this value together with the second equation of (64) and integrating on y 6 , it is possible to compute is an integration function. In result, we find from the above formula ( 1 ρ) 2 the coefficient The first part of the N-connection coefficients are found by integrating two times on y 6 in (44) written in the form for the coefficient 1 γ defined in (48). This formula can be integrated two times on y 6 for explicit values of h a 1 which results in 1 n k 1 (x i , y 4 , y 6 ) = 1 1 n k 1 + 1 2 n k 1 dy 6 h 6 | h 5 | 3/2 = 1 1 n k 1 + 1 2 n k 1 dy 6 , also containing a second set of integration functions 1 1 n k 1 (x i , y 4 ) and redefined second integration function 1 2 n k 1 (x i , y 4 ). The second set of N-connection coefficients on s = 1 can be found from the linear algebraic equations (45) and express 1 w i 1 = ∂ i 1 1 / 1 * 1 . Summarizing the above formulas for the d-connection and N-connection coefficients on s = 1, we write the data for generating such generic off-diagonal metrics as stationary solutions for the system (43)- (45), The quadratic elements for off-diagonal 6-d stationary configurations in heterotic supergravity with nonholonomically induced torsion are constructed using above N-adapted coefficients We can repeat the method of constructing solutions for s = 1 to next shells s = 2 and s = 3 and integrate the systems (46) and (47), respectively. It is necessary to extend the formulas for generating and integration functions and effective sources recurrently as we provided above for redefinition of respective values for s = 0 to s = 1. The quadratic elements are parametrized in corresponding forms: and The generic off-diagonal d-metrics (67) define exact stationary solutions with Killing symmetries on ∂ t and ∂ 9 in effective nonholonomic 10-d gravity with sources (34) being determined by the equations of motion in heterotic supergravity (12).

A nonlinear symmetry of generating functions and effective sources
There is an important nonlinear symmetry relating nontrivial generating functions and effective sources considered in the above classes of solutions, used in differen forms in [42,43]. It allows us to redefine the generating functions and introduce effective cosmological constants instead of effective sources. Let us study these properties in the case of heterotic string gravity.
Changing the generating data ( , ϒ) ↔ (˜ , = const), where , which is equivalent to 2 = dy 4 ϒ −1 ( 2 ) * and/or up to certain classes of integration functions depending on coordinates x i , we express the solutions (56) and (57) of the system of nonlinear PDEs (49) and (51) in two equivalent forms, Having defined two equivalent nonlinear formulas for h a [ ] = h a [ ], we can express in two equivalent forms the Nadapted coefficients for the d-metric and N-connection (58).
In the second case, the N-connection coefficients are computed using integrals on dy 4 for certain values determined by (˜ , ) via (68) and integration functions, We observe that if the effective source ϒ = ϒ(x k ) does not depend on y 4 , we have the same expression for w i in terms of generating functions, w i = ∂ i * = ∂ i * . The nonlinear symmetry reflects the property of changing the effective sources mentioned in (34), We can identify with , or any other value H , F , R , and their sums, depending on the class of models with effective gauge interactions we consider in our work. In a similar form, using recurrent formulas, we can prove the existence of such nonlinear symmetries generalizing (68), We consider the convention that ϒ 9 9 = ϒ 10 10 = 3 ϒ(x k , y a ,  y a 1 , y a 2 , y 10 and identify 3 with 3 . Such nonlinear transforms can be used for simplifications of formulas for generic off-diagonal solutions. Prescribing certain effective matter field initial distributions, we can redefine the constructions equivalently, via new classes of generating functions, for effective cosmological constants. Finally, we note that Eq. (67) simplifies for Such generic off-diagonal configurations are described by generic off-diagonal metrics with effective cosmological constants, generalized generating functions and integration functions. The contributions of 3 ϒ(x i , y 4 , y a 1 , y a 2 ) can be encoded into generating and integrating functions.

The Levi-Civita conditions
In general, a solution constructed for a generic off-diagonal metric (24) and canonical d-connections s D is characterized by nonholonomically induced d-torsion coefficients T γ s α s β s (7) completely defined by the N-connection and dmetric structure. The zero torsion conditions (8) can be satisfied by a subclass of nonholonomic distributions determined by corresponding parametrizations of the generating and integration functions and sources. By straightforward computations (see the details in Refs. [35,[38][39][40]), we can verify that if the coefficients of N-adapted frames and s vcomponents of d-metrics are subject to the respective conditions all d-torsion coefficients are zero. The n-coefficients solve the conditions (71) if The explicit form of solutions of constraints on w k derived from (71) depends on the class of vacuum or non-vacuum metrics we try to construct; see for details [43]. For instance, if we choose a generating function =ˇ (x i , y 4 ), for which (∂ iˇ ) * = ∂ iˇ * , we solve the conditions for w i in (71) Imposing respective conditions from (72) and (73) on the coefficients of (58), we define such a class of quadratic ele-ments for off-diagonal 4-d stationary configurations in heterotic supergravity with zero induced torsion, We can impose similar conditions and generate exact offdiagonal solutions with respective data ( sǧ , sŇ , s∇ ) for which all nonholonomic torsions are zero, but this will impose very strong restrictions on the dynamics of effective matter fields on extra shells in heterotic supergravity. In order to consider realistic solutions in string gravity with 6-d interior almost-K ähler configurations, the torsion is positively not zero. Such 10-d solutions are conventionally parameterized as (s = 0 :ǧ, N , ∇; s = 1, 2, 3 : s g, s N, s D), for stationary configurations with Killing symmetry on ∂ t and ∂ 9 . Hereafter we shall work with solutions with nontrivial nonholonomically induced torsions considering that it is always possible to state additional constraints resulting in LC-configurations in the 4-d case or extra dimensions.

Small N-adapted nonholonomic stationary deformations
We can construct very general classes of generic off-diagonal stationary solutions in heterotic supergravity. It is not clear what physical meaning these configurations may have for generating and integration functions with arbitrary data. Using small off-diagonal deformations of some known physically important solutions we can understand physical properties of such solutions characterized by locally anisotropic polarization/running of constants and nonlinear off-diagonal gravitational interactions determined by (super-) string corrections. We consider a "prime" pseudo-Riemannian metric of type +g a 1 (x k , y 4 , y 6 ) e a 1 2 +g a 2 (x k , y 4 , y a 1 , y 8 ) e a 2 2 +g a 3 (x k , y 4 , y a 1 , y a 2 , y 10 ) wheree a s are N-elongated differentials. Such a metric is diagonalizable if there is a coordinate transform u α s = u α s (u α s ) when ds 2 =g i (x k )(dx i ) 2 +h a s (x k )(dy a s ) 2 , with sẘ i s = sn i s = 0. To construct exact solutions with non-singular coordinates it is important to work with "for-mal" off-diagonal parameterizations when the coefficients sẘ i s and/or sn i s are not zero but the anholonomy coeffi-cientsW α s β s γ s (u μ s ) = 0; see (4). We suppose that some data (g i ,h a s ) may define a known physically important diagonal exact solution in GR or heterotic string gravity (for instance, a black hole, BH, configuration of Kerr or Schwarzschild type). Our goal is to study certain small generic off-diagonal parametric deformations of the prime d-metric and N-connection coefficients (75) into certain target metrics where the coefficients (g α s = η α sg α s , w η i sẘ i s , n η i s n i s ) define, for instance, a d-metric s g (24) as a solution of nonholonomic equations of motion in heterotic string gravity (12) rewritten as a nonlinear system of PDEs in nonholonomic 10-d gravity (39)- (47). Let us construct exact solutions, for instance, of type (67). For certain well-defined conditions, we can show an expression using d-metric and N-connection deformations stated in explicit form on all shells for the d-metric and n-and wcoefficients: For the coefficients of d-metrics, (77) εχ a (x k , y 4 ), x k , y 4 , y 6 ), η a 2 = 1 + εχ a 2 (x k , y 4 , y a 1 , y 8 ), η a 3 = 1 + εχ a 3 (x k , y 4 , y a 1 , y a 2 , y 10 ); and for the coefficients of N-connection, , y 4 , y a 1 , y a 2 , y 10 ), for a small parameter 0 ≤ ε 1, when (76) transforms into (75) for ε → 0 (which in turn, can be diagonalized). In general, there are no smooth limits from such nonholonomic deformations which can be satisfied for arbitrary generation and integration functions, integration constants, and general (effective) sources on the corresponding shells. The goal of this subsection is to analyze such conditions when εdeformations with nontrivial N-connection coefficients for the prime and target d-metrics can be related to new classes of solutions of motion heterotic string equations.
We denote nonholonomic ε-deformations of certain prime d-metric (75) into a target one (76) with polarizations (77) in the formg → ε g = ( ε g i , ε h a s , ε N a s−1 b s ). The goal of this subsection is to compute the formulas for ε-deformations of prime d-metrics resulting in stationary solutions in equivalent nonholonomic 10-d gravity with Killing symmetries on ∂ t and ∂ 9 .
The geometric constructions will be provided in detail for the 4-d configurations and then extended to higherdimensional shells. Deformations of h-components are characterized by ε g i =g i (1 + εχ i ) = e q(x k ) resulting in a solution of the 2-d Laplace equation (39).
At the next step, we compute ε-deformations of vcomponents on the s = 0 shell, Parametrizing the generation function, and introducing this value in (79), one obtains In result, we can compute χ 3 for any deformation χ from a 2hypersurface y 4 = y 4 (x k ). Such a hypersurface, in general, is defined in non-explicit form from˚ =˚ (x k , y 4 ) when the integration function h [0] 3 (x k ), the prime valueg 3 (x k ), and the fraction (˚ 2 ) * /ϒ satisfy the condition (81).
We can find the formula for hypersurface˚ (x k , y 4 ) prescribing a value of ϒ. Introducing (64) into (79 ), one obtains i.e. we compute χ 4 for any data (˚ ,g 3 , χ). The formula for a compatible source is h ϒ = ±˚ * /2 |g 4 h [0] 3 |. It transforms (81) into a 2-d hypersurface formula y 4 = y 4 (x k ) defined in non-explicit form from The ε-deformations of N-connection coefficients w i = ∂ i / * for nontrivialẘ i = ∂ i˚ /˚ * are found following Eqs. (64) and (77), where there is no summation on index i. We can compute the deformations of the n-coefficients (we omit such details; see below for the necessary formulas).
In a similar way, the ε-deformations of 1 v-components on the s = 1 shell are computed, The first shell generation function is parameterized 1 → , resulting in It is possible to compute χ 5 for any deformation 1 χ from a 3-hypersurface y 6 = y 6 (x k , y 4 ) (for stationary solutions with Killing symmetries on ∂ t and ∂ 5 ). We note that, in general, such a hypersurface is defined in non-explicit form from 1˚ = 1˚ (x k , y 4 , y 6 ) when the integration function h [0] 5 (x k , y 4 ), the prime valueg 5 (x k , y 4 ) and the fraction ( 1˚ 2 ) * 1 / 1 ϒ satisfy a condition similar to (81).
We can find the formula for hypersurface 1˚ (x k , y 4 , y 6 ) by prescribing a first shell value of effective source 1 ϒ. Similarly to (64) and (79), it is possible to generalize and compute i.e. we compute χ 6 for any data ( 1˚ ,g 5 , 1 χ). One defines a first shell compatible source if 1 ϒ = ± 1˚ * 1 /2 |g 6 h [0] 5 |. It generalizes (81) into a 2-d hypersurface formula y 6 = y 6 (x k , y 5 ) which has to be computed in non-explicit form from dy 6 1˚ = ±(h [0] On the first shell, the ε-deformations of N-connection coefficients , where there is no summation on index i 1 . We omit computations of deformations of the n-coefficients but we shall present necessary formulas below, see similar details in [43]. Summarizing the above formulas and extending for all shells s = 0, 1, 2, 3, we obtain such coefficients for εdeformations of a prime metric (75) into a target stationary metric satisfying the equations of motion in heterotic supergravity: where n i (x k ) is a redefined integration function including contributions from the prime metric. On a shell s, these formulas are defined recurrently (we omit parameterizations of functions on shell coordinates because they can be determined in compatible form with indices and labels for shells), These formulas with s = 1, 2, 3 allow us to parametrize all coefficients of vertical components of d-metrics and Nconnections. For small parametric deformations, the values χ, s χ and ϒ, s ϒ have to be considered as generating functions. The values with a "circle" are prescribed by a chosen prime solution (in our case, we can chose the 4-d Kerr metric embedded into a 10-d gravity string spacetime). Fixing a small value ε, we compute such deformations and prove their stability (see [38] and the references therein) for any stable prime solution.

Nonholonomic heterotic string deformations of the Kerr metric
In this section, we study generic off-diagonal deformations and generalizations of the 4-d Kerr metric to new classes of exact solutions of equations of motion in heterotic string theory; see [31,32,52,71]. We prove that using the AFDM extended to models with almost-Kähler internal spaces, the Kerr solution can be constructed as a particular case by prescribing a corresponding class of generating and integration functions. In general, such solutions are with nontrivial cosmological constants and nonholonomically induced torsions. Imposing additional nonholonomic constraints, we can generate effective vacuum solutions and extract Levi-Civita configurations. A series of new classes of small parametric solu-tions when the Kerr metrics are nonholonomically deformed into general or ellipsoidal stationary configurations in fourdimensional gravity and/or extra dimensions are considered. We provide and study examples of generic off-diagonal metrics encoding nonlinear interactions with 3-form and gaugelike fields, nonholonomically induced torsion effects and instanton configurations.

Preliminaries on the Kerr vacuum solution and nonholonomic variables
A 4-d ansatz written in terms of three functions (h, Y, A) on coordinates x i = (ρ, z) defines the Kerr solution of the vacuum Einstein equations (for rotating black holes) if we choose Some values M = const and ρ = 0 result in a horizon x 1 = 0 and the "north/south" segments of the rotation axis, x 2 = +1/ − 1. For further applications of the AFDM, we can write this prime solution in the form This is possible if the coordinates x 1 ( x 1 , x 2 ) and x 2 ( x 1 , x 2 ) are defined for any x 2 , t). We can consider an N-adapted basis e 3 = dt + (∂ i y 3 )dx i , e 4 = dy 4 + (∂ i y 4 )dx i , for some functions y a , a = 3, 4, with ∂ t y 4 = −A(x k ). The Kerr metric was intensively studied in the so-called Boyer-Lindquist coordinates (r, ϑ, ϕ, t), for r = m 0 (1 + p x 1 ), x 2 = cos ϑ, which can be considered for applications of the AFDM. Such coordinates are expressed via parameters p, q which are related to the total black hole mass, m 0 and the total angular momentum, am 0 , for the asymptotically flat, stationary, and antisymmetric Kerr spacetime. The formulas m 0 = Mp −1 and a = Mqp −1 when p 2 + q 2 = 1 imply m 2 0 − a 2 = M 2 . In these variables, the metric (85) can be written for any coordinate functions x 1 (r, ϑ), x 2 (r, ϑ), y 3 = t + y 3 (r, ϑ, ϕ)+ϕ B/A, y 4 = ϕ, ∂ ϕ y 3 = −B/A, for which (dx 1 ) 2 +(dx 2 ) 2 = −1 dr 2 + dϑ 2 , and the coefficients are A = − −1 ( − a 2 sin 2 ϑ), We refer the reader to [31,52,71] for main results and methods for stationary black hole solutions (the coordinates x 1 , x 2 introduced above correspond to x, y, respectively, in chapter 4 of the first book). The primed quadratic linear elements (85) (or (86)) define solutions of vacuum Einstein equations parametrized in the form (16) and (17) with zero sources. A straightforward application of the AFDM is possible if we consider a correspondingly N-adapted system of coordinates instead of the "standard" prolate spherical, or Boyer-Lindquist system. Parametrizations (88) are most convenient for a straightforward application of the AFDM. This way we can generalize the solutions for coefficients depending on more than two coordinates, with non-Killing configurations and/or extra dimensions.
Working with general classes of stationary solutions generated by the AFDM, the Kerr vacuum solution in GR can be considered as a "degenerate" case of 4-d off-diagonal vacuum solutions determined by primary metrics with data (88) when the diagonal coefficients depend only on two "horizontal" N-adapted coordinates. Such a metric contains off-diagonal terms induced by rotation frames in a form when the nonholonomically induced torsion is zero. In N-adapted frames, further generic off-diagonal and extra dimension generalizations can be performed following standard geometric methods (see the following sections).

Off-diagonal deformations of 4-d Kerr metrics by heterotic string sources
Let us consider the coefficients (88) for the Kerr metric as the prime metricg (in general, a prime metric may or may not be an exact solution of the Einstein or other modified gravitational equations). Our goal is to construct nonholonomic deformations, see the sources (34), (37), and (19), respectively, with (21), (69). The main condition is that the target metric g is of type (74), which positively defines a torsionless off-diagonal solution of field equations in the 4-d gravity sector with sources determined from the heterotic string theory (12) with ansatz for the sources (38). The N-adapted deformations of coefficients of metrics, frames, and sources are parametrized in the form The primesg i ,h a ,ẘ i ,n i (88) are given by coefficients depending only on (x k ). The general deformations of the Kerr solution determined by generating function and extradimensional string sources with ansatz resulting in cosmological constants described in terms of polarization functions, where |η 4 |, may not have a smooth limit to 1, where the coefficient 1/4 was introduced inˇ and the func-tionǍ has any value for which ∂ i ˇ /ˇ * = ∂ i Ǎ .
Summarizing the above formulas, we obtain the quadratic element where , n k =n k + η n k = ∂ k (− y 3 + ϕ B/A + η n) and q is a solution of We used the important relationh 3 h 4 = AC − B 2 and emphasize that it is possible to take any function η n(x k ).
The solutions (90) are for stationary LC-configurations determined by off-diagonal heterotic string gravity effects on Kerr black holes when the new class of spacetimes come with Killing symmetry on ∂/∂y 3 and generic dependence on three (from maximally four) coordinates, (x i (r, ϑ), ϕ). Similar solutions were constructed and studied in massive gravity with extra dimensions [43]. Off-diagonal modifications are possible for any nontrivial values ofˇ and any small constants H , F , R . The solutions depend on the type of generating functionˇ (x i , ϕ) we fix in order to suit certain experimental/observational data for fixed systems of reference/coordinates. These can be re-parameterized for an effective K , which should also be compatible with experimental data. In such variables, we can mimic stationary heterotic string gravity effects by off-diagonal configurations in GR with integration parameters which should also be fixed by additional assumptions on symmetries of interactions. See Sect. 3.3 for ellipsoid configurations and details on parametric Killing symmetries in Refs. [33][34][35].

Nonholonomically string induced torsion for Kerr metrics in the 4-d sector
If we do not impose the LC-conditions (17), a nontrivial source ϒ → from heterotic string gravity induces stationary configurations with nontrivial d-torsion (7). The torsion coefficients are determined by metrics of type (59) with h ϒ = ϒ = K as in (89) and parametrizations of coefficients and coordinates distinguishing the prime data for a Kerr metric (88). Such solutions can be written in the form where the generating function (x i , ϕ) is not subject to any integrability conditions. Nontrivial stationary off-diagonal torsion effects may result in additional effective rotations if the integration function 2 n k = 0. Considering two different classes of off-diagonal solutions (91) and (90), we can study if heterotic string corrections to GR can have a nonholonomically induced torsion or if such effects are characterized by additional nonholonomic constraints as in GR (for zero torsion).
It should be noted that configurations of type (91) can be constructed in various theories with noncommutative and commutative variables. We can consider warped and trapped brane type variables in string, Finsler-like and/or Hořava-Lifshitz theories [24,38,51,55,56] when nonholonomically induced torsion effects play a substantial role.

Small modifications of Kerr metrics and effective string sources
It is not clear what physical meaning general deformations of the Kerr metric described by metrics of type (91) may have. We can choose certain subclasses of nonholonomic distributions describing stationary ε-deformations described by Eqs. (83). Using the Kerr solution (88) as a primary metric with assumptions for the string sources, we compute small deformations into d-metric and N-connection coefficients, The ε-deformations are computed thus: as a solution of 2-d Poisson equations (39); are kept as above but a coordinate transform is performed in order to haveN 4 i =ẘ i = 0. We use polarization functions linearized on ε, and for the coefficients of N-connection, Summarizing the above for the 4-d ε-configurations, we obtain the quadratic element In general, these heterotic string deformations of the Kerr metric are nonholonomically induced torsion coefficients, linear in ε.
We can impose additional constraints on χ, n i and˚ which allows us to extract LC-configurations. The corresponding ε-deformed analog of the metric (90) can be written where χ defines a deformed generating function ε = (x k , y 4 )[1 + εχ(x k , y 4 )] as in Eq. (80) but subjected to the condition which together with χ n(x i ) are chosen to result in zero nonholonomically induced torsion.

String induced ellipsoidal 4-d deformations of the Kerr metric
We provide some examples of how the Kerr primary data (88) are nonholonomically deformed by heterotic string interactions into target generic off-diagonal solutions of vacuum and non-vacuum Einstein equations for the canonical dconnection and/or the Levi-Civita connection. Generic offdiagonal metrics of type (92) can be parameterized as ellipsoidal deformations of the Kerr metric on a small eccentricity parameter ε.
As regards the result, the condition h 3 = 0, i.e. ϕ (r, ϕ, ε) = a 2 sin 2 ϑ, states an ellipsoidal "deformed horizon" If a = 0, we obtain the parametric formula for an ellipse with eccentricity ε, r + = 2m 0 1+εζ sin(ω 0 ϕ+ϕ 0 ) . Such configurations correspond to the generating function determined by effective heterotic string source K , as follows from (93). If the anholonomy coefficients (4) computed for (92) are not trivial for w i and n k = 1 n k , the generated solutions cannot be diagonalized via coordinate transforms.
The corresponding nonholonomically deformed 4-d spacetimes have one Killing symmetry on ∂/∂y 3 . For small ε, the singularity at = 0 is "hidden" under ellipsoidal deformed horizons if m 0 ≥ a. Similarly to the Kerr solution, there are ϕ-deformed cases for both the event horizon, r + = m(ϕ) + (m 2 (ϕ) − a 2 sin 2 ϑ) 1/2 , and the Cauchy horizon, which are effectively embedded into an off-diagonal background determined by N-coefficients. Using an ellipsoid type generating function (95) in (92), we construct a class of generic off-diagonal solutions of effective Einstein equations with heterotic string gravity effective cosmological constant K , which in its turn can be related to arbitrary sources via a redefinition of the generating functions (see Eqs. (70) and (37)) adapted to ε-deformations. The corresponding quadratic line elements are The new classes of ε-deformed solutions determine Kerrlike black hole solutions with additional dependencies on the variable ϕ of certain diagonal and off-diagonal coefficients of the metric. There is an obvious anisotropy in the angle ϕ. The values ζ and ω 0 have to be chosen in accordance with experimental data. The function˚ depends on corresponding frame distributions for the prime metric. Fixing a = 0 for a ε = 0, we get ellipsoidal deformations of the Schwarzschild black holes. We studied these constructions in detail in [38]; see also the references therein on the stability and interpretation of such solutions with both commutative and/or noncommutative deformation parameters. In general, a black hole/ellipsoid interpretation is not possible for "non-small" N -deformations of the Kerr metric. For certain embeddings, we can generate black hole-like configurations with deformed horizons and locally anisotropic polarized physical constants.

Ellipsoid Kerr-de Sitter configurations in R 2 and heterotic string gravity
The first examples of generic off-diagonal ellipsoid-solitonic deformations of similar Kerr-Sen black holes were constructed in [72]. Recently, asymptotically de Sitter solutions with spherical symmetry for R 2 gravity were studied in [14] and the nonholonomic geometric off-diagonal evolution of such metrics was analyzed in [70]. In this section, we show that those constructions can be related to the 4-d part of heterotic string MGTs. We consider a prime 4-d metric which for e √ 1/3φ = 3λ 2ς 2 = 1 8ς 2 R and g μν = e √ 1/3φ g μν = R 8ς 2 g μν , R = 0, defines exact solutions with spherical symmetry in R 2 gravity, for equations R μν = 2ς 2 g μν . The effective cosmological constant ς 2 is usually related to nonlinear scalar fields/dilaton-like interactions in effective Einstein gravity resulting from R 2 gravity. In our model of heterotic string gravity, we can choose and study quadratic gravity 4-d models determined by heterotic string effective sources. The metric (97) describes asymptotically de Sitter solutions with λ > 0 and R = 0.
Introducing new 4-d coordinates, the metric (97) is written as a "prime" metric for some constants M, λ and u α = ( x k , y a ). In order to work with a "formal" off-diagonal metric of type (75) with nontrivial valuesh * a ,ẘ i andn i , butW α βγ ( u μ ) = 0, see (4), we consider a coordinate transform u α = u α (u α ) with ϕ = ϕ(y 4 , x k ) and t = t (y 3 for ∂ i ϕ = ∂ϕ/∂ x i and ∂ a ϕ = ∂ϕ/∂y a . Choosing we express (97) as The prime d-metric (100) allows us to apply the AFDM and construct ε-deformations of geometric/physical objects and physical parameters as shown in Sect. 2.4.
The parametric equation of an ellipse with radial parameter r + = M and eccentricity ε, r + M 1−ε cos(ω 0 ϕ+ϕ 0 ) , can be determined in a simple way for λ = 0. We have to find solutions of a third order algebraic equation in order to determine possible horizons for nontrivial λ.
We construct ellipsoidal deformations of d-metric (100) if χ = ς χ = 8 M r ς 2˚ −2 cos(ω 0 ϕ + ϕ 0 ), with identification (98). Following the same method as in the previous subsection but for ς χ used for d-metric coefficients (83), we compute where n i (x k ) is a redefined integration function including contributions from the prime metric (100). The generating functions ς χ and 0 q can be determined for an ellipsoid configuration induced by the effective cosmological constant ς 2 in R 2 gravity. Finally, we generate a class of generic off-diagonal metrics for ellipsoid Kerr-de Sitter configurations related to the cosmological constant in heterotic string gravity, Such metrics have Killing symmetry on ∂/∂y 3 and are completely defined by generating functions 1 q and ς χ and effective source 2ς 2 = K induced from string theory. They define ε-deformations of Kerr-de Sitter black holes into ellipsoid configurations with effective (polarized) cosmological constants determined by constants in string theory and equivalent MGTs. If the LC-conditions are satisfied, such metrics can be modeled in GR with nontrivial cosmological constant.

Extra-dimensional off-diagonal string modifications of the Kerr solutions
Various classes of exact solutions in heterotic string gravity can be constructed, which depend on which type of effective sources (19) and (38) are chosen. As a result, various classes of generic off-diagonal deformations of the Kerr metric into higher-dimensional exact solutions can be constructed. In this subsection we shall construct and analyze a series of 6-d and 10-d solutions encoding possible higher-dimension interactions with effective cosmological constants, warping configurations, f -modifications and certain analogies to almost-Kä hler gravity models.

6-d deformations with nontrivial cosmological constant
Solutions are determined by certain configurations of the NSfields H α 1 β 1 μ 1 which are nontrivial on the first shell s = 1, see Eqs. (21) and the respective term in (38), where the coefficient 3 is used for an effective 6-d space time with trivial extension on four other internal space coordinates. We introduce an effective source in the quadratic element (96) extended with one nontrivial extra 2-d shell.
Such a family of generic off-diagonal stationary solutions is described by × dy 6 + 1+ε In the above formulas, we consider a generating function H 1 χ( x i , y 4 , y 6 ), which results, in general, in nontrivial nonholonomic torsions encoding contributions of H α 1 β 1 μ 1 via an effective cosmological constant ( H s) 2 . One considers a summation on i 1 = 1, 2, 4 for terms like Off-diagonal extra-dimensional gravitational interactions modify a Kerr metric for any nontrivial cosmological constant determined by a corresponding ansatz for H α 1 β 1 μ 1 in 6-d. In a similar form we can generalize the constructions in 8-d and 10-d with nontrivial extra shell components H α s β s μ s and R = 0 in (19).

10-d deformations with NS 3-form and 6-d almost Kähler internal spaces
The class of solutions (103) can be generalized for 10-d spacetimes with indices of fields running over 10-d values and nontrivial sources (34) redefined in (38), with nontriv-ial H = −5( H s) 2 , F = −5n F ( F s) 2 and R = −5trn R ( R s) 2 in = H + F + R . The corresponding quadratic line element is × dy 5 + 1 + ε 1 n i 1 dy 6 1 χ + 3 χ g 9 (x k ( x k ))[dy 9 + 1 + ε 3 n i 3 × dy 10 3 χ + This generic off-diagonal stationary metric for 10-d spacetimes defines a very general class of stationary solutions of the nonholonomic equations of motion in heterotic gravity (12)- (15) with effective scalar field encoded into the Nconnection structure. It is important to use "shell by shell" nonholonomic variables for the quadratic element (104).
Only in such cases we can understand the nonlinear symmetries and classify the types of generating and integration functions corresponding to horizontal and higher order vertical conventional sources. Such properties cannot be encoded in a minimal form if notations for indices running coordinate values from 0 to 9 are used as in standard papers on heterotic supergravity. For small ε-deformations, the 4-d component of such metrics describes a Kerr type black ellipsoid with eccentricity ε and ϕ-anisotropic polarization of physical constants and horizons. The nonholonomic deformations also encode sources from all gauge-like and interior space gravitational fields via a redefinition of generating functions. Similar classes of non-vacuum solutions can also be modeled for Einstein-Finsler spaces if extra-dimensional coordinates are treated as velocity/momentum coordinates [10,24,26,43,44,53,54]. The metrics possess a respective Killing symmetry on ∂ t and ∂/∂y 9 . Using the AFDM, we can construct solutions depending on all 10-d coordinates which may describe geometric evolution and time propagating Kerr black holes determined by heterotic string gravity effects. They define ε-deformations of Kerr-de Sitter black holes into ellipsoid configurations with effective cosmological constants determined by constants in GR, possible fmodifications and extra-dimension contributions [72]. By nonholonomic frame transforms and distortions of the canonical d-connection, the above metric can be rewritten in canonical almost-Kähler variables on 6-d internal space, as we prove in [48]. We omit such constructions (being very important in various models of deformation and brane quantization of MGTs and geometric evolution theories [10,11,63,70]) in this work.

Outlook and concluding remarks
In this work, we have applied the anholonomic frame method, AFDM, for constructing new classes of stationary solutions of motion equations in heterotic supergravity. Such solutions have generic off-diagonal metrics for effective ten-dimensional, 10-d, spacetimes enabled with generalized connections and depend on all possible 4-d and extradimensional space coordinates. They admit subclasses of solutions with warping on coordinate y 4 , nearly almost-Kähler 6-d internal manifolds in the presence of nonholonomically deformed gravitational and gauge instantons. The almost-K ähler structure is necessary if we want to generate the Kerr metric with possible (off-) diagonal and nonholonomic deformations to black ellipsoid type solutions with locally anisotropic polarized physical constants, small deformations of horizons, embedding into nontrivial extradimension vacuum gravitational fields and/or gauge configurations. These solutions preserve two real supercharges corresponding to N = 1/2 supersymmetry from the viewpoint of four non-compact dimensions and various nonholonomic deformations.
Following the AFDM, we can integrate the equations of motion in heterotic supergravity in very general forms with dependence on all 10-d spacetime coordinates. Such constructions are possible to more general classes of almost-Kähler nonholonomic variables and so-called canonical nonholonomic variables. This allows us to decouple string modified Einstein equations with effective sources in general form which are similar to Einstein-Yang-Mills-Higgs, EYMH, systems in higher-dimensions and with generalized gaugelike interactions. It is possible to consider associated SU (3) structures and solve generalized BPS equations and Bianchi identities. A crucial difference from former approaches is that our geometric methods allow us to work with generating and integration functions for off-diagonal metrics and connections transforming equations of motion into nonlinear systems of partial differential equations, PDEs. In particular, we can reproduce former results for a diagonalizable ansatz with dependence on radial and warping coordinates as solutions of ordinary differential equations, ODEs.
To illustrate the power and importance of the AFDM as the most general geometric method of constructing analytic solutions of (modified) motion/gravitational and field equations, we show how this formalism can be applied for generating Nadapted (i.e. adapted to nonlinear connection structures) YM and instanton configurations with possible associated SU (3) nonholonomic structures; see the associated work [48]. New classes of exact solutions describing small parametric modifications of Kerr metrics with effective string sources are provided. We show that in a certain sense, a large class of physical effects in modified gravity models like R 2 can be equivalently modeled/explained by nonholonomic constraints and effective sources in heterotic string gravity. In explicit form, exact/parameteric extra dimension deformations of the black hole metrics in 6-d and 10-d gravity with NS-3 form and 6-d almost-Kähler internal spaces are constructed and analysed.
Finally, we emphasize that there are a plethora of future directions which can be pursued using our methods and results as starting points. This includes the construction of cosmological solutions in the heterotic string gravity and/or the study of smooth compact nonholonomic varieties in both heterotic and geometric flow context. Similar analyses can be performed in type II string theory in particular, including the Ramond-Ramond sector and/or considering geometric flows on internal spaces. We worked to the lowest order of α but possibilities in the AFDM exist to extend the constructions to higher orders.