Higher-spin initial data in twistor space with complex stargenvalues

This paper is a supplement to and extension of arXiv:1903.01399. In the internal twistor space of the 4D Vasiliev’s higher-spin gravity, we study the star-product eigenfunctions of number operators with generic complex eigenvalues. In particular, we focus on a set of eigenfunctions represented by formulas with generalized Laguerre functions. This set of eigenfunctions can be written as linear combinations of two subsets of eigenfunctions, one of which is closed under the star-multiplication with the creation operator to a generic complex power — and the other similarly with the annihilation operator. The two subsets intersect when the left and the right eigenvalues differ by an integer. We further investigate how star-multiplications with both the creation and annihilation operators together may change such eigenfunctions and briefly discuss some problems that we are facing in order to use these eigenfunctions as the initial data to construct solutions to Vasiliev’s equations.


Introduction
Vasiliev's equations of higher-spin gravity [1,2] are known as an interacting theory with infinitely many higher-spin fields, which can be seen as a non-linear extension of (Fang)-Fronsdal equations of free higher-spin fields [3,4]. General relativity, which in some sense can be seen as a non-linear extension of the spin-2 Fronsdal equation, has much richer physical content than merely spin-2 particles. Likewise, Vasiliev's equations are not just about higher-spin particles, and in order to comprehend various aspects of their physical implications, efforts have been made to find and interpret their solutions [5][6][7][8][9][10][11][12][13][14][15][16]. 1 A frequently used method to solve Vasiliev's equations, which was proposed in [5,20], is that we can first construct solutions in the absence of ordinary spacetime (or more often referred to as "initial data" of solutions), and then turn on spacetime by doing gauge transformations (see [9,21] for recent review and development).
To illustrate this method, let us look at the Vasiliev's equations in 4D spacetime for all integers spins, written in their component equations:

JHEP06(2020)066
where all the fields can depend on three sets of coordinates: x µ for the ordinary 4D spacetime, Y α and Z α for the two 4-dimensional internal and external symplectic manifolds (often referred to as "twistor spaces"), which can be decomposed as Y α = (y α ,ȳα) and Z α = (z α ,zα) with (y α ) † =ȳα and (z α ) † = −zα. The underlined indices are Sp (4,R) or USp(2,2) indices for the AdS or dS background, and {α, β, · · · } are SL(2,C) indices raised or lowered by the Levi-Civita symbols ε αβ , ε αβ , εαβ or εαβ. The SL(2,C) indices that appear in (1.1) refer to those of Z-coordinates, while the Y -coordinates are building blocks of symmetry algebra generators, which are implicit in (1.1) and whose indices are contracted with the ones of field components. In the equations, Φ is a zero-form, U µ and V α (with its hermitian conjugateVα) are the spacetime and twistor-space components of a one-form field, functioning as gauge fields. The star-product between the fields can be formally defined as f 1 (y,ȳ, z,z) f 2 (y,ȳ, z,z) Furthermore, κ is the inner Klein operator satisfying κ = κ y κ z , κ y = 2πδ 2 (y) , κ z = 2πδ 2 (z) , (1.3) (idem.κ), where δ 2 is the two dimensional Dirac delta function, and π andπ are operations that flip the signs of Y and Z coordinates: π x µ ; y α ,ȳα; z α ,zα = x µ ; −y α ,ȳα; −z α ,zα , π x µ ; y α ,ȳα; z α ,zα = x µ ; y α , −ȳα; z α , −zα . (1.4) For the bosonic model, the fields should satisfy ππ (Φ, U µ , V α ) = (Φ, U µ , −V α ), so that non-integer-spin degrees of freedom are projected out, and depending on the background being AdS or dS, different reality conditions should be imposed on these fields (see [6,15] for details, which we will skip here). The simplest solutions to (1.1) are vacuum solutions. As can be easily seen, if we set all fields to zero except the spacetime gauge field U µ , then U µ becomes a pure-gauge independent of Z-coordinates, and thus with an arbitrary gauge function L (x, Y ) is a solution. Whether such a solution is AdS or dS or something else depends on the expression of L.
To obtain solutions other than the vacua, we can do things in the opposite way: if we trivialize spacetime i.e. set U µ to zero and let Φ and V be independent of x µ , then the first three equations in (1.1) are directly solved and thus we can focus on solving the last three equations without involving spacetime (similar to the method of separating JHEP06(2020)066 variables for solving differential equations). In other words, (1.1) can be solved by using the following ansatz: where the primed fields are independent of spacetime. After we obtain the solution for the primed fields, we can turn on spacetime by properly doing a gauge transformation 2 with gauge function G that enables the extraction of Fronsdal fields. We often write G = L H, where L is the vacuum gauge function in (1.5), so that we can do the gauge transformation in two steps, because doing the L-gauge transformation is much easier than the full gauge transformation G, and with the L-gauge we can already study some spacetime properties of the solution at the linearized level (see e.g. [8,15,21]). Now by substituting (1.6) into (1.1), the first three equations of (1.1) are directly solved and the last three are converted to the same equations with all the fields primed: The last two equations of (1.7) can be directly solved by the ansatz 3 Φ = Ψ (y,ȳ) κ y , and (1.7a) after being substituted with this ansatz can also be solved for a (n) α (see [10] for details, which we will not involve in this paper). Here the V α field is expressed as a star-product Taylor series of Ψ, where we label the coefficient of the nth-order term by a superscript (n), and it is a holomorphic function of the z-coordinates (thus this particular gauge for finding solutions is often called the holomorphic gauge). In this gauge the purely Y -dependent expression of Ψ or equivalently Φ is like the "initial data", i.e. once Ψ or Φ is properly chosen, the solution is certain. 4 With the above framework settled, finding a solution with desirable properties is to a large extent the art of constructing the initial data. A very useful and systematic way of constructing the initial data was given in [10], where the authors constructed on the 2 Due to the π,π-automorphisms of the symmetry algebra, two different adjoint representations exist in Vasiliev's equations, namely the adjoint and the twisted adjoint representations. Φ lives in the twisted adjoint representation whose gauge transformation is modified with a π as shown in (1.6).
3 Ψ lives in the adjoint representation and Φ lives in the twisted one. The κy can be used to convert between them. [7] This ansatz was first explicitly used in [8]. See also [9,15] for details. 4 The solution is certain under the assumption that the expression of the gauge function G can be fixed by demanding a good physical interpretation of the solution. See [21] for a recent discussion.

JHEP06(2020)066
Y -coordinates two Fock spaces, whose number operators are the sum and difference of a pair of Cartan generators of the AdS 4 isometry algebra. Linear combinations of the starproduct eigenfunctions (stargenfunctions) of the number operators are used as the initial data, which after a spacetime-dependent gauge transformation can turn into large sets of solutions with spacetime symmetries corresponding to the Cartan generators. To make this short paper self-contained, let us briefly go through some technical details of such stargenfunctions. Let us define the creation and annihilation operators as some constant linear combinations of the Y -coordinates, i.e.
so that [a − , a + ] = 1. The number operator in principle should be defined as a + a − , which is equal to a + a − − 1 2 , but for convenience we use as in [10] the shifted number operator instead (we drop the word "shifted" from now): (1.10) We denote the stargenfunction of the number operator f λ L |λ R (a + , a − ) with λ L and λ R being the left and the right stargenvalues and assume that the function depends on the Y -coordinates only via the creation and annihilation operators, and thus We can make two copies of the above set of creation, annihilation and number operators and stargenfunctions (label them with "1" and "2"), such that w 2 ± w 1 are a pair of Cartan generators 5 of the (A)dS isometry algebra, then a convenient way of constructing the initial data Ψ or Φ , is to let them be equal to the product between the two copies of stargenfunctions i.e. 6 (1.12) For example, in the case that the energy and angular momentum operators are E = 1 2 (w 2 + w 1 ) and J = 1 2 (w 2 − w 1 ) and the stargenfunctions live in the adjoint representation, if we use diagonal stargenfunctions i.e. set λ 1L = λ 1R = λ 1 and λ 2L = λ 2R = λ 2 , which leads to vanishing commutators E, f λ 1 ,λ 2 |λ 1 ,λ 2 and J, f λ 1 ,λ 2 |λ 1 ,λ 2 , then we expect the solution after switching on spacetime by a gauge transformation should exhibit the symmetries of the time-translation and the rotation along a certain spatial axis.
The paper [10] only allowed the stargenvalues λ's to be half-integers and mainly focused on diagonal cases, which seemed to be too much constraint later when we studied 5 The generators are selected among E, J, iB and iP , which are time-translation, rotation, boost and spatial transvection respectively, and for non-compact generators B and P the imaginary factor i has to be multiplied for constructing the Fock spaces. 6 One can prove that the star-product on the r.h.s. is equal to the ordinary product, due to commutativity between the creation, annihilation and number operators from different copies.

JHEP06(2020)066
fields on BTZ-like backgrounds. It is well-known that the 3D BTZ black hole [22] can be obtained from AdS 3 by compactifying the non-compact spatial direction corresponding to a spatially transvectional isometry [23]. Field fluctuations in the context of Vasiliev's equations on such a 3D background were studied in [18], and in [16] we tried to study them in 4D instead with a BTZ-like background obtained from AdS 4 by the same kind of compactification [24][25][26]. Such compactification leads to a periodicity condition on the fields, which leads to the requirement in the initial data that the stargenfunctions of the spatial transvection generator iP should in general be non-diagonal with complex stargenvalues, and furthermore, the excitation of the (angular) momentum along the compactified direction corresponds to adding quantized imaginary numbers to the stargenvalues in the initial data. This motivates us to systematically study the solutions to (1.11) in a more general sense with arbitrary complex left and right stargenvalues and to investigate how far we can go to organize these stargenfunctions into Fock spaces. In section 2, we present solutions to (1.11) in terms of special functions, which are converted into integral representations in section 3 for convenience of doing star-products. Then in section 4 we focus on the set of stargenfunctions expressed with generalized Laguerre functions and investigate how they transform by star-multiplying creation and annihilation operators to generic complex powers. In section 5 we summarize the results and briefly discuss some problems of using these stargenfunctions to further construct valid solutions to Vasiliev's equations.

Solve for the stargenfunctions
In order to solve (1.11), we first use the lemmas and (1.9) to convert the star-products with number operators into derivatives w.r.t. the creation and annihilation operators, which converts (1.11) into or equivalently by taking the sum and the difference: The solution to (2.3) can be written as

JHEP06(2020)066
Note that the two equations stay the same by exchanging {a + ↔ a − , λ L ↔ λ R }, corresponding to the exchange of the solutions f + and f − above. By substituting (2.5) into (2.4) we get (2.6) For generic choices of the stargenvalues, solutions to (2.6) can be written as where the L is the generalized Laguerre function, and the C's are integration constants.
Then we can substitute either one of (2.7) into (2.5), which gives the solution to the original equations (2.8) Because the Tricomi confluent hypergeometric function we can alternatively write (2.8) as (2.10) There exist subtleties that, for special stargenvalues, the two branches of solutions presented above may degenerate. For example, the two terms in (2.10) degenerate when λ R + 1 2 ∈ Z because of the factor cos (πλ R ) in (2.9). For another example, the two terms where (·) · is the Pochhammer symbol. Note that when both stargenvalues λ L and λ R are half-integers, (2.8) and (2.10) reduce to the situation discussed in [10], and in [16] we only focused on the situation that JHEP06(2020)066 either stargenvalue is complex with the other still being a half-integer. This was because we wanted to limit ourselves to the small closed contour integral representation of the stargenfunctions. However, in this paper, we let both stargenvalues be complex numbers in general, and use different integral representations as shown in the next section.

Integral representation
To enable star-product computation of special functions, we very often need to represent them by integrals with integrands where Y -coordinates only appear in exponents, which we shall do in this section. Note that usually the integrals cannot cover the whole parameter space, and the prescription of analytic continuation is needed for the rest of the parameter space after the integrals are done.
We first rewrite the generalized Laguerre function and the Tricomi confluent hypergeometric function in terms of integrals on the real axis To be precise, (3.1) holds when Re (λ R ) < 1 2 and Re (λ L ) > − 1 2 , and (3.2) holds when Re (λ L ) > − 1 2 and Re (w) > 0. We can further convert them into contour integrals on the complex s-plane. Let us define the contour integral (figure 1) where is an infinitesimally small positive number, and where Ln uses the principal branch of the phase angle (−π, π] of its parameter. With this convention we can prove that Furthermore, for the half-circle integral we can do a change of the variable s = e iθ −1, then which, by counting the power of , obviously vanishes in the limit → 0 when Re (λ R ) < 1 2 . Therefore, using (3.5) and (3.2) we can derive Let us also define the contour integral (figure 2) where 1+i 1−i is a half-circle surrounding the right side of 1 on the complex plane. Here we adopt the branch cuts that are consistent with Ln (s − 1) and Ln (s + 1) . (3.9) Then we can derive |s + 1| λ R + 1 2 ds , JHEP06(2020)066 ds . (3.10) Again, we can see with a change of the variable s = e iθ + 1 that the half-circle integral vanishes in the limit → 0 when Re (λ L ) > − 1 2 . Therefore, combining (3.10), (3.1) and (3.2) we can derive The second term on the r.h.s. is the result of −1+i , and it vanishes when λ L − λ R ∈ Z, which can be explained by the cancellation of the branch cuts where s < −1 on the real-axis.
To summarize, we can rewrite the solution (2.10) as and it is important to notice that the convergence of the above contour integrals requires only Re(w) > 0 and has no restrictions on λ L,R . In the diagonal case λ L = λ R , due to the cancellation between two branch cuts, C −∞,1 can be replaced by C −1,1 , which is a contour integral surrounding the section of the real axis between −1 and 1 ( figure 3). At the end of this section, let us briefly compare the contours in this paper and the previous ones in [10,16]. The stargenfunctions in [16] (or [10]) are just special cases in section 2 with one (or two) of the left and right stargenvalues set to half-integers. Therefore, we should be able to reproduce these special cases also from the contour integrals. 7 Take the first term on the r.h.s. of (3.14) for example, if we set λ L = λ R ∈ Z+ 1 2 , all brunch cuts can be dropped, then the contour reduces to two circles around the points s = ±1, one of which has zero residue, and thus the integral reproduces stargenfunctions with Laguerre polynomials as in [10]. Take the second term for another example, if we set λ R + 1 2 ∈ Z + , the left-side brunch cut in figure 1 goes away, then the contour reduces to a circle around s = −1, and thus by keeping λ L ∈ C with the other brunch cut, the integral reproduces stargenfunctions with generalized Laguerre functions as in [16].

Creation and annihilation operators
The set of stargenfunctions (2.8) is the most direct generalization from [10,16]. In this section, we investigate how they change by star-multiplying creation and/or annihilation operators to complex powers. Let us call the stargenfunctions represented by the two terms of (2.8), respectively, the "plus" and the "minus" branches. For convenience of star-product computation, we write the two branches by the integral representation (3.1), i.e. 8 where the constant overall factors have been dropped hereafter for simplicity. Note that the two branches are in general different, although they have the same stargenvalues, and that the two degenerate into one when λ R − λ L ∈ Z as already discussed around (2.11). In this section, we will investigate how these functions (as if they were eigenstates of the number operator in quantum mechanics) can be related by the creation and annihilation operators.

The plus and minus branches
We start from the diagonal case λ L = λ R = λ We would like to show that by doing a left or a right star-product with a ± to a complex power can shift the stargenvalue accordingly, and to this end we first express a ± with the power p using the Mellin transform (as done in [16])

JHEP06(2020)066
so that we can pick out the factor e −τ a ± for the convenience of doing star-products. 9 Then by using (4.1), (4.3) and the following set of star-product results e −τ a + e 2sw = e 2sw−τ (1−s)a + , (4.4a) e 2sw e −τ a + = e 2sw−τ (1+s)a + , (4.4c) we can easily prove that Thus the shift of the left (right) stargenvalue is indeed equal to the power of the creation (annihilation) operator acting from the left (right); or equal to the opposite power if it is the annihilation (creation) operator. Furthermore (4.5) shows that, starting from a diagonal stargenfunction, the creation / annihilation operator always brings it to the plus / minus branch, no matter whether the star-product is done from the left or from the right, and any non-diagonal stargenfunction of the two branches can be reached in this way, up to an overall constant, from a diagonal stargenfunction. Naively we can further use (a ± ) p 1 (a ± ) p 2 = (a ± ) p 1 +p 2 to conclude that the plus / minus branch is always closed under the star-multiplications with creation / annihilation operators to generic complex powers, but here we encounter some subtlety about associativity. Star-products are not well-defined unless associativity holds, and associativity is bound to work if all the formulas to be star-multiplied are Taylor-expandable in the Y -coordinates. However, a ± to complex powers cannot be Taylor-expanded. To circumvent the problem, we need to use integral representations, so that all the a ± 's stay in exponents, which makes the formulas Taylor-expandable. In other words, the integrals that we have used in this paper are not only technical tools to simplify the computation, but also very often the prescription to guarantee associativity. To ensure associativity during the computation, we need to keep in mind that all star-products should be computed before doing the integrals. 10 For example, in order to show we need to express all the factors by the integrals: (1 + s) λ+ 1 2 e −τ 1 a + e −τ 2 a + e 2ws , 9 We impose the restrictions Re λ (L,R) < 1 2 , Re a ± > 0 and Re (p) < 0 for convergence of the above integrals, but after we have finished the star-product computation and completed the integrals, we obtain back formulas like a ± p and the generalized Laguerre functions, then by analytic continuation we can lift these restrictions. 10 See also [16,21] for relevant discussions.

JHEP06(2020)066
then we should do the integrals only after finishing the star-products within the integrand. By doing the star-products, the above formula equals (4.8) and by rescaling τ 1,2 it further equals In this way all star-products have been done between exponential functions of Y , thus the problem of associativity has been circumvented, and we can then do the integrals in (4.9), which gives f + λ+p 1 +p 2 |λ as expected. Using this method we can prove that indeed the plus (minus) branch is closed by star-multiplications with arbitrarily many factors of creation (annihilation) operators.

Mixing the branches
Now naturally we want to ask the question: if both the creation and the annihilation operators together act on a stargenfunction, then what happens?
The first thing to check is whether the result is still a stargenfunction. As a test, we now examine whether is a stargenfunction of w.
Using the lemma we obtain where in the last step we have simplified the formula by rescaling the integration variables τ ± . We denote
In the same way as deriving (2.2a), we convert the w on the l.h.s. of (4.14) into derivatives w.r.t. a ± and hence obtain which leads to We can convert the terms (4.17a) (4.17b) and (4.17c) by integrations by parts, so that the curly bracket above can be replaced with which is exactly equal to the factor (λ + λ + − λ − ), and thus (4.14) is proven. Moreover (4.15) can also be proven in the same manner, whose details we will skip. In (4.10), we can use (4.5) to move one or both of the creation and annihilation factors to the right side, and thus it is clear that a diagonal stargenfunction of the type (2.8), starmultiplied with two factors respectively of creation and annihilation operators to arbitrary complex powers, gives a new stargenfunction with the expected stargenvalues. The new stargenfunction in general should be a linear combination of different branches, but we are not yet able to decompose it. What makes the problem more complicated is that, when the plus and minus branches degenerate, there exists an extra branch of stargenfunctions. Such an extra branch is very likely to be activated by some combinations of creation and annihilation operators. To illustrate this, in (4.10) we set λ + = λ − =λ, which leads to a diagonal stargenfunction, and by performing the τ ± integrals in (4.13) we get Note that the extra branch can be represented (except for the special cases with halfinteger stargenvalues) by the second term in (2.10), which contains a factor U 1 2 − λ, 1, 4w that logarithmically diverges as w → 0. In (4.18), we can find a similar factor U −λ, 1, 2(1−s 2 ) c+s w whose expansion w.r.t. w also gives log w, and we therefore speculate that the extra branch is involved.

Conclusion and discussion
In this paper, motivated by the construction of initial data in the holomorphic gauge to solve Vasiliev's equations, we have studied the stargenfunctions of the number operator, as defined in (1.11) with generic complex stargenvalues. We have re-written these stargenfunctions using integral representations and also expressed the creation and annihilation operators to complex powers by integrals via the Mellin transform, so that their star-products can be computed conveniently. We have in particular picked out a set of stargenfunctions (2.8) to investigate how they are changed by creation and annihilation operators. This set of stargenfunctions can be written as linear combinations of two subsets of stargenfunctions (the "plus" and "minus" branches), except that the two subsets intersect when left and the right stargenvalues differ by an integer. The whole plus (minus) branch of stargenfunctions can be related (up to overall constants) by -and is closed under -the star-multiplications with creation (annihilation) operators to generic complex powers. However, by mixing the usage of creation and annihilation operators, (2.8) does not seem to close, although the stargenvalues are likely to be what we expect. We have JHEP06(2020)066 confirmed that a diagonal stargenfunction of the type (2.8), star-multiplied by two factors respectively of creation and annihilation operators to complex powers, yields another stargenfunction with the expected stargenvalues. We speculate that the resulting stargenfunction should be a linear combination of the plus and minus branches when they are non-degenerate and that another extra branch should be involved in the degenerate case.
We are not yet able to use the above results to construct valid solutions to Vasiliev's equations. One of the difficulties happens when we go beyond linear combinations of the stargenfunctions. As in (1.8) the twistor space gauge field V is expressed in terms of starproduct series of the initial-data field Ψ. Thus, for example if we use the stargenfunctions f ± λ L |λ R above to construct Ψ, we will encounter star-multiplications to all orders between these stargenfunctions, and we have to make sure that such star-multiplications make sense. Let us take a look at the diagonal stargenfunctions to illustrate the problem. When the stargenvalues are half-integers, according to [10], f λ|λ should act like a projector: f λ 1 |λ 1 f λ 2 |λ 2 ∝ δ λ 1 λ 2 f λ 2 |λ 2 . Naively we expect the same thing to happen for generic complex stargenvalues, which however, does not seem to be true.
If we multiply the expression (4.2) by itself, by using the lemma e 2sw e 2s w = 1 1 + ss e 2 s+s 1+ss w , (5.1) and denoting ∆λ ≡ λ 1 − λ 2 , we get where in the second last step the integration variable s has been changed by s = s−u su−1 or equivalently u = s+s 1+ss . In the last line of (5.2), obviously the integral is problematic: let be a small positive number, then we have and as the factors −∆λ and ∆λ indicate, if is set to 0, for any ∆λ the integral does not converge. 11 In the paper [10] for half-integer stargenvalues, a prescription has been given: the realaxis integral JHEP06(2020)066 (both types of integrals can produce the same Laguerre polynomial at the beginning). In that case, while the real-axis integral is problematic, the contour one converges and leads to the expected projector result. On the contrary, in this paper for complex stargenvalues, we have to prescribe with a different contour C −1,1 as shown in figure 3, i.e. we now have the factor: Within a given branch of the function on the r.h.s. of (5.5), 12 by going through a full cycle along the contour C −1,1 , both (s − 1) and (s + 1) acquire the same phase factors that cancel each other, i.e. the value of the primitive function is not changed by a full cycle along C −1,1 . Therefore, the integral (5.4) gives 0. For the special case ∆λ = 0, (5.4) becomes the sum of two cancelling residues around s = 1 and −1, which also gives 0. Therefore, with such a prescription, f λ 1 |λ 1 f λ 2 |λ 2 = 0 always holds, even if λ 1 = λ 2 . 13 It is not yet certain to us how this should be interpreted -perhaps we should for consistency exclude half-integer stargenvalues when we use such stargenfunctions to construct the initial data, or perhaps we should find a better prescription that we are not yet aware of. 14 Another complication is the star-product with the inner Klein operators. As introduced around (1.12), to construct the initial data, we need to make two copies of the stargenfunctions. If each copy has two branches e.g. (+, −), then two copies give four branches (++, +−, −+, −−) in the initial data. A stargenfunction of any particular branch starmultiplied by κ y orκȳ in general may lead to combinations of all branches. Thus the computation is very much involved, and we will continue these investigations in our future work. JHEP06(2020)066 particular, the author is grateful to C. Iazeolla and P. Sundell for pointing out important details about the contour integrals. The author would also like to thank all the organizers and hosts of the APCTP-KHU Workshop on Higher Spin Gravity in Pohang, South Korea, for hospitality and for providing this excellent opportunity for discussions that helped the improvement of this paper.
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.