On Exact Solutions and Perturbative Schemes in Higher Spin Theory

We review various methods for finding exact solutions of higher spin theory in four dimensions, and survey the known exact solutions of (non)minimal Vasiliev's equations. These include instanton-like and black hole-like solutions in (A)dS and Kleinian spacetimes. A perturbative construction of solutions with the symmetries of a domain wall is described as well. Furthermore, we review two proposed perturbative schemes: one based on perturbative treatment of the twistor space field equations followed by inverting Fronsdal kinetic terms using standard Green's functions; and an alternative scheme based on solving the twistor space field equations exactly followed by introducing the spacetime dependence using perturbatively defined gauge functions. Motivated by the need to provide a higher spin invariant characterization of the exact solutions, aspects of a proposal for a geometric description of Vasiliev's equation involving an infinite dimensional generalization of anti de Sitter space is revisited and improved.


Introduction
Higher spin (HS) theory in four dimensions, in its simplest form and when expanded about its (anti-)de Sitter vacuum solution, describes a self-interacting infinite tower of massless particles of spin s = 0, 2, 4.... The full field equations, proposed long ago by Vasiliev [1,2,3] (for reviews, see [4,5]), are a set of Cartan integrable curvature constraints on master zero-, one-and two-forms living on an extension of spacetime by a non-commutative eight-dimensional twistor space. The latter is fibered over a four-dimensional base, coordinated by a Grassmann-even SL(2, C)-spinor oscillator Z A = (z α ,zα), and the fiber is coordinatized by another oscillator Y A = (y α ,ȳα); the master fields are horizontal forms on the resulting twelve-dimensional total space, valued an infinite-dimensional associative algebra generated by Y A , that we shall denote by A, and subject to boundary conditions on the base manifold.
A key feature of Vasiliev's equations is that they admit asymptotically (anti-)de Sitter solution spaces, obtained by taking the HS algebra A to be an extension of the Weyl algebra, with its Moyal star product, by involutory chiral delta functions [6,7], referred to as inner Klein operators, relying on a realization of the star product using auxiliary integration variables [4]. Introducing a related class of forms in Z-space, that facilitates a special vacuum two-form in twistor space, the resulting linearized master fields can be brought to a special gauge, referred to as the Vasiliev gauge, in which their symbols defined in a certain normal order are real analytic in twistor space, and the master zero-and oneforms admit Taylor expansions in Y at Z = 0 in terms of Fronsdal fields on the mass shell and subject to physical boundary conditions. Although Vasiliev equations take a compact and elegant form in the extended space, their analysis in spacetime proceeds in a weak field expansion which takes an increasingly complicated form beyond the leading order. Indeed, they have been determined so far only up to quadratic order. In performing the weak field expansion, a number of challenges emerge. Firstly, obtaining these equations requires boundary conditions in twistor space, referring to the topology of Z space and the classes of functions making up A [2,8,9]. Proper way to pin down these aspects remains to be determined. Second, the cosmological constant, Λ, which is necessarily nonvanishing in Vasiliev's theory (as the transvection operators of the isometry algebra are realized in A as bilinears in Y ), appears in the effective equations to its first power via critical mass terms, but also to arbitrary negative powers via non-local interactions [4,10]. Thus, letting φ denote a generic Fronsdal field, it follows that ∂φ ∼ |Λ|φ on-shell, and hence interactions with any number of derivatives are of equal relevance (at a fixed order in weak field amplitudes). This raises the question of just how badly nonlocal are the HS field equations, the attendant problem of divergences arising even at the level of the amplitudes [11,12] and what kind of field redefinitions are admissible. One guide available is the holographic construction of the bulk vertices [13,14,15]. Clearly, it would be desirable to find the principles that govern the nonlocal interactions, based on the combined boundary conditions in twistor space as well as spacetime, such that an order by order construction of the bulk vertices can proceed from the analysis of Vasiliev equations. The simple and geometrical form of Vasiliev equations, in turn, may pave the way for the construction of an off-shell action that will facilitate the computation of the quantum effects.
In an alternative approach to the construction of HS equations in spacetime, it has been proposed to view Vasiliev's equations as describing stationary points of a topological field theory with a path integral measure based on a Frobenius-Chern-Simons bulk action in nine dimensions augmented by topological boundary terms, which are permitted by the Batalin-Vilkovisky formalism, of which only the latter contributes to the on-shell action [16,17].
This approach combines the virtues of the on-shell approach to amplitudes for massless particles flat spacetime with those of having a background independent action, in the sense that the on-shell action is fixed essentially by gauge symmetries and given on closed form, which together with the background independence of Vasiliev's equations provides a machinery for perturbative quantum computations around general backgrounds.
In this context, it is clearly desirable to explore in more detail how the choice of boundary conditions in the extended space influences the classical moduli space of Vasiliev's equations, with the purpose of spelling out the resulting spaces, computing HS invariant functionals on-shell, and examining how the strongly coupled spacetime nonlocalities are converted into physical amplitudes using the aforementioned auxiliary integral representation of star products in twistor space.
The aim of this article is to review three methods that have been used to find exact solutions of Vasiliev equations, and to describe two schemes for analyzing perturbations around them. In particular we will describe the gauge function method [18,19] for finding exact solutions and summarize the first such solution found in [20], as well as its generalization to de Sitter spacetime studied in [21] together with the solutions of for a chiral version of the theory with Kleinian (2, 2) signature. As we shall see, this method uses the fact that the spacetime dependence of the master fields can be absorbed into gauge functions, upon which the problem of finding exact solutions is cast into relatively manageable deformed oscillator problem in twistor space. The role of different ordering schemes for star product as well as gauge choices to fix local symmetries in twistor space will also be discussed.
Next, we will describe a refined gauge function method proposed in [22], where the twistor space equations are solved by employing separation of twistor variables and holomorphicity in the Z space in a Weyl ordering scheme and enlarging the Weyl algebra in the fiber Y space by inner Kleinian operators. This approach provides exact solution spaces in a particular gauge, that we refer to as the holomorphic gauge, after which the spacetime dependence is introduced by means of a sequence of large gauge transformations, by first switching on a vacuum gauge function, taking the solutions to what we refer to as the L-gauge, where the configurations must be real analytic in Z space, which provides an admissibility condition on the initial data in holomorphic gauge. The solutions can then be mapped further to the Vasiliev gauge, where the linearized, or asymptotic, master fields, are real analytic in the full twistor space and obey a particular gauge condition in Z space which ensures that they consist of decoupled Fronsdal fields in a canonical basis; the required gauge transformation, from the L gauge to the Vasiliev gauge, can be implemented in a perturbation scheme, which has so far been implemented mainly at the leading order. We will describe a black hole-like solution in some detail and mention other known solutions obtained by this method so far, including new solutions with six Killing symmetries [23].
We shall also outline a third method, in which the HS equations are directly tackled without employing gauge functions. In this method, solving the deformed oscillators in twistor space also employs the projector formalism, though the computation of the gauge potentials does not rely on the gauge function method. The black hole-like solution found in this way in [24] will be summarized.
We shall also review two approaches to the perturbative treatment of Vasiliev's equations. One of them, which we refer to as the normal ordered scheme, is based on a weak field expansion around (anti-)de Sitter spacetime [3,4,25]. It entails nested parametric integrals, introduced via a homotopy contraction of the de Rham differential in Z space used to solve the curvature constraints that have at least one form index in Z space, followed by inserting the resulting perturbatively defined master fields into the remaining curvature constrains with all form indices in spacetime. In an alternative scheme, the equations are instead solved exactly in at he aforementioned L-gauge, and perturbatively realized large HS gauge transformation is then performed to achieve interpretation in terms of Fronsdal fields in asymptotically (anti-)de Sitter spacetimes in Vasiliev gauge [7]. The advantages of the latter approach in describing the fluctuations around more general HS backgrounds will be explained.
A word of caution is in order concerning the usage of 'black hole' terminology in describing certain types of exact solutions to HS equations. This terminology is, in fact, misleading in some respects since the notion of a line interval associated with a metric field is not HS invariant. Indeed, the apparent singular behaviour at the origin may in principle be a gauge artifact. This point is discussed in more detail in Section 4.2. Moreover, given the nonlocal nature of the HS interactions, the formulation of causality, which is crucial in describing the horizon of a black hole, is a challenging problem without any proposal for a solution yet in sight; in fact, a more natural physical interpretation of the black hole-like solutions may turn out to be as smooth black-hole micro states [7,26]. Another aspect of the known black hole-like solution in HS theory is that they activate fields of all possible spins, and apparently it is not possible to switch of all spins except one even in the asymptotically AdS region.
So, what is meant by a black hole solution in HS theory? Firstly, the SO(3) × SO(2) symmetry of the solution (which is part of an infinite dimensional extended symmetry forming a subgroup of the HS symmetry group) is in common with the symmetry group arising in the asymptotically AdS BH solution of ordinary AdS gravity. Second, the solution contains a spin-two Weyl tensor field which takes the standard Petrov type D form, with a singularity at the origin; more generally, the spin-s Weyl tensors are of a generalized Petrov type D form, given essentially by the s-fold direct products of a spinone curvature of the Petrov type D form. The BH terminology is thus used in the context of HS theory with the understanding that it is meant to convey these properties, albeit they do not constitute a rigorous definition of a black hole in HS theory.
The use of HS invariants for exact solutions to capture their physical characteristics has been considered and in some cases they have been computed successfully. These particular invariants alone do not, however, furnish an answer to the question of whether it makes any sense to think about event horizons in HS theory at all, and if so, how to define them; in fact, their existence rather supports the aforementioned micro state proposal, wherein the HS invariants can be interpreted as extensive charges defining HS ensembles.
Motivated by the quest for giving a physical interpretation of the exact solution in the context of underlying HS symmetries, a geometrical approach to HS equations was proposed in [27]. We shall summarize this proposal in which the HS geometry is based on an identification of infinite dimensional structure group in a fibre bundle setting, and the related soldering phenomenon that leads to a HS covariant definition of classes of (non-unique) generalized vielbeins and related metrics, and as such an infinite dimensional generalization of AdS geometry. In doing so, we will improve the formulation of [27] by dispensing with the embedding of the relevant infinite dimensional coset space into a larger one that involves the extended HS algebra that includes the twistor space oscillators.
Finally, we are not aware of any exact solutions of HS theories in dimensions D > 4 [28,29], while in D=3, assuming that the scalar field is coupled to HS fields, we can refer to [8,30,31] for the known solutions. Purely topological HS theory, which has no dynamical degrees of freedom, and which allows a more rigorous definition of black holes, is known to admit many exact solutions whose description goes beyond the scope of this review.

Bosonic model in (A)dS
Vasiliev's theory is formulated in terms of horizontal forms on a non-commutative fibered space C with four-dimensional non-commutative symplectic fibers and eight-dimensional base manifold equipped with a non-commutative differential Poisson structure. On the total space, the differential form algebra Ω(C) is assumed to be equipped with an associative degree preserving product ⋆, a differential d, and an Hermitian conjugation operation †, that are assumed to be mutually compatible. The base manifold is assumed to be the direct product of a commuting real four-manifold X 4 with coordinates x µ , and a non-commutative real four-manifold Z 4 with coordinates Z A ; the fiber space and its coordinates are denoted by Y 4 and Y A , respectively. The non-commutative coordinates are assumed to obey where C AB is a real constant antisymmetric matrix. The non-commutative space is furthermore assumed to have a compatible complex structure, such that where the complex doublets obey [y α , y β ] ⋆ = 2iǫ αβ and [z α , z β ] ⋆ = −2iǫ αβ . The horizontal forms can be represented as sets of locally defined forms on X 4 × Z 4 valued in oscillator algebras A(Y 4 ) generated by the fiber coordinates glued together by transition functions, that we shall assume are defined locally on X 4 , resulting in a bundle over X 4 with fibers given by Ω(Z 4 )⊗ A(Y 4 ). The algebra A(Y 4 ) can be given in various bases; we shall use the Weyl ordered basis, and the normal ordered basis consisting of monomials in a ± = Y ± Z with a + and a − oscillators standing to the left and right, respectively. We assume that the elements in Ω(Z 4 ) ⊗ A(Y 4 ) have well-defined symbols in both Weyl and normal order. The normal order reduces to Weyl order for elements that are independent of either Y or Z, and in the cases where depend on both Y and Z, they can be composed using the Fourier transformed twisted convolution formula in normal ordered scheme as The model is formulated in terms of a zero-form Φ, a one-form and a non-dynamical holomorphic two-form with Hermitian conjugate J = (J) † , where b is a complex parameter and κ := κ y ⋆ κ z , κ y := 2πδ 2 (y) , κ z := 2πδ 2 (z) , (2.7) are inner Klein operators obeying κ y ⋆ f ⋆ κ y = π y (f ) and κ z ⋆ f ⋆ κ z = π z (f ) for any zero-form f , where π y and π z are the automorphisms of Ω(Z 4 ) ⊗ A(Y 4 ) defined in Weyl order by π y (x; z,z; y,ȳ) = (x; z,z; −y,ȳ) , π z (x; z,z; y,ȳ) = (x; −z,z; y,ȳ) .
It follows that dJ = 0, J ⋆ f = π(f ) ⋆ J and π(J) = J, idem J, with π := π y • π z ,π := πȳ • πz . (2.9) It is useful to note that the inner Kleinian takes the following forms in different ordering schemes: κ = e iy α zα in normal ordering scheme (2π) 2 δ 2 (y)δ 2 (z) in Weyl ordering scheme (2.10) The nonminimal and minimal models with all integer spins and only even spins, respec- where τ is the anti-automorphism It follows that τ (J,J ) = (−J, −J ). Models in Lorentzian spacetimes with cosmological constants Λ are obtained by imposing reality conditions as follows [21]: Basic building blocks for Vasiliev equations are the curvature and twisted-adjoint covariant derivative defined by respectively, where the π-twisted star commutators is defined as and for parameters obeying the same kinematic conditions as the connection.
It remains to be a challenging problem to determine if these equation can be derived for a suitable tensionless, or critical tension, limit followed by a consistent truncation of string field theory on a background involving AdS 4 . It will be very interesting to also determine if this equations follow from a consistent quantization of a topological string field theory. For a more detailed discussion and progress in this direction, see [32,33]. The component fields of V A do not transform properly under the Lorentz transformations generated by To remedy this problem and achieve manifest Lorentz covariance, one introduces the field-dependent Lorentz generators [4,25] M αβ = M 20) and their complex conjugates, where where (ω αβ µ , ωαβ µ ) is the canonical Lorentz connection. It is defined up to tensorial shifts [27] that can be fixed by requiring that the projection of W ′ onto M (0) αβ and its complex conjugate, vanish at Z = 0, that is The above redefinitions ensure that under the Lorentz transformations with parameters This is the form of the equations typically used to seek exact solutions as it displays the role of the deformed oscillator algebra in the last two equations. Here one may exploit the technology developed in study of noncommutative field theories and construction of projection operators in a suitably defined oscillator space.

The nonminimal chiral model in Kleinian space
In this model the spinor oscillators are now representations of SL(2, R) L × SL(2, R) R , and as such their hermitian conjugates are now given by The field equations are now given by These models are referred to as being chiral in view of the half-flatness condition on the twistor space curvature, namely Fαβ = 0. These models admit the coset space H 3,2 = SO(3, 2)/SO(2, 2) as a vacuum solution, which has the Kleinian signature (2,2). For a detailed description of these spaces, including the curved Kleinian geometries, see [34].
Our motivation for highlighting this case is due to its being the first exact solution of Vasiliev equation in which all HS fields are nonvanishing, and that the Kleinian geometry is relevant to N = 2 superstring as well as integrable models.

The method
In order to construct solutions to Vasiliev's equations, one may consider the approach [19] in which they are homotopy contracted in simply connected spacetime regions U to deformed oscillator algebras in twistor space at a base point p ∈ U ; the constraints are thus integrated in U using a gauge function g = g(Z, Y |x) obeying and initial data subject to The fields in U can then be expressed explicitly as after which the Lorentz covariant HS gauge fields can be obtained from (2.22) subject to (2.23), which serves to determine the spin connection ω µab . Thus, the deviations in the spacetime HS gauge fields away from the topological vacuum solution, that is the solution with W µ = 0, thus come from the gauge function g as well as the non-linear shift on the account of achieving manifest Lorentz invariance. The deformed oscillator algebra requires a choice of topology for Z 4 , initial data for B ′ and a flat background connection. In what follows, we shall assume that Z 4 has the topology of R 4 with suitable fall-off conditions at infinity [17,7], and impose for nontrivial flat connections on Z 4 , that are not pure gauge, see [21]. The gauge function represents a gauge transformation that is large in the sense that it affects the asymptotics of gauge fields so as to introduce additional physical degrees of freedom to the system, over and above those contained in the twistor space initial data and flat connection; strictly speaking, in order to define such transformations, one should first introduce a set of classical observables forming a BRST cohomology modulo a set of boundary conditions on ghosts, after which a large gauge transformation is a gauge transformation that does not preserve all the classical observables. In particular, in order to describe asymptotically maximally symmetric, or Weyl flat, solutions, one may take where L = L(Y |x) is a metric vacuum gauge function, to be described below. In order to obtain exact solutions, we shall choose g = L for all C ′ , that we refer to as the L-gauge. However, in order to extract Fronsdal fields in the asymptotic region, one has to impose a gauge condition in twistor space to the leading order in the weak field expansion in the asymptotic region, which introduces a dressing of the vacuum gauge function by an additional perturbatively determined gauge function; see Section 6.

Vacuum solutions
In order to obtain solutions containing locally maximally symmetric asymptotic regions, one may take the gauge function L(Y |x) to be corresponding coset representatives. In what follows, we shall focus on the spaces AdS 4 = SO(3, 2)/SO(3, 1), dS 4 = SO(4, 1)/SO (3,1) and H 3,2 := SO(3, 2)/SO(2, 2), which can be realized as the embeddings These spaces can be conveniently described in a unified fashion using the stereographic coordinates x a ± (a = 0, 1, 2, 3) obtained by means of the parametrization where U ± denotes the two stereographic coordinates charts, each covering one half of the space (3.9); on the overlap one has λ 2 x 2 ± = −1, and the coordinate transition function where the reflection map The boundary is given by λ 2 x 2 ± = 1, which has the topology of S 2 × S 1 in the case of AdS 4 and H 3,2 , and S 3 ∪ S 3 in the case of dS 4 . Instead of covering the vacuum manifold with two charts, one may extend either one of the charts to R 4 \ {x a : λ 2 x 2 = 1}, which provides a global cover using a single chart, with the understanding that {x a : (3.14) On |λ 2 |x 2 < 1, the corresponding vacuum gauge function where A global description can be obtained using two gauge functions , then the two locally defined solutions on U ± can be glued together using the gauge transition function T + − := L −1 + ⋆ L − = 1 defined on the overlap region where λ 2 x 2 ± = −1. For later purposes, it is convenient to introduce alternative coordinate systems which are defined by the embeddings (with |λ| 2 = 1) x 2 > 0 : X 0 = sinh τ sin ψ , X i = n i cosh τ sin ψ , X 5 = cos ψ , The metrics for (A)dS in these coordinate systems are given in (3.41)- (3.44). In the case of H 3,2 , we will find the following coordinate system to be useful (with |λ 2 | = 1) where n i n j δ ij = 1, 0 ≤ r ≤ ∞ and 0 ≤ t ≤ 2π.

Instanton solutions of minimal model in (anti) de Sitter space
Having obtained a vacuum gauge function, the next task is to solve the deformed oscillator problem (3.4) subject to the initial data in twistor space. To this end, it is helpful to constrain the primed configurations further by assuming that they preserve a nontrivial amount of HS symmetries. This can be achieved by imposing symmetry conditions by seeking a subspace k ′ of HS gauge parameters ǫ ′ that are unbroken, i.e.
for all ǫ ′ ∈ k ′ ; upon switching on the gauge function g, the resulting full solution is invariant under gauge parameters in the space k = g −1 ⋆ k ′ ⋆ g. For example, one may require that an n-dimensional subalgebra g n of the maximal finite dimensional subalgebra g 10 of the HS algebra remains unbroken, which implies that k ′ is given by the intersection of Env(g n ) and the HS algebra. In particular, taking n = 10 yields the vacuum solution which preserves the HS algebra itself. Taking n < 10, the first distinct cases with nontrivial Weyl zero-form arise for n = 6; the space k ′ is then given exactly, as we shall describe below for a particular realization of g 6 , or perturbatively. In the latter case the g 6 will be realized as an algebra of the HS algebra in the leading order.
In [20], asymptotically anti-de Sitter solutions with g 6 = o(1, 3) were constructed by taking k to be generated by the full Lorentz generators M αβ from (2.20). Thus, in the primed basis, the corresponding symmetry conditions read , which are given by y α y β plus perturbative corrections, and whose star product commutators close modulo Lorentz transformations acting on the component fields; thus, consistency of the invariance conditions implies that all canonical Lorentz tensors that are not singlets must vanish. Alternatively, it is possible to use other embeddings of o(1, 3) into the algebra of primed HS gauge transformations; for example, one can simply take y α y β , as we shall comment on in Section 5, though it remains an open problem whether the resulting solutions are gauge equivalent to those that will be presented below. Taking g n to be generated by unperturbed functions of Y is useful, however, in considering unbroken symmetry algebra involving transvection operators, as we shall spell out in further detail in Section 5 in the case of domain walls and related time dependent solutions.
Turning to (3.24), the simplest possible ansatz for B ′ is a constant, viz.
which leads to a deformed oscillator problem with exact solutions closely related to those of the 3D HS theory constructed by Prokushkin and Vasiliev [8]. Adapting to the 4D Type A model, for which b = 1, the following solution for the twistor space connection was found in [20] S Expanding exp(itu/2) results in integrals of the degenerate hypergeometric functions times positive algebraic powers of t, which improve the convergence at t = 0. Thus V α is a power series expansion in u with coefficients that are functions of ν that are well-behaved provided this is the case for the coefficient of u 0 . This is the case for ν in some finite region around ν = 0, as discussed in detail in [20,35]. Indeed, as we shall see below after carrying out the integration over t, ν must lie in the interval −3 ≤ ν ≤ 1.
The solution in spacetime is obtained in the two regions λ 2 x 2 < 1 and λ 2 x 2 < 1 using the stereographic gauge functions L ≡ L(x|Y ) and L ≡ L(x|Y ) where x andx are related by the reflection map in the overlap region where λ 2 x 2 < 0 and λ 2 x 2 < 0. From (3.16), one finds [20] This shows that the physical scalar field is given in the x a -coordinate chart by while the Weyl tensors for spin s = 2, 4, . . . vanish. Using instead L, the physical scalar field in the x a -coordinate chart is given by As a result, the two scalar fields are related by a duality transformation in the overlap region Thus, if the transition takes place at λ 2 x 2 = λ 2 x 2 = −1, then the amplitude of the physical scalar never exceeds 2ν.
The master fields S A and W ′ µ are obtained from (3.5) and (2.22). The generating functional for the spacetime gauge fields is given by [20] 1 subject to (2.23), which serve to determine the spin connection, and where a αα is defined in (3.16), and This gives [20] exhibiting branch cuts for Re ν ≤ −3 and Re ν ≥ 1 2 . For ν ≪ 1, and in the interval −1 ≤ a 2 ≤ 1, this function can be approximated by [20] Using (2.23), one determines ω αβ and e αα from (3.32), while the HS Fronsdal potentials φ µa 1 ...as vanish for s > 2: Even though the HS fields vanish, it is to be noted that the solution of the metric and scalar field constitute a solution of a highly nonlinear system of equations in which all 1 An Euclidean version of this solution has been obtained in [21], and as the spin connection plays an eminent role in this solution and assuming that the action, as proposed in [17], is finite on this solution, we use the terminology of "instanton solution". 2 The unitary representations of Wigner's deformed oscillator algebra can obtained starting from the standard Fock space and factoring out ideals that depend on integer part of (1 + ν)/2, that is, the ideal jumps for odd values of ν [36,37,38,39]. It would be interesting to examine to what extent it is possible to extend the solution to general ν properly taking into account the branch points in Q at odd ν.
higher derivatives play a role. One can reverse engineer a two derivative action describing the coupling of gravity to scalar field that admits the same solution [20] but such an action is clearly of limited use in the context of HS theory.
An advantage of presenting the solutions in stereographic coordinates is that it facilitates their unified description for (A)dS. In these coordinates the solution for the scalar field is given by (3.29) and the metric by and it is understood that the integration variable in It is also convenient to give the result in the coordinate system defined in (3.19). In these coordinates, the solution takes the form [35] AdS 4 : where we have set |λ| 2 = 1 and In addition to the SO(3, 1) symmetry generated by M ′ αβ , the solution is also left invariant by additional transformations with rigid HS parameters where the ℓ'th level is given by [20] ǫ The full symmetry algebra is thus a higher-spin extension of SO(3, 1) ≃ SL(2, C), that we shall denote by where sl(2, C) is generated by M ′ αβ and its hermitian conjugate, and we have indicated that in general the structure coefficients may depend on the deformation parameter ν.
Turning back to the solutions given above, a holographic and cosmological interpretation of (3.41) has been discussed in [35], where a bouncing cosmology scenario was observed, and its comparison with a similar phenomenon occurring supergravities [40] was made. In this context, it is useful to examine the behaviour of the solutions, both with AdS and dS asymptotics, near the boundary as well as distant future.
For |ν| ≪ 1 and near the boundary, where λ 2 x 2 → 1, or equivalently a 2 → 1, the scale factors η and Ω behaves as 3 , which means that the solutions are asymptotically maximally symmetric spacetimes with undeformed radius.
Another interesting limit to consider is η → 0; for |ν| ≪ 1 this takes place for a 2 + 1 ≪ 1, that is, for a 2 is close but not equal to −1, which corresponds to τ → ±τ crit in the AdS case, and ψ → ψ crit in the dS case. In the former case, we have 4 In Einstein frame 5 , it takes infinite proper time to reach the critical surface, which means that one may interpret the future region of the solution as a singularity free SO(3, 1) invariant cosmology with a finite asymptotic scalar field, as In the dS case, one may instead interpret the critical limit as a domain wall at an infinite proper space-like distance from the center of the solution.

Solutions of the non-minimal chiral model in Kleinian space
In obtaining the solutions described above, symmetries on the master fields were imposed. In [21] projection operators were used as well. In the case of the non-minimal chiral model 6 in Kleinian space, it is possible to use projectors to build solutions with non-vanishing Weyl zero-form and HS fields. They are where P ⋆ P = P ,P ⋆P =P , and satisfy the conditions π •π(P,P ) = (P,P ). The projectors P andP are independent. Consider the simplest case in which where, in terms of constant spinors (λ α , µ α ) we have defined Upon the L-dressing, and expanding the result in Y -oscillators, the following component results are found in [21] while the anti-self-dual Weyl tensors take the form In stereographic coordinates, the Kleinian space is covered in two charts with 0 ≤ h 2 ≤ 2, and hence the Weyl tensors blow up in the limit h 2 → 0 preventing the solution from approaching H 3,2 in this limit. In coordinate system introduced in (3.21), the metric reads ds 2 = −(dr 2 + r 2 dt 2 ) + (1 + r 2 )dΩ 2 2 , and the pre-factor in this solution reads which, indeed, diverges at the boundary r → ∞ as noted above. In [21], it was also found that where For the spin connection, the result is (3.63) Note that in the last term of the second equation the full spin connection ω γδ arises. The metric g µν = e a µ e b ν η ab takes the form The vierbein has potential singularities at h 2 = 0 and h 2 = 2. The limit h 2 → 0 is a boundary at which e µ a ∼ h −2 x µ x a , i.e. a scale factor times a degenerate vierbein. In the limit h 2 → 2 one approaches the boundary of a coordinate chart. Also in this limit, the vierbein becomes degenerate, viz. e µ a ∼ h −2 (Jx) µ (Jx) a .

Perturbative construction of domain-wall solution
As mentioned earlier, if we wish to construct a solution of the HS equations that has a symmetry group that includes any translation generators P a , given that there is no known realization of these generators that would form a closed algebra with the full Lorentz generators M αβ , the symmetry conditions (3.22) need to be imposed in terms of the undeformed generators that are bilinear in the oscillators. The deformation of this symmetry to accommodate nonlinear corrections can then be computed perturbatively in a weak field expansion scheme. This is the framework which was pursued in considerable detail in [20], where the perturbative construction of solutions with 3, 4 and 6 parameter isometry subgroups of the AdS group were considered. Here we shall outline the key aspects of this constructions by describing the example of a domain-wall solution having ISO(2, 1) symmetry and its appropriate HS extension [20].
Thus, the symmetry conditions to be imposed are As shown in [20], these conditions are solved by Denoting the ISO(2, 1) transformations discussed above by ǫ ′ (0) , we can seek its nonlinear deformation by expanding where ǫ ′ (n) is constant in spacetime but may depend on Y and Z. The symmetry condition at first order is satisfied by C ′ given in (3.69). To establish the symmetry at second order, we need to satisfy where B ′ (2) is obtained from the normal ordered perturbative scheme (see section 6.1 below) to be [20] B ′ This condition (3.71) is solved in [20] where it is found that where λ αβ , λ αβ are arbitrary constant parameters. For a more detailed discussion of the procedure outlined above, see [20].

Other known solutions
The solutions described in Section 3.3 were generalized in [21] to find new Lorentz-invariant vacuum solutions, in which in addition to the continuous parameter ν, an infinite set of independent and discrete discrete parameters θ k = {0, 1}, each turning on a Fockspace projector P n (u), were activated. Should they be proved to be gauge-inequivalent to AdS 4 , as they seem to be, they would represent monodromies of flat but non-trivial connections (V α ,Vα) on Z. An interesting limit of this solution arises upon setting ν = 0 and (θ k − θ k+1 ) 2 = 1, leading to the degenerate metric . (3.74) The methods above can also be extended to the Prokushkin-Vasiliev theory in D = 3, giving rise to Lorentz-invariant instanton solutions (with additional twisted sectors of the theory excited, and the characteristic extra deformation parameter λ that allows to vary the mass of the scalar), as well as to the above projector vacua [31].
Finally, in [41] a different class of exact solutions was constructed by means of the gauge function method coupled with a different choice of gauge in twistor space, there referred to as axial gauge. As for the perturbative construction of solutions, it is worth mentioning the plane wave solution of [19,4], whose elevation to an exact solution remains to be investigated, to our best knowledge.

The method
The method developed in [6,22,7] for finding the exact solution of Vasiliev equations also exploits the gauge function method to solve for (B, V, W ) in terms of (B ′ , V ′ ) from (3.5) and (2.22), which, in turn, are to be determined by solving (3.4). It is here where the method differs from the method described above, by making the following factorized ansatz for (B ′ , V ′ ): Note that this ansatz for V ′ α is holomorphic in z, and in the following we shall refer to the solution in this form as given in holomorphic gauge [7]. In this section we shall consider the nonminimal model in which we recall that conditions (2.35) apply. It is shown in [7] that this ansatz solves the fully non-linear equations (3.4) provided that and that obeys the deformed oscillator algebra Note that only the first power of ν in (4.4) survives in the commutator (4.5). One class of solutions is given by [6] n≥1 where the constant spinors v ± (0)α are used to defined the projected oscillators The presence of z + z − term breaks manifest Lorentz covariance, which will be restored when we consider the field dependent gauge transformation in Section 7 that is needed to cast the results into a form that can be interpreted in terms of Fronsdal fields that obey the standard boundary conditions.
At this point, B ′ and V ′ α are determined, with Ψ(Y ) representing an arbitrary initial datum. One can proceed to compute (B, V, W ′ ) from (3.4), (3.5) and (2.22). However, one needs to ensure that the star products involving Ψ are well defined. The analyticity properties of the resulting (B, V, W ′ ) also require special care. The strategy adopted in [6,22,7] is to employ projection operators with well defined group theoretical origin, and easily deducible symmetry properties. We shall illustrate aspects of this procedure below with a relatively simple example, namely the black hole-like solution [7] closely related to that of [24], which we shall also describe in a subsequent section below.

Black hole solution
In seeking an exact solution of Vasiliev equations which has the symmetries of 4D static black hole solution that is asymptotically AdS 4 , namely spatial rotations and time translations generated by

SO(3) × SO(2) :
M rs (r, s = 1, 2, 3) , and M 05 = E = 1 4 (σ 0 ) αα y αȳα . (4.8) The first instance of such a solution was found in [24] in a different approach that will be summarized later. In both approaches the following projector plays an important role: where P ′ is given by and the reality conditions dictate that with M ∈ R. This projector clearly has the desired symmetry property since [M rs , P ′ ] ⋆ = 0 and [E, P ′ ] ⋆ = 0. In performing the L-dressing, we need the result where K AB are the Killing parameters taking the form In an x-dependent eigenspinors of u αβ , these parameters can be expressed as (4.14) The Kerr-Schild vector k αα associated with E is obtained via a projection of v αα with the eigenspinors of u αβ ,ūαβ, and can be written as 7 In obtaining the above results, the computation is first done in stereographic coordinates, and then coordinate change is made to go over to the spherical (global) coordinate system in which the AdS metric reads ds 2 = −(1 + r 2 )dt 2 + (1 + r 2 ) −1 dr 2 + r 2 dΩ 2 2 . From the point of view of the factorized Ansatz (4.1)-(4.2), the fact that Ψ is proportional to a projector can be seen as a way of enforcing the Kerr-Schild property of a black-hole solution, as it effectively causes a collapse of all non-linear correction in V ′ α ,V ′ α down to the linear order, at least from the point of view of the oscillator dependence. Coupled with the gauge freedom on S α this allows to effectively reach a gauge in which the full solution only contains first order deformations in ν [6]. Another noteworthy fact is that, due to the factorized dependence on Y and Z, one can effectively separately rotate the two oscillators by means of a factorized gauge function As discussed in [6,22,7], turning on the second factor is useful since by choosing it appropriately one can achieve collinearity between the spin-frame (v + , v − )(x) on Z (which is obtained by pointwise rotation of v (0) given in (4.7) by L) and the eigenspinors (u + α , u − β ) of u αβ in order to remove singularities that appear in the solution for the master one-form, that are gauge artifacts. We note that the factor L(x, Z), being purely Z-dependent, does not affect B and only acts non-trivially on V α ,Vα.
As noted earlier, this solution is closely related to that of Didenko and Vasiliev [24]. Indeed, the solution for B is the same, while the relationship between the solutions for (S A , W ′ ) is more subtle and it will be discussed in the next section.
Note that B does not depend on Z. Thus, B| Z=0 = B = C(Y ), and its holomorphic components give the Petrov type -D Weyl tensors 24) and similarly for the anti-holomorphic components givingCα 1 ...α 2s . The singularity of individual Weyl tensors does not necessarily imply a physical singularity in HS gravity for the following reasons. At the level of the master fields, r = 0 also appears as the only point at which V α , as well as W ′ , acquire a pole on a plane in Z × Y defined by a α | r=0 = z α + i(σ 0ȳ ) α = 0, due to the zero at t = −1 that the denominator of the integrand in (4.18) develops at r = 0. The master-field curvature is however given by B ⋆ κ and one can argue that at the master-field level, which is the only sensible way to look at such solution in the strong-field region, B remains, in fact, regular at r = 0. Qualitatively this can be understood as follows. r appears in (4.17) as the parameter of a delta sequence: away from the origin one has a smooth Gaussian function, approaching a Dirac delta function on Y as r goes to zero [6]. However, unlike the delta function on a commutative space, the delta function in noncommutative twistor space, thought of as a symbol for an element of a star product algebra, is smooth. Indeed, it is possible to show [6] that by changing ordering prescription one can map the delta function to a regular element, and the smoothness of such change of basis manifests itself in the fact that the solution of the deformed oscillator problem obtained in the new ordering can be mapped back smoothly to the solution above. In this sense, the singularity in r = 0 may be an artifact of the ordering choice for the infinite-dimensional symmetry algebra governing the Vasiliev system.
In order to extract the x-dependence of even just the spin-2 component of W ′ , one still needs to evaluate the complicated parametric integrals in (4.22a). However, as the Weyltensors take the simple form in x-space given in (4.24), we expect that a suitable gauge transformation exists that will give the metric in the standard Kerr-Schild form, namely g µν = g AdS µν + 2M k µ k ν /r. Finally, let us note that the basic black hole-like solution reviewed above has a generalization in which infinite set of projection operators P n and twisted projection operators P n are introduced 8 . The twisted projectors P n are invariant under SO(3) × SO(2) discussed above while the projectors P n are invariant under SO(3) and the (B ′ , S ′ A ) sector of the 8 Even though P ′ n ⋆ P ′ n = P ′ n , we refer to P ′ as twisted projector to emphasize the fact that it is related to the projector P ′ n by the relation P ′ n = P ′ n ⋆ κy.
solution takes the form [6,22,7] B ′ = n=±1,±2,... ν n P ′ n + ν n P ′ n , P ′ n = P ′ n ⋆ κ y , (4.25) where the contour integrals are performed around a small contour C(ε) encircling ε. The expressions for (V n,A , V n,A ), and their g = L ⋆ L dressing can be found in [7]. It turns out that the solutions with only ν n parameters switched on correspond to black hole like solutions, of which the case summarized above arises for n = 1. Solutions with only nonvanishing ν n parameters correspond to massless particle modes, and surprisingly black hole modes as well entering from second order onwards in a perturbative treatment of the solution; for a detailed description of this phenomenon see [7,26]. Note that the SO(3) invariant projectors are associated with spin-0 modes. To extend this construction to spin-s particle modes, one needs the spin-s generalization of the projectors P ′ n discussed above. A particular presentation of such projectors can be found in [42].

Other known solutions
By means of the same factorized Ansatz (4.1)-(4.2) black-hole-like solutions with biaxial symmetry have also been found [6,22,43], some of which being candidate HS generalizations of the Kerr black hole [43]. The separation of variables in holomorphic gauge was also instrumental to finding solutions with g 6 isometries of cosmological interest, that we shall present in a forthcoming paper [23].

Direct Method and the Didenko-Vasiliev Solution
While, as we have seen so far, the gauge function method is in general of great help in constructing exact solutions, it is sometimes possible to attack the equations directly, by virtue of some other simplifying Ansatz or gauge condition. We shall generically refer here to any method which does not rely on the use of gauge function as direct method. One such solution has been found so far in this way by Didenko and Vasiliev [24], which has nonvanishing HS fields, and which contains the Schwarzschild black hole solution in the spin 2 sector. Indeed, motivated by the phenomenon that solutions of Einstein equations that can be put in Kerr-Schild form solve the linearized as well as the nonlinear form of the equations, the Authors of [24] thought of an Ansatz that would generalize some distinctive features of black hole solutions in gravity. First, it is based on an AdS timelike Killing vector, in the sense that the Weyl tensor will be a function of some element K AB as (4.13). More precisely, if f (K) satisfies the Killing vector equation, a proper ansatz for a solution of the linearized twisted adjoint equation will be given by f (K) ⋆ κ y . Second, they chose f (K) in such a way that the Ansatz linearize the Vasiliev equations. For the latter purpose it is important that the function f (K) is a projector -in fact, a Fock-space vacuum projector that coincides with (4.12). Such choice in particular reduces equations (2.31) and (2.32) to two copies of the 3D (anti)holomorphic deformed oscillator problem that arises in Prokushkin-Vasiliev HS theory in 3D [8] in terms of the oscillators in (4.21) 9 .
The ansatz [24] where g is another function to be determined is then employed to deal with the remaining equations that involve W , namely (2.28), (2.29) and (2.30). The resulting exact solution is given by [24] B = 4M r exp 1 r 2 1 2 y α u αβ y β + 1 2ȳαūαβȳβ + y α u αβ v βαȳα , (5.3a) (5.4) 9 In this section we shall use the conventions of [24] which differ form ours.
Note also that S A does not satisfy the Vasiliev gauge S(Z, Y )| Z=0 = 0, and that W above has not been redefined as in (2.22). Nonetheless, it has been noted in [24] that a HS transformation of the form δW = D 0 ǫ (1) with ǫ (1) = − 1 2 1 0 dtz α S α | z→tz + h.c. + f (Y |x) and arbitrary f (Y |x) maps W to W phys given by whose spin 2 component gives the frame field and the associated metric where k µ = e αα 0 k αα . This is the metric of a black hole of mass M in AdS 4 in Kerr-Schild form. The terminology of black hole in HS context requires caution as discussed in the introduction. In addition to the SO(3) × SO(2), the solution summarized above has been shown to also have 1/4 of the N = 2 supersymmetric HS symmetry of the model, and their infinite dimensional extension thereof [24].
The solution (5.3) differs from (4.17)-(4.19) both in the form of the internal connection (V ± ,V ± ) and in that of the gauge field generating functions. As for the internal connection, the difference can be ascribed to the two choices employed in [6] and [24] for solving the deformed oscillator problem (referred to as "symmetric" and "most asymmetric", respectively, in [6,22]). One can show that the resulting internal connections can be connected via a gauge transformation (see [6]), although the small or large nature of this transformation is yet to be investigated. The comparison of W ′ in (4.19) and W in (5.3) is technically more complicated, as the two differ also by the shift of the Lorentz connection (2.22), and it will be postponed to a future work, but we note that having the same B identical in both solutions strongly suggests that the physical gauge fields should be equivalent, and in particular equivalent to (5.5).
It is worth mentioning that even by working without the gauge function method, with a specific choice of gauge the Didenko-Vasiliev solution can be simplified in such a way that the W connection is reduced to the vacuum one W 0 . This simplification was studied in [44], along with the embedding of the solution in the N = 2 and N = 4 supersymmetric extensions of the bosonic Vasiliev equations.
In the normal order, defined by the star product formula (2.4), the inner Klein operators become real analytic in Y and Z space, viz. κ = κ y ⋆ κ z = exp(iy α z α ) ,κ = κȳ ⋆ κz = exp(iȳαzα) . (6.1) Assuming that the full field configurations are real-analytic on Z 4 for generic points in X 4 , one may thus impose initial conditions where the x-dependence is understood. In order to compute the Z-dependence of the fields one may choose V A | C=0 , which is trivial flat connection, to vanish, and a homotopy contractor for the de Rham differential on Z 4 , which entails imposing a gauge condition on V A . One may then solve the constraints on D A B, F AB and F Aµ on Z 4 in a perturbative expansion in C. This procedure gets increasingly complex with increasing order in the expansion, which schematically can be written as where B (n) , V (n) and W (n) (a µ ) are n-linear functionals in C, and W (n) (a µ ) is a linear functional in a µ . These quantities, which are constructed using the homotopy contractor on Z 4 , depend on Z, and are real-analytic in Y 4 × Z 4 provided that C and a µ are real analytic in Y -space and all star products arising along the perturbative expansion are well-defined. As for the remaining equations, that is, that F µν = 0 and D µ B = 0, it follows from the Bianchi identities that they are perturbatively equivalent to F µν | Z=0 = 0 and D µ B| Z=0 = 0, which form a perturbatively defined Cartan integrable system on X 4 for C and a µ .
To Lorentz covariantize, one imposes with W ′ from (2.22), that is where w does not contain any component field proportional to y α y β andȳαȳβ in view of (2.23). Upon substituting the above relation into W (n) (a µ ; C, . . . , C), it follows from the manifest Lorentz covariance that the dependence of F µν | Z=0 and D µ B| Z=0 on the Lorentz connection arises only via the Lorentz covariant derivative ∇ and the Riemann two-form (r αβ , rαβ). Thus, the resulting equations in spacetime take the form [6] ∇w + w ⋆ w + 1 4i r αβ y α y β + h.c. + i αβ , (6.10) Alternatively, in order to stress the perturbation expansion around the maximally symmetric background, including the spin-two fluctuations, it is more convenient to work in terms of the original one-form a µ . The perturbative expansion up to 3rd order in Weyl curvatures reads where (V, V 2 , U ) are functionals that can be determined from (6.8) and (6.9). One next assumes that the homotopy contraction in Z-space is performed such that which we refer to as the Vasiliev gauge, and expands a = W 0 + a 1 + a 2 + · · · , C = C 1 + C 2 + · · · . (6.15) where W 0 is the maximally symmetric background; a 1 a and C 1 are linearized fields; a n and C n are nth order fluctuations. The resulting linearized field equations on X-space provides an unfolded description of a dynamical scalar field and a tower of spin-s Fronsdal fields where we use the convention that repeated indices are symmetrized. Computing the functional V(W 0 , W 0 , C), the linearized unfolded system is given by [4] where The oscillator expansion of (6.18) furnishes the definition of the spin−s Weyl tensors, and gives the field equations for spin-s fields which remarkably do not contain higher than second order derivatives, and indeed they are the well known Fronsdal equations for massless spin s-fields in AdS. As for (6.19), its oscillator expansion gives the AdS massless scalar field equation in unfolded form.
Perturbative expansion around AdS at second order is rather complicated but still manageable. Schematically they take the form Even though these terms have been known in the form of parametric integrals for some time, their detailed structure and consequences for the three point functions were considered much later in [11,48], where the field equation for the scalar field was examined, and used for computing the three-point amplitude for spins 0 − s 1 − s 2 . If s 1 = s 2 , only the first source term in (6.22b) contributes and gives a finite result, in agreement with the boundary CFT prediction. However, if s 1 = s 2 , only the second source term in (6.22b) contributes and gives a divergent result [11,48]. This divergence was confirmed later [12,49] in the three point amplitude for spins s − 0 − 0, resulting from the last term in (6.22a). Soon after, it was shown that a suitable redefinition of the master zero-form cures this problem [50,47], as has been also confirmed with the computation of relevant three-point amplitudes [51,52]. A similar redefinition in the one-form sector has also been determined so that the divergence problem arising in the last term in the first equation above has also been removed [53].
In determining the higher order terms in the perturbative expansions of Vasiliev equations, it remains to be established in general what field redefinitions are allowed in choosing the appropriate basis for the description of the physical fields. In wrestling with this problem, the remarkable simplicity of the holographic duals of this highly nonlinear and seemingly very complicated interactions may provide a handle by means of their holographic reconstruction. Such reconstruction has been achieved for the three and certain four-point interactions [13,14,15]. Putting aside the analysis and interpretation of the nonlocalities [14,54,55], which are present, and nonetheless in accordance with holography by construction, the issue of how to extract helpful hints from them with regard to the nature of the allowed field redefinitions in perturbative analysis of Vasiliev equation remains to be seen.
Of course, ultimately it would be desirable to have a direct formulation of the principles that govern nonlocal interactions, based on the combined boundary conditions in twistor space as well as spacetime, as we shall comment further below.

A Proposal for an Alternative Perturbation Scheme
In what follows, we shall show that for physically relevant initial data Ψ(Y ) given by particle and black hole-like states, the solutions obtained using the factorization method can be mapped to the Vasiliev gauge used in the normal ordered perturbation scheme at the linearized level. Whether the Vasiliev gauge is compatible with an asymptotic description in terms of Fronsdal fields to all orders in perturbation theory, or if it has to be modified, possibly together with a redefinition of the zero-form initial data, is an open problem. In finding the proper boundary conditions in both spacetime and twistor space it may turn out to be necessary to require finiteness of a set of classical observables involving integration over these spaces.
One can define formally the aforementioned map to all orders of classical perturbation theory by applying a gauge function to the holomorphic gauge solution space, where H = 1 + n 1 H (n) is a field dependent gauge function to be fixed as to impose the Vasiliev gauge in normal ordering. To begin with, let us consider fluctuations around AdS for which Ψ(Y ) consists of particle states, focusing, for concreteness, to the case of scalar particle states worked out in [7]. Upon switching on the gauge function L, the field V (L) α := L −1 ⋆ V ′ α ⋆ L develops non-analyticities in the form of poles in Y -space in the particle sector [7]. Applying H removes these poles and expresses the results in terms of Fronsdal fields (at the same time restoring the manifest Lorentz covariance, as we shall see), at least in the leading order. This can be see as follows. We want to obtain V (G) α such that The leading order gauge transformation reads Contracting by z α and using the fact that by definition Z A V (G) A = 0, one finds In particular, activating only the scalar ground state (and its negative-energy counterpart) via the parameters ν ±1 in the projector ansatz for B (see (4.25)), this was computed in [7] with the result where we have set b = ν 1 = ν −1 = 1, and u := y α z α = y + z − − y − z + , y α = y α + M αβ (x, η)ȳβ .
The matrix M αβ (x, η) can be found in Section 5.2.1 of [7]. This result for H (1) is regular in Z but has a pole in Y. It follows that is now real-analytic everywhere on C. Furthermore, it was shown in [7] that the above expressions for (B, V α ) lead to the relation (7.9) in agreement with the result obtained in the standard perturbative analysis of Vasiliev equations at leading order. The emergence of z α in (7.7) shows that manifest Lorentz covariance is restored, in comparison with the expression for V (L)(1) α . The prefactor in (7.7) is a consequence of the fact that we are considering the lowest mode alone in the solution for Fronsdal equation the scalar field, as opposed to summing all the full set of modes.
Expressing the exact solutions obtained in holomorphic gauge in terms of Fronsdal fields amounts to setting up an alternative perturbation scheme in which one constructs the higher orders of the gauge function H (n) subject to dual boundary conditions, that is, to conditions restricting the both the twistor-space dependence of the master fields and their spacetime asymptotic behaviour. Indeed, one requires that after having switched on H, the master fields have symbols in normal order that are real-analytic at Y = Z = 0, and symbols in Weyl order that belong to a (associative) star-product subalgebra with welldefined classical observables defined by traces over twistor space and integrations over cycles in spacetime. The proposal is that this problem admits a non-trivial solutions, and that it fixes H (n) up to residual small HS gauge transformations and the initial data for the master zero-form B to all orders in classical perturbation theory. It would be interesting to see whether these type of field redefinitions are related to those recently proposed by Vasiliev in order to obtain a quasi-local perturbation theory in terms of Fronsdal fields [50,53]. An important related issue is whether the gauge function G is large in the sense that it affects the values of the HS zero-form charges, which are special types of classical observables given by traces over twistor space defining zero-forms in spacetime that are de Rham closed [20,27,56].
A nontrivial test of the factorization approach is first to show that the solution is finite after performing the higher order H (n) gauge transformations, and second, to show that the resulting n-point correlators are in agreement with the result expected from holography. The corrections beyond the leading order remain to be determined, while the computation of correlators has been performed in which the second order solution in standard perturbative scheme has been used. It has been shown that a naive computation of B (2) in the standard perturbative scheme leads to divergences [11,12], and later it was shown that these divergences can be removed by a suitable redefinition of C(Y |x) [50,47]. Whether there exists a principle based on any notion of quasi-locality in spacetime that governs the nature of such redefinitions to all orders in perturbation is not known, to our best knowledge.
An advantage of the factorization method is that here we start from a full solution to Vasiliev's theory, defined as a classical field theory on the product of spacetime and twistor space (not referring a priori to the conventional perturbative approach). This provides a convenient framework for the description of the solutions with particles fluctuating around nontrivial backgrounds. The key principle here is that linear superposition principle holds for the zero-form initial data Ψ(Y ). For example, if we want to describe the solution to Vasiliev equations for particles propagating around BH solution of Section 4, one simply takes Ψ = ν n P n + ν n P n . The exact solution for the combined system is obtained in this way, but in small fluctuations of a particle propagating in a fixed and exact black hole background, one may treat the parameter ν n and ν n as small and large, respectively. A very interesting open problem is thus to combine this scheme with the aforementioned proposal for dual boundary conditions in order to work out new types of generating functions for HS amplitudes.

Aspects of Higher Spin Geometry
HS theory has been mostly studied at the level of the field equations and in terms of locally defined quantities in ordinary spacetimes or twistor space extensions thereof. It is clearly desirable to develop a globally defined framework for the classical theory, in order to provide geometric interpretations of the exact solutions, and shed further light on the important issue of the choice of boundary conditions on the master fields on the total space required for Vasiliev's equations to produce physically meaningful anti-holographic duals of boundary conformal field theories. An attempt at a global description of HS theory was made in [27] where bundle structures based on different choices of structure groups, soldering forms and classical observables were considered 10 . Here, we shall highlight a particularly interesting choice of structure group, and the resulting infinite dimensional coset description involving tensorial coordinates, and associated generalized frame field and resulting metrics.

Structure group
Noting that Vasiliev equations are Cartan integrable in arbitrary number of commuting dimensions, yet without changing the local degrees of freedom associated to the zeroforms (more on this below), we consider the formulation of Vasiliev's theory in terms of the master fields (W, B, S α ,Sα) thought of as horizontal forms on a noncommutative fibered space C with eight-dimensional fibers given by Y 4 ×Z 4 and base given by an infinitedimensional commuting real manifold M, consisting of charts coordinatized by X M . In each chart, the one-form field W = dX M W M thus takes its values in the HS Lie algebra hs(4) := P (Y, Z) : τ ( P ) = ( P ) † = − P , (8.1) and the deformed oscillators (S α ,Sα) and the deformation field B are thought of as zeroforms on M, valued in representations of the HS Lie algebra as described in section 2, where the τ -map is also defined. The extended equations of motion are given by Eqs.
(2.28)-(2.32), with the only difference being that the manifold χ 4 is replaced by M.
In order to define the theory globally on M, we glue together the locally defined master fields using transition functions from a structure group, which by its definition is generated by a structure algebra h given by a subalgebra of hs(4). One interesting choice is [27] h = hs + (4) ⊕ sl(2, C) , hs + (4) := 1 2 (1 + π) hs(4) , where sl(2, C) is the algebra of canonical Lorentz transformations. The corresponding connection, also referred to as the generalized Lorentz connection, is given by where ω is the canonical Lorentz connection, while the projection is assumed to form a section 11 . The Vasiliev equations and gauge transformations in terms of these master fields are spelled out in [27].

Soldering mechanism
The horizontal differential algebra on C, which is quasi-free in the sense that the curvature constraints are cartan integrable modulo zero-form constraints, can be projected to a free horizontal differential algebra on the reduced total space C| Z=0 . To this end, one first solves the constraints on Z given initial data where hs + (4) := 1 2 (1 + π)hs(4), we assume that M is soldered by E, that is, the tangent space of M is assumed to be identified with the coset hs(4) / hs + (4) via E. Thus, expanding where P A is a basis for the π-odd elements of hs(4), and denoting the inverse of the frame field E M A by E M A , the local translations with gauge parameters ξ = ξ A P A can be identified by usual means [27] as infinitesimal diffeomorphisms generated by globally defined vector fields ξ = ξ A E M A ∂ M combined with local generalized Lorentz transformations with parameters (ı ξ Γ, ı ξ ω).

Generalized metrics
Taking traces of ⋆-products of generalized vielbeins E = dX M E M and adjoint operators on twistor space, one can construct structure group invariants that are tensor fields on M [27]; in particular, we may consider symmetric rank-r tensor fields on Σ ′ equipped with local coordinates σ m . The area functional is structure group invariant, and its minimum, if well-defined, is Diff(M) invariant, hence serving as a classical observable for the HS theory. Clearly, the dressings by the vertex operator-like operators result in a large number of inequivalent metrics for each r, but this is not a novelty in HS theory, as it is possible to consider similar dressings of the Einstein frame metric in ordinary gravity. The tensorial calculus pertinent to examining the variational principles for r = 2 is well understood, whereas for r > 2 it remains to be investigated further.
In particular, one may ask whether there exists any principle for singling out a preferred metric (possibly of rank two), using, for example, calibrations based on abelian p-forms of the type that will be discussed below. The application of these ideas to a geometrical characterization of the exact solutions of HS theory remains to be investigated.

Abelian p-form charges
Another type of intrinsically defined classical observables facilitated by the introduction of soldering one-forms are the charges of on shell closed abelian p-forms, viz.

On-shell actions
Actions have been proposed that imply the Vasiliev equations [16,17] upon applying the variational principle. Their on-shell evaluations involve subtleties stemming from their global formulation and the crucial role played by the boundary conditions. Thus, it remains unclear whether these actions vanish on-shell. Nonetheless, one can still employ the abelian p-form charges discussed above as off-shell topological deformations. To this end, one set of candidates that have been considered are [27] S Tr κ ⋆ τ 4 R ⋆ R +τ 4 (E ⋆ E + r (S)+ ) ⋆ R + 1 2 (E ⋆ E + r (S)+ ) ⋆2 , which is the curvature of Ω. Using the field equation [27] R + E ⋆ E + r (S)+ = 0 , H [4] , (8.27) where H [4] given in (8.22). Infinities arise from the integration over Σ 4 as well as twistor space. Assuming that the divergence from the AdS vacuum has a definite reality property, such that it can be removed by choosing τ 4 appropriately, one is left with an integral over perturbations that may in principle be finite modulo a prescription for integration contours; for related discussions, see [48,57]. Provided that one considers perturbations corresponding to boundary sources, it would be natural to interpret S top [Σ 4 ] as the generating functional for the boundary correlation functions.
One may also construct topological two-forms given on-shell by with H [2] given in (8.22). One application of this surface operator is to wrap Σ 2 around a point-like defect or singularity such as the center of the rotationally symmetric and static solution of [24]. In this case, the leading order contribution, which comes from the AdS vacuum, is a divergent integral over twistor space. If the divergence has a definite reality property, one can cancel it by choosing τ 2 appropriately. One would then be left with an integral over the perturbations. As the latter involve nontrivial functions in twistor space, the integral may be finite. It would be interesting to seek an interpretation of the resulting value of S top [Σ 2 ] as some form of entropy of the black hole solutions reviewed above.
An alternative framework in which certain 2-and 4-forms in x-space, referred to as the Lagrangian forms, are introduced was proposed in [58], where their possible application for the computation of HS invariant charges and generating functional for the boundary correlations functions were discussed; see also [59] for the computation of such HS charges on black-hole solutions at the first order in the deformation parameters. Asymptotic charges in HS theories have also been described in [60,61,62]. Their relation to the HS charges 12 as well as their evaluation on certain exact solutions will be discussed elsewhere [26].