Off-Diagonal Deformations of Kerr Metrics and Black Ellipsoids in Heterotic Supergravity

Geometric methods for constructing exact solutions of motion equations with first order $\alpha ^{\prime}$ corrections to the heterotic supergravity action implying a non-trivial Yang-Mills sector and six dimensional, 6-d, almost-K\"ahler 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)-d and/or (1+3)-d domain wall configurations, with possible warping nearly almost-K\"ahler 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 4-d 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 gravity 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 -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-Worker (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 motion equations are transformed into certain systems of nonlinear second order ordinary equations (ODE), 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 motion equations 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 can not 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 mater 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 3-4 spacetime coordinates transforms the Einstein equations into systems of nonlinear PDEs, which can not 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 higher dimensional 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 form, certain classes of N-elongated frames of reference, 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 motion equations 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 zerotorsion conditions at the end, i.e. after we have found a class of generalized solutions. We can not 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 nonhlonomic 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 section 4.4 [1] 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 back-grounds 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 modelled 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 ( [24,51,8]) 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 section 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 higher dimension spacetimes is outlined. We develop the AFDM as a geometric method for constructing exact solutions of motion equations 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 4-d 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-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 4-d, 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 motion equations in heterotic string theory.
Finally (in section 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 modelled in the low-energy limit of heterotic string theory as a N = 1 and 10-d supergravity coupled to super Yang-Mills theory, see details in [48]. The nonholonomic variables will be parameterized in forms which allow a general decoupling of motion equations and generating exact solutions depending, in principle, on all spacetime coordinates. Such a higher dimension spacetime is modelled 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,44,42,43]. 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 lables and shell indices. This minimizes the procedure of separating the motion equations 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 works, it is not possible to understand how the AFDM can be applied for generating off-diagonal solutions in (super) string/gravity theories.

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 follow an approach re-defined for nonholonomic distributions on (super) manifolds [43,44,40]. 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.
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, Nconnection, 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 as is when s N = N as is ( s u)dx is ⊗ ∂/∂y as . 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 N-adapted local bases, with N-elongated partial derivatives, e νs = (e is , e as ), and cobases with N-adapted differentials, e µs = (e is , e as ), For s = 0, ∂ ∂y a 0 , e i 0 = dx i , e a = e a 0 = dy a + N a i dx i ; or/ and (3) ∂y as , e as = ∂ ∂y as , e is = dx is , e as = dy as + N as is dx is for s = 1, 2, 3.
The anholonomy relations [e αs , e βs ] = e αs e βs − e βs e αs = W γs αsβs e γs (4) are computed W bs isas = ∂ as N bs is and W as jsis = Ω as isjs , where the curvature of N-connection is defined as the Neijenhuis tensor, Ω as isjs := e js N as is − e is N as js . Any metric structure on M can be parametrized via nonholonomic frame transforms as a distinguished metric (d-metric, in boldface form), s g = {g αsβs }, s g = g isjs ( s u) e is ⊗ e js + g asbs ( s u)e as ⊗ e bs (5) = g ij (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 asbs ( s u)e as ⊗ e bs .
Introducing dual coframes e µs = (e is , e as ) decomposed with respect to N-elongated differentials du as 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 liniear 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 s D completely defined by s g following certain geometric principles. The canonical d-connection s D is characterized by a nonholonomically induced torsion d-tensor which is completely defined by s g for any chosen s N = {N as is }. The N-adapted coefficients are parameterized "shell by shell" by formulas −Ω as jsis , T cs asjs = L cs asjs − e as (N cs js ), T as bscs = C as bscs − C as csbs }.
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 LC-configurations with zero torsion if we additionally impose for (7) the conditions T γs 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 details in [48]. The Ricci d-tensor Ric = { R αsβs := R τs αsβsτs } is a respective contracting of the coefficients of curvature tensor, R αsβs = { R hsjs := R is hsjsis , R jsas := − P is jsisas , R bsks := P as bsksas , R bscs = S as bscsas }.
Consindering the inverse d-metric to s g, we define and compute the scalar curvature of s D, s R := g αsβs R αsβs = g isjs R isjs + g asbs R asbs with respective h-and v-components of scalar curvature, R = g ij R ij , S = g ab S ab ,

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 noholonomic 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 'Alambert 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 N-adapted operator s A D and curvature F = F( 1 ψ) via a map constructed in [1,48]).
The equations (12) can be written as effective Einstein equations for the canonical d-connection s D, L cs asjs = e as (N cs js ), C is jsbs = 0, Ω as jsis = 0.
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 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 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 effective matter fields φ, H αsβsµs , interior space fields R µsαsβsγs and gauge fields F µsαs .
The nonholonomic motion equations 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 above sources, and respective effective cosmological constants, are related via condition (13).

Decoupling & integration of nonholonomic motion equations
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-confiugrations with s ∇ will be analysed 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 necessary shells, we can extend the constructions for non-Killing configurations. Such a procedure is described for arbitrary finite shell s in Refs. [43,40,41,42]. 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 where e a = dy a + N a . 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 as = h(y 4 )ȟ as (x i , y bs ) 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. 5 we shall put a left lable K in order to emphasize that this is a d-metric with Killing symmetry 2.3.2
This imposes certain N-adapted symmetries on the Einstein d-tensor for the ansatz (24), see details in [43,40].

N-adapted sources and nonholonomically modified Einstein equations
We will be able to integrate nonholonomic motion equations 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 5 independent effective sources, respectively, on h-and s v-subspaces, which will be related to certain effective cosmological constants corresponding to 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 non-zero, for instance, φ h Λ = 0 but H 2 Λ = 0. This depends on the type of vacuum, non-vacuum, or effective vacuum model we study. It is possible to compensate 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 d-metrics 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 0 φ(x i , y a , y as ) . 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 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 motion equations (12)-(15) written with respect to N-adapted frames if the conditions (36) are satisfied. Nontrivial coefficents with φ and its partial derivatives will appear if certain physical equations are written, for instance, in coordinate frames. We have chosen special type 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 H Υ µsνs , F Υ µsνs and R Υ µsνs for which Υ µsνs = H Υ µsνs + F Υ µsνs + R Υ µsνs , can be nonholonomicaly transformed into an N-adapted diagonal form Nonholonomic deformations of string gauge fields can be parametrized in the form when, 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 effective source. Such parametrizations are possible by considering the ansatz with effective cosmological constants H αsβsµs = H s |g βsµs |ǫ αsβsµs , F µsαs = F s |g βsµs |ǫ µsαs , R µsαsβsγs = R s |g βsµs |ǫ µsαsβsγs , with absolute anti-symmetric ǫ-tensors (we refer readers to similar details for 4-d to [42] and references therein). For such an ansatz, we obtain effective energy-momentum tensors were we considered formulas (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 motion equations and effective gravitational field equations
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 is instead of s w τ s−1 , s n is 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 Φ-vaules) and α-, β-, γ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).
The 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 [43,40]. In brief the decoupling property of the AFDM can be explained for 4-d configurations: 1. The 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 ̟ in (48), four values h 3 (x i , y 4 ), h 4 (x i , 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 ) another two such functions. We note that h a are coefficients of a d-metric, ̟ is a generating function and Υ is related to extra dimensional string contributions 3. Using h 3 and ̟, or Φ, from the previous point, we compute the coefficients α i and β, see (48). This allows 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).

2.3.6
Integration of nonholonomic motion equations 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) transform respectively into 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) as 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 h 3 (x k ) is an integration function. Considering (55), (53) and formula (56), we compute The first part of the N-connection coefficients are found by integrating two times on y 4 in (51) written in the form n * * i = (n * i ) * = − n * i (ln |h 3 | 3/2 /|h 4 |) * 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 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), 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 details in [14,70].
Examples of 4-d -10-vacuum configurations: The configurations with zero source can not 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)-(32) are zero but the Ricci and torsion d-tensor are not trivial. The first equation is a typical example of 2-d Laplace equation with solutions expressed in the form g i = e q(x k ,Υ=0) .
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 4-d spacetimes.
2. In such cases, h * 3 = 0 and h * 4 = 0. It is possible to solve the equation (26) and (40) in a self-consistent form for Υ = 0 if ̟ * = 0 for respective coefficients in (48). Fixing ̟ = ̟ 0 = const, we can consider arbitrary functions w i (x k , y 4 ) because β = α i = 0 for such configurations. The conditions (48) are satisfied by any where 0 h 3 (x k ) is an integration function and h 3 (x k , y 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 functions 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 zero effective source and adding respective generating functions gener. functs. h 9 (x i , y 4 , y a1 , y a2 , y 10 ) 3 w i3 (x i , y 4 , y a1 , y a2 , y 10 ) .
Finally, we note that we can construct generic off-diagonal vacuum solutions in heterotic supergravity with different types on different shells 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 4-d 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 time like coordinate y 3 = t and do not depend on the coordinate y 5 . For non-trivial 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 5 (x k , y 4 ) 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 1 n * 1 * 1 i 1 = (n * 1 i 1 ) * 1 = − n * 1 i 1 (ln |h 5 | 3/2 /|h 6 |) * 1 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 ds 2 K10d = g α 2 β 2 (x k , y 4 , y a 1 , y a 2 , y a 2 )du α 3 du β 3 = ds 2 K8d [ see (66)] + [h [dy 9 + 3 1 n k 3 + 3 2 n k 3 dy 10 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 motion equations 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 above classes of solutions, used in differen forms in [42,43]. It allows us to re-define 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 (68) up to certain classes of integration functions depending on coordinates x i , we express the solutons (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 equvalent forms the N-adapted 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, The nonlinear symmetry reflects the property of changing the effective sources mentioned in (34), We can indentify Λ 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 ) → 3 Λ = φ 3 Λ+ H 3 Λ+ F 3 Λ+ R 3 Λ 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 re-define the constructions equivalently, via new classes of generating functions, for effective cosmological constants.
Finally, we note that formula (67) simplifies for 3 , 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 Nconnection and d-metric structure. The zero torsion conditions (8) can be satified by a subclass of nonholonomic distributions determined by corresponding parametrizations of the generating and integration functions and sources. By straightforward computations (see details in Refs. [38,35,39,40]), we can verify that if the coefficients of N-adapted frames and s v-components of d-metrics are subjected to respective conditions, The n-coefficients solve the conditions (71) if . The explicit form of solutions of constraints on w k derived from (71) depend on the class of vacuum or non-vacuum metrics we try to construct, see details in [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) in explicit form if Υ = const, or if such an effective source can be expressed as a functional Υ(x i , y 4 ) = Υ[Φ]. The third condition for s = 0, ∂ i w j = ∂ j w i , can be satisfied for any generating functionǍ =Ǎ(x k , y 4 ) for which w i =w i = ∂ iΦ /∂ 4Φ = ∂ iǍ . Following similar considerations for shells s = 1, 2, 3, we formulate the LC-conditions for generating functions Υ(x i , y 4 ) = Υ[Φ], or Υ = const; s = 1 : Imposing respective conditions from (72) and (73) on coefficients of (58), we define such a class of quadratic elements for off-diagonal 4-d stationary configurations in heterotic supergravity with zero induced torsion, We can impose similar conditions and generate exact off-diagonal 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 4-d 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 solutons 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 typeg = [g i ,h as ,N js bs ] when +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 ) (e a 3 ) 2 , wheree as as N-elongated differentials. Such a metric is diagonalizable if there is a coordinate transform To construct exact solutions with non-singular coordinates it is important to work with "formal" off-diagonal parameterizations when the coefficients sẘ is and/or sn is are not zero but the anholonomy coefficientsW α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 ors Schwarzshild 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 e 3 = dt + n η ini dx i , e 4 = dy 4 + w η iẘi dx i , e 5 = dy 5 + n η i 1n i 1 dx i 1 , e 6 = dy 6 + w η i 1ẘ i 1 dx i 1 , e 7 = dy 7 + n η i 2n i 2 dx i 2 , e 8 = dy 8 + w η i 2ẘ is dx is , e 9 = dy 9 + n η i 3n i 3 dx i 3 , e 10 = dy 10 where the coefficients (g αs = η αsgαs , w η isẘis , n η is n is ) define, for instance, a d-metric s g (24) as a solution of nonholonomic motion equations 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 express using d-metric and N-connection deformations stated in explicit form on all shells for the d-metric and n-and w-coefficients: For the coefficients of d-metrics, y 4 ), η a 1 = 1 + εχ a 1 (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 coeffisients of N-connection, n η i = 1 + ε n η i (x k , y 4 ), w η i = 1 + ε w χ i (x k , y 4 ), n η i 1 = 1 + ε n η i 1 (x k , y 4 , y 6 ), w η i 1 = 1 + ε w χ i 1 (x k , y 4 , y 6 ), n η i 2 = 1 + ε n η i 2 (x k , y 4 , y a 1 , y 8 ), w η i 2 = 1 + ε w χ i 2 (x k , y 4 , y a 1 , y 8 ), n η i 3 = 1 + ε n η i 3 (x k , y 4 , y a 1 , y a 2 , y 10 ), w η i 3 = 1 + ε w χ i 3 (x k , 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 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 as , ε N a s−1 bs ). 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 detalis for 4-d configurations and then extended to higher dimensional 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 v-components 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 2-hypersurface 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 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 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 formulas (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 necessary formulas).
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 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 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 motion equations in heterotic supergravity: where n i (x k ) is a re-defined 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), ε h 3+2s = [1 + ε χ 3+2s ]g 3+2s = 1 − ε 1 4g 3+2s dy 4+2s ( sΦ2 s χ) * s s Υ g 3+2s ; ε h 4+2s = [1 + ε χ 4+2s ]g 4+2s = 1 + ε 2( s χ + sΦ sΦ * s s χ * s ) + 1 4g 3+2s dy 4+2s ( sΦ2 s χ) * 1 s Υ g 4+2s ; These formulas with s = 1, 2, 3 allow us to parametrize all coefficients of vertical components of d-metrics and N-connections. 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 references therein) for any stable prime solution.
We can impose additional constraints in order to extract LC-configurations as we considered in section (2.3.8).

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 motion equations 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 solutions when the Kerr metrics are nonholonomically deformed into general or ellipsoidal stationary configurations in four dimensional gravity and/or extra dimensions are considered. We provide and study examples of generic off-diagonal metrics encoding nonlinear interactions with 3-form and gauge like fields, nonholonomically induced torsion effects and instanton configurations.

Preliminaries on the Kerr vacuum solution and nonholonomic variables
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 and y 3 = t + y 3 (x 1 , x 2 ), y 4 = ϕ + y 4 (x 1 , 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 ).
We send readers to [71,31,52] for main results and methods for stationary black hole solutions (the coordinates x 1 , x 2 introduced above correspond to x, y respectively from chapter 4 of the first book). The primed quadratic linear elements (85) (or (86)) (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-Linquist 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.

solutions of vacuum Einstein equations parametrized in the form
Working with general classes of stationary solutions generated by the AFDM, the Kerr vacuum solution in GR can be considered as a "degenerated" 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 nonolonomically 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 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 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 extra dimensional string sources with ansatz resulting in cosmological constants described in terms of polarization functions, where |η 4 ′ |, may have not a smooth limit to 1, where the coefficient 1/4 was introduced inΦ and the functionǍ is any value for which ∂ i ′Φ/Φ * = ∂ i ′Ǎ.
Summarizing the above formulas, we obtain the quadratic element We used an 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 are 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 section 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 subjected 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 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 [38,55,56,24,51] 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 formulas (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 ε g i =g i [1 + εχ i ] = [1 + εe 0 q 1 q/ K Λ]g i 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 coeffisients of N-connection, n η i = 1 + ε n η i (x k , y 4 ), w η i = 1 + ε w χ i (x k , y 4 ).
Summarizing the above for 4-d ε−configurations, we obtain the quadratic element In general, these heterotic string deformations of the Kerr metric are nonholonomically induced torsion coefficients, linear on ε.
We can impose additional constraints on χ, n i andΦ which allows us to extract LC-configurations. The corresponding ε-deformed analogue of the metric (90) can be written where χ defines a deformed generating function ε Φ =Φ(x k , y 4 )[1 + εχ(x k , y 4 )] as in formula (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 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 d-connection and/or the Levi-Civita connection. Generic off-diagonal metrics of type (92) can be parameterized as ellipsoidal deformations of the Kerr metric on a small eccentricity parameter ε.
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 can not 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 both the event horizon, r + = m(ϕ) + m 2 (ϕ) − a 2 sin 2 ϑ 1/2 , and the Cauchy horizon, which are effectively imbedded 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 re-definition of generating functions (see formulas (70) and (37)) adapted to εdeformations. The corresponding quadratic line elements are The new classes of ε-deformed solutions determine Kerr-like black hole solutions with additional dependencies on variable ϕ of certain diagonal and off-diagonal coefficients of metric. There is an obvious anisotropy in 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 details in [38], see also references inside on 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 The parametric equation of an ellipse with radial parameterr + = 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 (101) where n i (x k ) is a re-defined 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 sting 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 modelled 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 NS-fields H α 1 β 1 µ 1 which are nontrivial on the first shell s = 1, see formulas (21) and respective term in (38) where the coefficient 3 is used for an effective 6-d space time with trivial extension on 4 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 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 10-d values and nontrivial sources (34) re-defined in (38), with nontrivial 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 This generic off-diagonal stationary metric for 10-d spacetimes defines a very general class of stationary solutions of the nonholonomic motion equations in heterotic gravity (12)-(15) with effective scalar field encoded into the N-connection 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 can not 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 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 re-definition of generating functions. Similar classes of non-vacuum solutions can also be modelled for Einstein-Finsler spaces if extra dimensional coordinates are treated as velocity/momentum coordinates [10,24,26,53,54,43,44]. 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 propogating Kerr black holes determined by heterotic sting gravity effects. They define ε-deformations of Kerr -de Sitter black holes into ellipsoid configurations with effective cosmological constants determined by constants in GR, possible f -modifications 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.
The above quadratic elements transform into a trivial embedding of the 4-d Kerr solution into a 10-d spacetime if ε → 0 and primary data (labelled by circles) can be transformed into a diagonal 6-d internal space. The solution (105) is equivalent to (104) up to re-definition of coordinates and some integration functions. It is difficult to see "shell by shell" nonlinear symmetries and construct generic off-diagonal solutions in the variant with standard variables used in heterotic string gravity. In former non N-adapted variables, the formulas are much more cumbersome and less adapted for studying limits to well known solutions and generalizations with nontrivial backgrounds and extra dimensional contributions. To elaborate and apply a corresponding "shell by shell" N-adapted geometric techniques of constructing exact solutions in 4d and extra dimension theories is important both from a mathematical stand point and a physical point of view.

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 dimenisonal, 10-d, spacetimes enabled with generalized connections and depend on all possible 4-d and extra dimensional 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 extra dimension 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 motion equations 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 gauge like 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 motion equations 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 N-adapted (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 modelled/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.