New classes of bi-axially symmetric solutions to four-dimensional Vasiliev higher spin gravity

We present new infinite-dimensional spaces of bi-axially symmetric asymptotically anti-de Sitter solutions to four-dimensional Vasiliev higher spin gravity, obtained by modifications of the Ansatz used in arXiv:1107.1217, which gave rise to a Type-D solution space. The current Ansatz is based on internal semigroup algebras (without identity) generated by exponentials formed out of the bi-axial symmetry generators. After having switched on the vacuum gauge function, the resulting generalized Weyl tensor is given by a sum of generalized Petrov type-D tensors that are Kerr-like or 2-brane-like in the asymptotic AdS4 region, and the twistor space connection is smooth in twistor space over finite regions of spacetime. We provide evidence for that the linearized twistor space connection can be brought to Vasiliev gauge.


Introduction
Vasiliev's equations [1] (for a recent review, see [2]) provide a fully nonlinear description of higher spin gauge fields in four dimensions coupled to gravity and matter fields. The basic feature of Vasiliev's theory is that the full field configurations are captured by master fields that live on an extension of spacetime by a noncommutative twistor space. The equations admit an exact solution given by the direct product of anti-de Sitter spacetime and an undeformed twistor space. In a specific gauge, certain linearized perturbations of the noncommutative twistor space structure give rise to Fronsdal fields. This suggests a holographic relationship to threedimensional conformal field theories [3][4][5]; see also [6][7][8]. In [9,10] this relation was examined under the assumption that the Gubser-Klebanov-Polyakov-Witten (GKPW) prescription [11,12] for on-shell computations of Witten diagrams can be applied to classical field configurations obtained from Vasiliev's equations.
However, the Fronsdal fields embedded into Vasiliev's master fields have non-local interactions [13][14][15] 1 that belong to a functional class widely separated [17] from that of the quasi-local Fronsdal theory [18], which is built by applying the canonical Noether approach to Fronsdal fields in anti-de Sitter spacetime. The GKPW prescription applies to the quasi-local theory by construction, as its action has self-adjoint kinetic terms, and the resulting holographic correlation functions indeed correspond to free three-dimensional conformal field theories. 2 Recent work [22] shows that there exists an explicit field redefinition that maps Vasiliev's theory to a quasi-local theory on-shell, obtained by carefully fine-tuning the perturbative expansion on the Vasiliev side, though it remains to be seen whether it coincides with that of [18]. Moreover, as later shown in [23] the required field redefinition is large, and hence it is unclear to what extent the method can be used to actually compute any holographic correlation functions. Thus, to our best understanding, the issue of whether holographic amplitudes can be extracted by applying the GKPW prescription to the Fronsdal fields embedded into Vasiliev's master fields remains an open problem.
An alternative approach, pursued in [24], is to seek a weaker relation between the two theories, namely at the level of two distinct effective actions, derived in their own rights following different principles, and then evaluated subject to suitable dual boundary conditions. To this end, one starts from Hamilton's principle applied to a covariant Hamiltonian action formulated using Weyl order on a noncommutative manifold whose boundary is given by the direct product of spacetime and twistor space [24,25]; the Weyl order is required for the noncommutative 1 For a review, see [16]. 2 The functional class encountered in the quasi-local Fronsdal theory in [18] (within the AdS/CFT context) has not yet been identified completely; for a discussion, see [17], [19] and Section 7 of [20]. At the cubic order, the separation between the functional classes of this theory and the Vasiliev theory has been spelled out in [21]. version of the Stokes' theorem to hold and for the imposition of boundary conditions. The resulting variational principle yields Vasiliev's equations in Weyl order, that can be mapped back to Vasiliev's normal order for special classes of initial data in twistor space following the perturbative scheme set up in [26,27]. The resulting form of the higher spin amplitudes [28][29][30] is closely related to first-quantized topological open string amplitudes [31], but nonetheless reproduce exactly the same correlation functions as the Witten diagrams computed in the quasi-local theory. We would like to stress the fact that the Hamiltonian form of the action implies that the dependence of the classical Vasiliev master fields on classical sources are of a different type than for fields obeying equations of motion following from an action with selfadjoint kinetic terms. Indeed, instead of applying the GKPW prescription, the higher spin amplitudes are obtained from functionals given by topological boundary terms added to the Hamiltonian action [24,29,32], whose on-shell values are given by higher spin invariants, as we shall comment on further below.
In this paper, we shall construct new perturbatively defined solution spaces to Vasiliev's equations in Weyl order, by taking into account classes of functions that resemble closely those used in [26]. We shall then demonstrate explicitly that they can be mapped to Vasiliev's normal order, at least at the linearized level, thus providing further evidence in favour of the covariant Hamiltonian approach outlined above.
To this end, we recall that at the linearized level, the fluctuations in the master fields that are asymptotic to anti-de Sitter spacetime form various representation spaces of the anti-de Sitter isometry algebra, including lowest-weight spaces as well as spaces associated to linearized solitons [33] and generalized Petrov type-D solutions [26,33,34]. Nonlinear completions of various Type-D solution spaces were constructed in [26,33,34]; for a review, see [35]. Of direct relevance for the work in this paper is the subspace that contains the the black-holelike solutions, 3 including spherically symmetric solutions. In these solutions, each individual Fronsdal field has a point-like source at the origin, showing up as a divergence in its Weyl tensor. However, upon packing all curvatures into a master zero-form, one obtains the symbol of a quantum-mechanical operator that approaches a delta function distribution at the origin [26], which defines a smooth state as seen via classical observables given by zero-form charges [29,30,37]. In this sense, the black-hole-like Type-D solutions to Vasiliev's theory are source free at the origin. 4 Furthermore, it is possible to dress these solutions with lowest-weight space modes [27] at the fully nonlinear level; in doing so, the latter modes induce Type-D modes 3 This subspace is related to the massless spectrum by means of a Z2-operation [26], reminiscent of a U-duality transformation [36]. 4 It remains to be examined whether additional topological two-forms describing Dirac strings need to be activated in the dynamical two-form [24,38].
already at the second order of classical perturbation theory. 5 Clearly, the full extent of the moduli space of the theory yet remains to be determined.
In this paper, we shall present a new infinite-dimensional class of bi-axially symmetric exact solutions that are asymptotic to anti-de Sitter spacetime and singularity free at the level of zeroform charges. We shall furthermore propose a super-selection mechanism based on requiring that the solutions can be brought to Vasiliev gauge (where the asymptotic linearized fluctuations are in terms of Fronsdal fields).
Our construction method follows closely the one devised in [26] using gauge functions and separation of twistor space variables, which is in effect equivalent to starting from an Ansatz in Weyl order. The key difference is that we shall expand the master fields over a new set of elements in the associative fiber algebra, thus adding a branch to the existing moduli space.
In a generic gauge, the expansion coefficients are functions on the base manifold. However, in the holomorphic gauge of [26] the Weyl zero-form is a constant while the twistor space one-form is given by a universal set of functions, related to Wigner's deformed oscillators, originally derived within the context of three-dimensional matter coupled higher spin gravity [39]. The resulting solution space is then mapped to Vasiliev  Provided that the resulting twistor space connection is smooth at the origin of the base of the twistor space, Vasiliev gauge can be reached by means of a second perturbatively defined gauge transformation. As we shall see, the real-analyticity requirement constrains the initial data in the Weyl zero-form already at the linearized level. 7 More specifically, the new sector of the fiber algebra is isomorphic to the group algebra where Z × Z is generated by two elements in Sp(4; C) given by exponentials of a pair of Cartan generators of sp(4; R). These correspond to linear symmetries of the two-dimensional harmonic oscillator, and generate the Killing symmetries of the solutions (including higher spin symmetries). As we shall see, the aforementioned super-selection rule amounts to restricting the master fields to a subalgebra of the group algebra not containing the unity.
The paper is organized as follows: In Section 2 we review parts of Vasiliev's bosonic higher spin gravity model that we shall use in constructing and interpreting the exact solutions. Solu- 5 This phenomena resembles some of the scattering processes in U-duality covariant field theory [36]. 6 Whether a more general vacuum gauge function can introduce additional classical moduli remains an open problem. 7 An optional criterion is that the fiber algebra is a unitarizable representation of the higher spin algebra and hence the anti-de Sitter isometry algebra; we expect this property to arise at higher orders of classical perturbation theory by requiring positivity of a suitable free energy functional.
tion spaces based on (semi)group algebras are constructed in Section 3 using the aforementioned method; the singular nature of the contribution from the identity is pointed out in Section 3.4.
In Section 4, we show that the Weyl tensor is given by a sum of Petrov type-D tensors that are Kerr-like or 2-brane-like in the asymptotic AdS 4 region, and we compute higher spin curvature invariants. In Section 5, we show in special cases that the twistor space one-form is real-analytic in twistor space over finite regions of spacetime, and that its linearized part can be brought to Vasiliev gauge. We conclude in Section 6.

Bosonic Vasiliev model
In this section, we describe the non-minimal bosonic higher spin gravity model of  [13,26] and the review [2].
We first provide the formal definition in terms of master fields on the direct product of a commuting space and a noncommutative twistor space. We then spell out the component form of the equations, including their reformulation in terms of deformed oscillators. Finally, we remark on choices of bases for the internal algebra, and the Lorentz covariant weak field expansion scheme leading to Fronsdal fields, stressing the role of Vasiliev gauge and smoothness in twistor space.

Master field equations
Vasiliev's original formulation of higher spin gravity is given in terms of two master fields Φ and A of degrees 0 and 1, respectively, and two closed and twisted-central elements I and I of degree 2, all of which are elements of a differential graded associative algebra Ω(M) of forms on a non-commutative manifold M, valued in an internal associative algebra A. Letting ⋆ denote the associative product of Ω(M) ⊗ A, which is assumed to be compatible with d, the fully nonlinear master field equations read where and π denotes an automorphism of the differential graded associative algebra. The two-forms are characterised by the subsidiary constraints for any f ∈ Ω(M) ⊗ A, idem I. Finally, the star functions where b n ,b n ∈ C. It follows that Φ ⋆ π(Φ) and hence B is covariantly constant, viz. In the context of higher spin gravity, it is usually assumed that where X 4 is a four-dimensional real commuting manifold, with coordinates x µ , and Z 4 is a four-dimensional real non-commutative symplectic manifold, with canonical coordinates Z α .
The compatibility between the star product and the differential amounts to the Leibniz' rule The differential star product algebra is assumed to be trivial in strictly positive degrees, in the sense that dΞ M := (dx µ , dz α , dzα) are taken to be graded anti-commuting elements obeying which are consistent with associativity. The algebra Ω(M) ⊗ A is also assumed to be equipped with an anti-linear anti-automorphism †, for which we use the convention In case of the basic bosonic models, without internal Yang-Mills symmetries, the internal algebra A consists of classes of functions on yet one more four-dimensional real non-commutative symplectic manifold, that we shall denote by Y 4 , with canonical coordinates Y α . We shall refer to Y 4 × Z 4 as the full twistor space, and Y 4 and Z 4 , respectively, as the internal and external twistor spaces. 9 The Sp(4;R) quartets are split into SL(2;C) doublets, viz. 10 Y α = (y α ,ȳα) , Z α = (z α ,zα) , The automorphism π and its hermitian conjugateπ are defined by π x µ ; y α ,ȳα; z α ,zα = x µ ; −y α ,ȳα; −z α ,zα , (2.13) π x µ ; y α ,ȳα; z α ,zα = x µ ; y α , −ȳα; z α , −zα , (2.14) and The equations given so far provide a formal definition of the basic bosonic model.

Star product, twisted central element and traces
In what follows, we shall use Vasiliev's original realization of the ⋆-product given by (2.18) We shall encounter ⋆-product compositions leading to Gaussian integrals involving indefinite bilinear forms. To define these we use the fact that the auxiliary integration is a formal representation of the original Moyal-like contraction formula, which means that the integration must be performed by means of analytical continuations of the eigenvalues of the bilinear forms. 9 Taking the master fields to be smooth functions of Y4 yields an anti-de Sitter analog of the Penrose-Newman transformation; to our best understanding, the precise relation between Y4 × Z4 and the original (commuting) twistor space of Penrose remains to be spelled out in detail. 10 The doublet indices are raised and lowered using f α = ε αβ f β and f β = f α ε αβ idem fα.
Symbol calculus. The star product rule implies that [f 1 (y,ȳ) , f 2 (z,z)] ⋆ = 0 , (2.19) that is, the variables Y α and Z α are mutually commuting. Moreover, from it follows that (2.23) One also has 24) resulting in that and also Twisted central element. The condition (2.4) can be solved by where κ y is an inner Klein operator obeying

29)
idem κ z . Thus, one may write where thus By hermitian conjugation one obtains The two-forms j z and jz can be extended to globally defined forms on a non-commutative space Z 4 having the topology of a direct product of two complexified two-spheres [24,26], with nontrivial flux In this topology, it is furthermore assumed that Φ belongs to a section that is bounded at infinity, while the twistor-space one-form is a connection whose curvature two-form falls off at infinity.
We note that the form of I given in Eq. Trf := where j y is given by replacing z α by y α in j z defined in Eq. (2.28). An alternative trace operation, of relevance to higher spin gauge theory, can be defined if A admits the decomposition where A n,n consist of operators whose symbols in Weyl order are regular at the origin of Y 4 .
One may then define the trace operation Trf n,n ⋆ (κ y ) n ⋆ (κȳ)n (2.39) that is, Indeed, in several applications it turns out that Trf is ill-defined while Tr ′ f is well-defined, as for example in the case that f is a polynomial on Y 4 .

Equations in components and deformed oscillators
We decompose the master one-form into locally defined components as follows: The reality condition (2.15) and the bosonic projection (2.17) imply Decomposing master equations into components using inner derivatives ı ∂µ , ı ∂α and ı ∂α , where which are related by hermitian conjugation, and for which the reality condition and the bosonic projection take the form: In terms of the new fields, the aforementioned equations read as can be seen using As we shall see below, the deformed oscillators are useful in defining the field redefinition to Lorentz covariant basis. They also provide a useful basis for finding exact solutions as they convert the differential equations on Z 4 into algebraic equations that can be solved using Laplace transformation methods [39]; for related details, see [26].

Lorentz covariance, Fronsdal fields and Weyl tensors
To arrive at a perturbative formulation in terms of Fronsdal fields on X 4 , one first solves Eqs.
(2.47)-(2.50) subject to an initial datum for Φ and U µ at Z α = 0 in a perturbative expansion in the zero-form initial data in Vasiliev gauge 11 In this gauge, initial data for the zero-form given by generic smooth symbols on Y 4 yields twistor space configurations that are smooth functions on Y 4 × Z 4 . Letting ω αβ µ denote the canonical Lorentz connection, one can show that [42] Φ, V α and 12 11 At the linearized level, this gauge yields the canonical basis for unfolded linearized Fronsdal fields [1]; for further details, see [13] and the review [2]. Beyond the linearized approximation, it has been used in amplitude computations [9,10,28,29] and related recent works [38,40]. Most exact solutions found so far, however, have been given in other gauges argued to be equivalent to Vasiliev gauge; for example, see [26,37,41]. 12 The resulting manifestly Lorentz covariant form of the master field equations can be found in [26,28]. where have Taylor expansions in (Y α , Z α ) around Y α = Z α = 0 in terms of Lorentz tensors. The redefinition induces a shift symmetry that can be used to set the coefficient of y α y β in W µ to zero, such that where W ′ µ consists of a spin-one field and a tower of higher spin gauge fields with s = 3, 4, . . . . Proceeding by assuming that e αα µ defines a vierbein, and taking Φ| Z=0 and W ′ µ to be weak fields in which the couplings in Eqs. (2.45)-(2.46) can be expanded perturbatively, one can show that the resulting algebraically independent fields are given by the Lorentz scalar ..an is the coefficient in W ′ µ of (σ a 1 ) αα y αȳα · · · (σ an ) αα y αȳα . These fields obey equations of motion on the Lorentzian manifold (X 4 , g µν ) with second-order kinetic terms, critical masses and dynamical metric. 13 The virtue of Vasiliev gauge is that the metric and the gauge fields (2.66) are identical to the Fronsdal tensors that can be obtained at the linearized level by integrating the generalized Weyl tensor using the generalized Poincare lemma (for example, see [43][44][45]). In other words, an asymptotic observer who sources the bulk using a linearized spin-s Fronsdal field will activate the corresponding component field given above, whose boundary value can thus be identified with a dual conformal field theory source coupled to a conserved spin-s current.
The higher order couplings depend on the choice of gauge as well as the initial data for Φ and W µ ; as proposed by Vasiliev [22], these initial data can be fine-tuned at higher orders in order to obtain quasi-local equations of motion in the gauge (2.60).
An alternative approach, which we shall follow here, is to restrict the initial data for the zero-form to specific classes of functions on Y 4 , corresponding to associative subalgebras of A leading to well-defined field configurations obeying physical boundary conditions on M.
13 Whether the resulting system admit any consistent truncation to a pure higher-derivative gravity theory remains an open problem.

Internal star product algebras and solution spaces
where ν belongs to a parameter space and G is a gauge function, obeying the master field equations form an admissible solution space if they generate a free differential algebra together with I and I (for each fixed value of ν). To construct such spaces we use associative star product algebras 14 that are closed under the actions of π,π, † and star multiplication by κ y andκȳ, and whose basis elements T λ , labeled by λ in a discrete set S, have finite traces. We say that A S is the elements in A S can be expanded in terms of the elements in A S ′ in a way compatible with the trace operation.
Expanding the master fields over A S yields a set of modes on X 4 and Z 4 that forms a free differential algebra together with j z and its hermitian conjugate. Using Cartan integration methods, the modes can be expressed locally in terms of zero-form integration constants, which define the ν parameters, and gauge functions. These data can then be adapted to boundary conditions, which may require a change of basis from A S to a basis A S ′ containing A S ; for example, in asymptotically anti-de Sitter spacetimes, it makes sense to impose boundary conditions in a Lorentz covariant basis adapted to a dual conformal field theory. We shall say that a subalgebra A S yields a higher spin gravity solution space if the resulting Lorentz covariant master fields in Vasiliev gauge have symbols defined in normal order that can be expanded over finite regions of X 4 in terms of the set of monomials on Y 4 × Z 4 that vanish at the origin of Y 4 × Z 4 , i.e. they are real-analytic at this point.
The resulting moduli spaces can be coordinatized by higher spin invariant functionals, playing the role of classical higher spin observables [28,32,38,40]. By choosing a structure group [32] and fixing a topology for the base manifold, one may extend the locally defined solutions to globally defined higher spin geometries supporting various types of topologically nontrivial observables. Working locally on X 4 , the accessible observables are on-shell closed zero-forms on X 4 given by combined integrals over Z 4 and traces over Y 4 of adjoint constructs built from (Φ, V α ,Vα; I, I; κ,κ), referred to as zero-form charges. Evaluated on solutions that are asymptotical to anti-de Sitter spacetime, these observables have been shown to have a physical interpretation as generating functionals for correlation functions of holographically dual conformal field theories.
We remark that various subalgebras of A can be obtained from different quantum mechanical systems in four-dimensional phase space. It is an interesting problem to examine which of these are admissible in the above sense, and to furthermore distinguish between these systems using higher spin invariant observables.

New class of biaxially symmetric solutions
In this section, we construct a new class of exact solutions to Vasiliev's equations on a direct product manifold of the form (2.7) using a gauge function and expansions in terms of exponentials of two Cartan generators of sp(4; C), which leads to biaxial symmetry.

Gauge function
From Eq. (2.45) and the fact that X 4 is commuting, it follows that U µ can be expressed in terms of a gauge function G defined locally on X 4 × Z 4 . Thus, setting 15 Eqs. (2.46) and (2.47) reduce to

4)
i.e. the primed fields are constant on X 4 , and Eqs. (2.48)-(2.50) take the form In order to obtain solutions that are asymptotic to AdS 4 , we choose 16 where L, which we shall refer to as the vacuum gauge function, is a locally defined map from and H is determined by imposing the Vasiliev gauge condition (2.60), viz.
in a perturbative expansion where the superscript (n) denotes an n-linear function of Φ ′ . Thus, the master fields in Vasiliev gauge are given by perturbative corrections of where the Maurer-Cartan form U consists of the frame field and Lorentz connection on the anti-de Sitter background spacetime, for which we shall use the explicit form in stereographic coordinates given in Appendix C. As for H, its existence requires that V (L) α admits a power series expansion on Z 4 around Z α = 0, to be examined in more detail in Section 5.
Thus, the dependence on X 4 arises via the gauge function, leaving X 4 -independent equations (3.5)-(3.7), to which we turn next.

Exact solutions in holomorphic gauge from abelian group algebras
One class of solution spaces arise from star product algebras where λ = (λ 1 , . . . , λ N ) belongs to an N -dimensional lattice Λ and for c, π : Λ → Λ. The second relation, which is equivalent to the bosonic projection ππ T λ = T λ , makes it possible to decompose under where a σ; λ;α andǎ σ; λ;α are holomorphic functions on Z 4 and are constant over Y 4 , which may be viewed as a gauge choice (for given zero-form initial data). Expanding 17 and introducingμ where the element γ obeys we obtain two decoupled systems of the form that can be solved using the method of [26] (see also [35]), drawn from the original method devised in [39]. Omitting discrete moduli which arise via projector algebras on Z 4 , two particular solutions that we label by ς = ±1, are given by where and where u + and u − are a set of spinor basis vectors obeying (3.31) 17 We use a convention such that if B ′ = b0 then µ σ; λ = b0ν σ; λ andμ σ; λ = b0ν σ; λ .

Using (B.4), we can choose
The original twistor space connection can thus be obtained by expanding the confluent hypergeometric function in a power series, followed by identifying powers of ζ and γ, though in what follows we shall mainly work directly with the generating functions.

Twistor space connection in Weyl order in holomorphic gauge
The twistor space connection V ′ α is given in the holomorphic gauge by (3.19). From Eq. (2.27), it follows that that is, the symbol in Weyl order of V ′ α is given by the argument of the Wigner map on the right-hand side. This quantity contains singular distributions on Y 4 , which we shall examine in more detail later, and on Z 4 , which we shall examine in what follows. To this end, we observe that the integrand in (3.28) has potential divergences at τ = 0, where log(τ 2 ) goes to infinity, and at τ = −1, where denominators vanish.
As for the potential divergence at τ = 0, it does not lead to any non-real-analyticity in Z 4 to any finite order in perturbation theory as follows from the fact that 18 for τ ∈ [−1, 1], while the same quantity is bounded for τ ∈ [−1, 1] if Re(ςμ ± σ ) 0. Thus, at τ = 0 there is no singularity as long asμ ± σ lies inside the unit disc; indeed, forμ ± σ sufficiently close to zero, the power series expansion of the confluent hypergeometric function yields a basis of functions of τ that can be used to convert the integral equation, obtained by inserting Eq.
(3.28) into the deformed oscillator equation, into an algebraic equation for symbols (for details, see [26,39]). Thus, in order for (3.28) to provide a solution, there has to exist an annulus of convergence in the ζ-space for the Laurent series definingμ ± σ where its modulus is less than one, which can be achieved by tuning the overall strength of the ν-and b n -parameters. In other 18 The confluent hypergeometric function 1F1(a; b; x) := ∞ n=0 (a)nx n (b)nn! obeys 0 < 1F1(a; b; x) < e x for b > a > 0 and x > 0. Its asymptotic form for large |x| is given by 1F1(a; words, the contribution to (3.28) from the region around τ = 0 is real-analytic on Z 4 to any finite order in perturbation theory.
Turning to the divergence at τ = −1, it induces a simple poleå ± α | pole inå ± α at z α = 0, which can be extracted using the formula and analytical continuation of U βγ z β z γ . It follows that Indeed, taking the exterior derivative of the right-hand side one obtains a delta function on the holomorphic slice of Z 4 that cancels the linear source term in Eq. (3.27). As for the higher order corrections toå α in the ν-expansion, they are finite but not analytic at z α = 0, given by combinations of positive powers and logarithms of z α .
As we shall see in Section 5, the nature of the twistor space connection as a distribution on and those to Φ (L) are given by ν σ; 0 Π σ ⋆ κ y . They are hence singular at Z α = 0 and Y α = 0, respectively. Thus, in order for a star product algebra to give rise to proper higher spin gravity configurations, it cannot contain the constant symbol on Y 4 ; in the case of a group algebra this can be achieved by a truncation to a proper semigroup algebra (without the unity), as we shall analyse in more detail in Sections 4 and 5.
In the remainder of this section, however, we shall proceed with the construction of solution spaces in the holomorphic gauge without truncating the underlying group algebras.

Abelian group algebra from Cartan subalgebra of sp(4; R)
In what follows, we shall give an explicit example of a solution space of the type introduced above in the case when the lattice is two-dimensional, i.e. λ = (m,m) with m,m ∈ Z. The underlying group algebra C[Z × Z] is realized as respectively, using sp(4; R) valued matrices obeying 20 from which it follows that det δ α β + aE α γ J γ β = 1 − a 2 2 . (3.43) As for the parameters, we take 19 Inequivalent exact solution spaces can be obtained by replacing E and J by Cartan subalgebra generators in sp(4; C), which we leave for future work. 20 We have suppressed the dummy indices, which are contracted using the north-west to south-east convention.
In what follows, we make the convention that all boldfaced quantities depend on mθ andmθ.
The symbol of T m,m is thus given by By performing the Gaussian integrals, we obtain We also need the symbol of T m,m ⋆ κ yκȳ , which is given by where K αβ K −1 βγ := δ α γ . (3.52)

New exact biaxially symmetric solutions in holomorphic gauge
The above construction of A E,J thus allows us to solve equations (3.5)-(3.7) using the Ansatz (3.18)- (3.19). In order to keep matters simple, we shall assume thatν =α = 0, and work with the following reduced version: where thus ν σ;m,m ∈ C and we recall that the twistor space connection is (anti-)holomorphic, as indicated above. From it follows that the reality condition (Φ ′ ) † = π(Φ ′ ) implies that By definition, the symmetries of the solution are generated by generalized Killing gauge Locally, the space of such parameters is given by where the parameters are arbitrary star polynomials in E and J; and globally, a Killing parameter belongs to an adjoint section obeying suitable boundary conditions, and we shall assume that ǫ (G) is real-analytic on Y 4 × Z 4 and falls off at infinity of X 4 , such that they leave the background spacetime invariant. This implies that the solutions have time-translational and rotational symmetries generated by E and J, respectively. Furthermore, if the Ansatz is expanded over only T m,0 or T 0,m , respectively, then the symmetry is further enhanced to the enveloping algebras of so(2) E ⊕ so(3) or so(1, 2) ⊕ so(2) J , where so (3) is the subalgebras of sp(4; R) commuting to E idem so(1, 2) and J. Acting on the solutions with the full higher spin algebra leads to an orbit that forms a higher spin representation space. The trace operation Tr ′ equips this space with an indefinite sesqui-linear form, as we shall comment on below in the context of higher spin invariant functionals.

Weyl zero-form and Weyl tensors
In this section we compute the Weyl zero-form, Weyl tensors and higher spin invariants formed out of them.

The Weyl zero-form in L-gauge
The Weyl tensors in the L-gauge are contained in the zero-form master field. From (3.13) and To compute T L m,m we use the lemma where L α β is a matrix that depends on the spacetime coordinates (see Appendix C for an explicit expression). It follows from (3.49) that where Under Y α = y α ,ȳα , the above matrices decompose into E L αβ =: where v L 2 := 1 2 v L αβ v L αβ , and idemκ L , which are derived from general properties of any 2×2 symmetric matrix. Furthermore, from (3.41) it follows that which in its turn implies which will be useful later when we determine the Petrov type.

Petrov types of the Weyl tensors
In what follows, we analyze in a few special cases whether the Weyl tensor as the sum (4.21) is of Petrov type D.

Asymptotic behaviour of the Weyl tensors
By using the gamma matrix realization in Appendix B and the global coordinates in Appendix C, we can investigate the asymptotic behaviour of the Weyl tensors. When r → ∞, we have Then the terms in (4.21) of spin-s Weyl tensor at large radius, by a simple power counting, scale as 31) and hence each term is either Kerr-like (when B 2 = C 2 ) or 2-brane-like (when B 2 = C 2 ) in the asymptotic region. The Weyl tensor as the sum of these terms falls off as 1 r s+1 , which is the regular boundary condition of asymptotically AdS 4 solutions.

Zero-form charges
Although the separate spin-s Weyl tensors blow up at the origin of spacetime, the limit of the full Weyl zero-form remains well-defined as the symbol of an operator. From this operator, it 25 Note, however, the matrix K αβ ≡ BE αβ + CJ αβ in this special case has determinant B 2 − C 2 2 = 0, which has consequences for the twistor space connection; see Eq. (5.7). 26 See Appendix D for details on spin-2.
is possible to obtain higher spin gauge invariant quantities given by (4.32) which are referred to as zero-form charges [37] and that are related to higher spin amplitudes [28][29][30]. On our exact solutions, i.e. by substituting (3.13) and (3.53), these charges are given by I 2p | on-solution := 1 32 where A m,m := sech(mθ) sech(mθ) The simplest case is p = 1: In [24], this zero-form charge has been proposed to be one of the contributions to the effective action for higher spin gravity in asymptotically anti-de Sitter spacetimes. As noted at the end of Section 3.6, the resulting contribution to the free energy functional is not positive definite.

Twistor space connection
In this section, we first compute the twistor space connection V (L) α , and show in special cases that it admits a regular power series expansion on Z 4 around Z α = 0 over finite regions of spacetime provided that the group algebra C[Z × Z] is truncated down to a non-unital subalgebra. We then demonstrate the existence of the linearized gauge function H (1) taking the linearized twistor space connection to Vasiliev gauge in a special case.

Generating function for twistor space connection in L-gauge
In order to facilitate the analysis, we write where ρ α is an auxiliary commuting spinor, and we denote j σ (τ ) ≡ j σ (ςμ σ ; τ ). Thus, if these two integrals are finite for boundedμ σ and finite ρ α , then V α is real-analytic in Y 4 × Z 4 .

Singular twistor space connection in
From the discussion in Section 3.3 and the fact that T L 0,0 = 1, it follows that (5.1) contains a term given by (a σ,0,0 ) α , which is not real-analytic in Z 4 . Thus, real-analyticity of V In other words, in the original Ansatz we sum over m,m ∈ Z, but due to the requirement of real-analyticity, we instead sum over only positive m and/or positivem. Furthermore, as can be seen from the table (3.57), for compatibility with the reality condition, along with the truncation we must set θ and/orθ to be real.
θ andθ cannot be both zero.
To summarize, in the following table, we give notations to the consistent truncations, and "×" means that the situation either includes the unity or is inconsistent with the reality condition. 27 Removing the unity from the presentation of A also removes a singularity from Φ (L) .

Regularity of twistor space connection in L-gauge for non-unital A E,J
Under the assumption that A E,J does not contain the unity, we proceed by investigating (5.2).

From (3.51) it follows that
where Thus the case of n = 1 is equivalent to the case of n = 0 by replacing B and C with B B 2 −C 2 and −C B 2 −C 2 , respectively, and multiplying an overall factor 1 √ (B 2 −C 2 ) 2 , which is possible provided that B 2 = C 2 i.e. if mθ = ±mθ for all allowed values of m andm (in the non-unital case). This can be achieved by a suitable choice of θ andθ.
For n = 0 we performing the ⋆-products between T L m,m and the Z-dependent exponential in (5.2), which yields The integrand has potential divergencies at τ = 0, τ = −1 and any value for τ where F 2 or M 2 (τ ) vanishes. As analysed in Section 3.3, the potential divergencies in j σ (τ ) at τ = 0 do not spoil the convergence of the integral provided that the ν-and b n -parameters are sufficiently small. Furthermore, since M αβ (τ ) ∼ (τ + 1) −1 as (τ +1) → 0, it follows that both the prefactor and the exponent are bounded at τ = −1.
To facilitate the investigation of F 2 and M 2 (τ ), which are thus functions of mθ,mθ, X 4 and τ , we use the gamma matrix realization in Appendix B and the coordinates for L in Appendix C. We have not succeeded in a complete analysis, but we have been able to cover a few important special cases as follows: The case A +,0 . In this case, we have θ ∈ R\{0},θ = 0, m ∈ Z + , and hence A = sech 2 mθ, B = tanh mθ, C = 0. Using the explicit matrices and spherical coordinates defined in the appendices, we obtain From mθ = 0 it follows that B = 0, and hence F 2 does not vanish except at r = 0. Moreover, since c(τ ) is purely imaginary, the quantity ςc (τ ) Bλr is purely imaginary as well. Thus M 2 (τ ) vanishes iff Thus, in this case the twistor space connection is real-analytic everywhere away from the equatorial plane in the spherical coordinates. 28 The caseθ = 0. Whenθ = 0, we resort to case-by-case investigation. We will only show two examples below.
For example, if we consider the region of small r, i.e. a small spatial sphere around the origin point, we have A valid choice of the parameters is θ ∈ R\{0},θ ∈ iR\{0} and m ∈ Z + , i.e. the truncation A +,±i . With this choice, we have B ∈ R\{0}, C ∈ iR. Then for C = 0 i.e.m = 0, both first leading terms of F 2 and M 2 (τ ) are non-zero. Form = 0, the discussion is the same as the aboveθ = 0 case. To summarize, the one-form field in this case is real-analytic in the small sphere except on the equatorial plane.
For another example, we consider the region of small ϑ, i.e. a narrow cone around the axis of symmetry ϑ = 0, we have A valid choice of the parameters is θ ∈ R\{0},θ ∈ R\{0} and m ∈ Z + i.e. the truncation A +,± .
With this choice, we have B ∈ R\{0}, C ∈ R. Then for r = 0 we have C − iBλr / ∈ R, and thus, with iςc (τ ) ∈ R, both first leading terms of F 2 and M 2 (τ ) are non-zero. The one-form field in this case is real-analytic in the narrow cone around the axis ϑ = 0 excluding the origin point.

Linearized twistor space connection in Vasiliev gauge
Finally, let us check in a special case that it is indeed possible to bring the linearized twistor space connection to Vasiliev gauge by means of a linearized gauge transformation, as described in Section 3.1, viz.
where H (1) is formally given by Note that, as explained around Eq. (5.12), in the caseθ = 0 with the truncation A +,0 , the regularity of the twistor space connection at ϑ = π 2 has not yet been verified in the L-gauge. However, we expect that this problem would not exist in Vasiliev gauge.
This yields a well-defined limit, such that the twistor space connection is indeed real-analytic at z α = 0. If one instead takes the limit z α → 0 under the integral, one ends up with a divergent integral; this divergence cannot, however, be interpreted as any pole or other singularity at z α = 0. Thus, the prescription that we use is the unique one leading to a sensible result. 29 We note, however, that in the holomorphic gauge the corresponding operations commute, and, correspondingly, the twistor space connection is non-real-analytic at z α = 0 in this gauge; see Similarly, we can also calculate for n = 1: n=1;σ,m α λ=1; ϑ= π 2 , t=φ=0; y=ȳ=0 . (5.29) Finally, using the analog of Eq. (5.1) in Vasiliev gauge, i.e. replacing the label (L) with (G) (1) and substituting (5.28) and (5.29), we obtain We note that, unlike the solution for the twistor space connection obtained starting in the holomorphic gauge, which refers to a splitting of z α into z ± as in Section 3.2, the above expression does not refer to any auxiliary spinor frame in , (5.32) which one can readily identify with (5.30) upon using (4.3).

Conclusion
In We have then switched on the spacetime dependence using a vacuum gauge function L. In the resulting gauge, which we refer to as L-gauge, the spacetime connetion describes an anti-de Sitter spacetime. The terms containing the unity of the internal algebra need to be removed, in order for the Weyl zero-form in L-gauge to be real-analytic on twistor space. The resulting generalized spin-s Weyl tensor, which thus obeys the Bargmann-Wigner equation, is given by a sum of generalized Petrov type-D tensors that are asymptotically Kerr-like or 2-brane-like.
For special values of the parameters, including the symmetry enhanced cases, the spin-s Weyl tensor is of generalized Petrov type D.
We have also shown that the twistor space connection in L-gauge, provided that the group algebra is truncated to a non-unital semigroup algebra as summarized in the table ( Thus, more briefly, we have found families of exact bi-axially symmetric solutions in the holomorphic and L-gauges, and we have verified that they can be brought to Vasiliev gauge at the linearized level in a special case, leaving the more general case as well as higher order perturbation for future study.
We end our conclusions by commenting on future directions. We have left a number of technical details unattended, that we would like to examine more carefully. Besides the issues related to real-analyticity of the linearized master fields in Vasiliev gauge, there is the intriguing degenerate case B = C. Moreover, by taking limits for θ and ν-parameters it is possible to make contact with the solutions found in [26], and more general Kerr-like extensions thereof by expanding the fiber subalgebra using a combination of group algebra elements and endomorphisms in Fock spaces.
More generally, we recall that the importance of Vasiliev's gauge at linearized level is that, when combined with normal order, the linearized spacetine connection W has a Yexpansion at Z = 0 in terms of unfolded Fronsdal tensors and the initial data H (1) | Z=0 modulo gauge transformations. 30 Exact solutions, however, are easier to find in Weyl order using the gauge function method. As far as we can see from the results here and elsewhere, we expect there to be an agreement at the linearized level between the holomorphic and Vasiliev gauges for a fairly large class of linearized zero-form initial data Φ ′(1) (Y ), and it would be desirable to establish this correspondence more precisely, e.g. by expanding Φ ′(1) (Y ) in terms of twistor space plane waves.
Turing to higher order perturbations, the next step is to compute the first subleading corrections to all master fields in Vasiliev gauge, and examine whether real-analyticity in twistor space for generic spacetime points constrains the initial data Φ ′(n) (Y ) for the zero-form and H (n) | Z=0 for the gauge function, for n = 1, 2. This may lead to modified asymptotic boundary conditions in AdS 4 and corresponding corrections to the zero-form charges. In particular, as proposed in [24], the zero-form charge I 2 is a contribution to the free energy functional. The corresponding sesqui-linear form is not definite on the representation space of the underlying higher spin symmetry algebra containing the initial data of our solutions. There are additional contributions to the free energy, however, that may lead to an interesting phase diagram.
The above analysis can also be performed for the closely related Kerr-like solutions outlined above. More generally, one may consider relaxing the Vasiliev gauge as well as the smoothness conditions in twistor space, which may lead to more general noncommutive geometries with interesting properties.

A The ⋆-exponent
Let Y α , α = 1, . . . , N , be oscillator variables obeying where N is even and C αβ is invertible. Denote where K αβ is a constant matrix obeying where indices are raised and lowered using the conventions Y α = C αβ Y β , Y β = Y β C βα , and and to proceed we will compute the symbol of w ⋆ g. To do so we use the identity Hence (A.2) can also be written as Using this, we can show that Proof.
The last two terms can be further converted: Thus (A.9) is proven.
By substituting (A.9), (A.5) can be converted to This differential equation can be solved by substituting the Ansatz which gives 12) and this equation requires the following set of ordinary differential equations for a (t) and b (t) to be satisfied: The general solution is given by 16) where C 1 and C 2 are constants. These are determined by requiring that

B Van der Waerden symbols and gamma matrices
To simplify some of the calculations in this paper, one can use a set of explicit matrix expressions of Pauli matrices and gamma matrices, which for example is given in this appendix. where η ab = diag{−1, 1, 1, 1}.
We further use

D Determination of Petrov type of spin-2 Weyl tensor
In this appendix, we briefly explain how to check (only for spin-2) the Petrov type of a Weyl tensor by using the eigenvalue method. For more details on this topic one can check [47].
The restricted Lorentz group SO + (3,1,R) is isomorphic to SO(3,C), and a Weyl tensor can be converted into its equivalent form with SO(3,C) indices. We can convert the Weyl tensor The indices I, J = 1, 2, 3 should be raised or lowered by the Kronecker delta, so whether they are upper or lower indices does not make a difference.
If we treat Q as a 3×3 matrix, then observing its eigenvalues and eigenvectors is sufficient for determining its Petrov type. Below we list all Petrov types and their corresponding Q-matrix criteria: Petrov types Q-matrix criteria In the list, λ 1,2,3 , λ and − 1 2 λ are eigenvalues of Q, λ 1 +λ 2 +λ 3 = 0 and I is the identity matrix. In particular, being Petrov type D means the matrix Q has three independent eigenvectors while two of them correspond to equal eigenvalues.
Using the explicit matrices and coordinates provided in Appendices B and C, for spin s = 2 we can evaluate the Weyl tensor (4.21) at a given spacetime point with a chosen set of parameters, and then we can evaluate the corresponding Q matrix to check its Petrov type. We have found that in general the Q matrix has three distinct eigenvalues (type I) and thus is not of type D, unless we choose some special parameters or consider only some special spacetime locations.