Higher spin matrix models

We propose a hybrid class of theories for higher spin gravity and matrix models, i.e. which handle simultaneously higher spin gravity fields and matrix models. The construction is similar to Vasiliev's higher spin gravity but part of the equations of motion are provided by the action principle of a matrix model. In particular we construct a higher spin (gravity) matrix model related to type IIB matrix models/string theory which have a well defined classical limit, and which is compatible with higher spin gravity in $AdS$ space. As it has been suggested that higher spin gravity should be related to string theory in a high energy (tensionless) regime, and therefore to M-Theory, we expect that our construction will be useful to explore concrete connections.


Introduction
Vasiliev's higher spin gravity (HSGR) is a gauge theory whose spectrum contains an infinite tower of massless higher spin and matter fields which enjoy infinite dimensional higher spin symmetries. In the tensionless (high-energy) limit of string theory its spectrum is spanned by massless modes which exhibits higher spin symmetries [1] (also in the first quantized level [2]), expectedly equivalent to those of HSGR [3]. The latter arguments have led to the idea that higher spin gravity might describe a highenergy/tensionless-limit of string theory (see e.g. [4][5][6][7][8][9]), or reciprocally that string theory might appear as a broken symmetry phase of some version of higher spin gravity. This also suggests that HSGR is a manifestation of M-theory.
As for matrix models (MMs), it has been argued that M-theory should be described by a matrix model [10,11] in 0`1 dimensions. It has been also conjectured in [12] that certain matrix model in zero dimension should describe an effective theory of type IIB strings. Both statements suggest that matrix models might be useful for the description of non-perturbative aspects of strings and M-theory.
In this paper we show that Vasiliev's HSGR in D spacetime dimensions can be regarded as a type of relativistic MM in D dimensions and we shall argue that certain modification of Vasiliev's equations which incorporate the action principle of matrix models will permit the construction of a new type of higher spin gravities. The latter, dubbed here higher spin matrix models, will be straightforwardly related to matrix models. In particular we shall construct the type IIB higher spin matrix models which combines vanishing (higher spin gravity) curvature and constant (matter fields) covariant derivatives conditions with type IIB matrix model equations of motion [13].

Matrix models and higher spin gravity
There is an interesting connection between super-Yang-Mills and string theories in ten dimensions proposed in [12]. The pure interaction part of super-Yang-Mills, extended with a chemical potential term, obtained after compactification of all directions to size zero is given by the MM action, rA I , A J srA I , A J s´1 2ψ γ I rA I , ψs˙`β1 F , I, J " 0, 1, ..., D´1. (2.1) Here x¨y is the (super)trace of the theory, which respects the correspondent cyclic properties for boson (A) and fermion (ψ) fields and α and β are some constants. When the size of the matrices A and ψ goes to infinity they may be regarded also as Hilbert space operators and therefore they must possess a classical limit, where they become functions on a classical phase-space and their commutator product will become a Poisson bracket. The classical limit of the action (2.1) can be identified with a string action treating the worldsheet as a phase-space and the fields A and ψ as functions on this phase-space. The conjecture of the authors [12] says that their (IKKT) matrix model is related to the type IIB Green-Schwarz action in the Schild gauge [12,14], where Y I are the spacetime coordinate components of the string and ψ is an (anti-commuting) Majorana spinor. Here the Poisson bracket is defined as, and a´g pσq is the determinant of the induced metric on the worldsheet, which can be regarded also as an independent scalar field.
The equations of motion obtained from (2.2) are rY J , rY I , Y J s P s P " 0 , γ I rY I , ψs P " 0, and from the the variation of a´g pσq, Note that the currentψγ I ψ vanishes owing to anti-symmetry of the fermion product ψ α ψ β and the symmetry of the matrices pCγ I q αβ " pCγ I q βα , where C is the conjugation matrix, and γ I are in the Majorana representation. From (2.3) we can solve for a´g pσq, and replacing back its value in the action (2.2) it returns the Nambu-Goto form (see e.g. [15]), which is equivalent to This means, Note that in a flat background the Poisson bracket is equivalent to the determinant of the induced metric on the workdsheet, so that (2.4) can be turned int the standard form found in many textbooks.
A variation of (2.1) more directly related to the string action (2.2) in the classical limit was proposed in [13], (2.6) where now φpσq is the matrix associated to ?´g . The equations of motion obtained from (2.6) are, where tX, Y u :" XY`Y X, is the anti-commutator. The equations of motion of the IKKT model (2.1) are obtained when φ " constant.
The matrix models (2.1) and (2.6) are in particular related to string theory but more generally, matrix models can be regarded as describing subspaces of non-abelian algebras (of functions in noncommutative spaces) by means of constraints obtained from an action principle. The action of a generic matrix model has the form, SrY s " xLrY pyqsy, (2.8) where x¨y is the (super)trace operation, L is a functional of functions Y pyq defined on the basis of the algebra Y Q y. From the variation of Y , δSrY pyqs " xδY CrY pyqsy " 0 produces the equation of motion, It is easy to show that this equation enjoy the transformation symmetrỹ where gpyq and its inverse g´1pyq are functions of the generators y's. It is enough for this purpose to asume a polynomial form of the constraints (2.9).

From matrix models to higher spin gravity
Matrix models can be extended to fiber bundles locally equivalent to MˆUpYq, with local sections given by a set of functions Y py; xq, i.e. which at points x " tx 0 , x 1 , ..., x D´1 u P M are expanded in the basis of an associative algebra UpYq constructed from polynomial functions of the generators of a Lie (super)algebra Y. If the representation of the algebra Y is given in terms of finite dimensional matrices then UpYq will consists of a general linear algebra or some of its associative subalgebras. For more general representations, including infinite dimensional ones, UpYq will be equivalent to the universal enveloping algebra of Y. We shall focus in the latter case. This is, the local sections can be expanded as, where y α P Y is an element of the Lie (super-)algebra Y and their symmetric products y αpnq :" y pα1¨¨¨yαnq P UpYq are given in Weyl order, for definitness, and the sum is up to infinity. Note that considering formal Taylor expansions of the sections Y py; xq in terms of the basis generators y's we can extend the universal enveloping algebra UpYq to non-polynomial classes of functions/distributions (see [16,17] for their use in fractional spin gravity). If y's are finite dimensional the expansion (2.10) will be truncated.
The reader familiarized with HSGR will note already a similarity of (2.10) with the fields of HSGR. Actually, their equivalence is up to determination of the algebra UpYq, which for higher spin gravity it is typically a Weyl algebra of multiple oscillators and their products with Clifford algebras [18][19][20][21][22]. As in HSGR, we shall require local invariance of the matrix model in M-space, i.e. that the systems of constraints (2.9) must be constant with respect to a covariant derivative. To satisfy this requirement we have to introduce a gauge connection, say W , and a constant curvature equation. The complete (integrable) system of constraints is given therefore by where F W :" dW`W^W is the curvature of the gauge connection (one-form) W and Y should be regarded now as a zero-form in the cotangent space of M, and D W " d`W is the covariant derivative, with exterior derivative d " dx µ B µ . These equations can be split in two subsets, the "kinematic" ones (2.11) which involves spacetime (exterior) derivatives d, and the "rigid" ones (2.12) consisting of the algebraic constraints coming from the matrix model (2.8). The system (2.11)-(2.12) is integrable, i.e. it is closed under repeated action of the covariant derivative since, D W CpY q " 0 as consequence of D W Y " 0 and D W D W Y " rF W , Y s " 0 as consequence of F W " 0, etc. The system (2.11)-(2.12) is also gauge invariant under infinitesimal transformations, δW " Dǫ, δY " rY, ǫs.
Solving the kinematic equations can be done absorbing the spacetime dependency in gauge-group elements g such that: Thus the rigid constraints encode most of the dynamics of the system (2.11)-(2.12).
As we mentioned, the equations (2.11)-(2.12) look quite similar to Vasiliev's HSGR equations (see an alternative form of Vasiliev's equations in [23]). The differences are in the details, the algebras involved and the explicit form of the constraints (2.12). Another important detail is the preference for the use of phase-space (deformation) quantization techniques [24,25], i.e. the use of classical functions endowed with a ›-product to construct the associative (non-commutative) enveloping algebras UpYq, instead of working with matrices or Hilbert-space operators (see the review [20] for further details). Though in this section we have not referred explicitly to the ›-product the reader may assume that the non-commutative algebras Y and UpYq are constructed in terms of classical functions and a ›-product.
The fusion of HSGR and MMs is given by (2.11)-(2.12). This means, while the constraints (2.12) encode the dynamics of a MM the kinetic extensions (2.11) will describe HSGR dynamics, i.e. the emergence of interacting generalized Lorentz connections (W ) and matter fields (Y ) with arbitrary spin which will extend standard gauge gravity.
Supporting these ideas, in reference [26] it was noticed that the physical degrees of freedom of higher spin gravity coupled to matter fields in 2`1 dimensions are encoded in a coordinates-free action of the rigid type (2.12), which is actually a matrix model equivalent to a Yang-Mills theory in even noncommutative dimensions. Actually, the equations of higher spin gravity in three dimensions are given by [27,28]: where compared to (2.11)-(2.12) the kinetic constraints are given by (2.15) and the rigid ones by (2.16), while the zero forms are given by Y " tB, S α u. Here α's are spinor indices in three dimensions. Now, according to the authors [26] the system (2.15)-(2.16) can be reduced, using the gauge-function method (2.13), to the spacetime-coordinates-free equations (2.14) now given by the deformed oscillator algebra [29,30], related to (2.16) up to spacetime dependent similarity trasformation (2.13), S α " gs a g´1, B " gbg´1.
Another form of writing (2.17) is, rq, ps " ip1`νkq, tq, ku " tp, ku " 0, which for k 2 " 1 becomes equivalent to the Wigner deformed oscillator algebra [29]. Here q and p are the coordinate and the conjugated momentum of a deformed oscillator and the b field has been factorized in the product of a scalar ν and the Klein operator k. It is worth mentioning that (2.17) implies ospp1|2q (anti-)commutation relations with spp2q " tq 2 , p 2 , qp`pqu and supercharges q and p. This can be extended to ospp2|2q treating k as an internal up1q generator and´ikq and´ikp as additional supercharges. The matrix model action proposed in [26] is given by, S PSV rs, bs " xis α s α b´2b´b 2 y, (2.18) from which (2.17) can be obtained variating b (even degree) and s α (odd degree) and using the correspondent (anti-)cyclic property of the supertrace x¨y. Thus, the equations of higher spin gravity in three dimensions are somehow a covariantization in M space of the Prokushkin-Segal-Vasiliev (PSV) matrix model (2.18). The emergence of higher spin fields on non-commutative geometries has been also observed in references [31,32].
More generally we shall refer as to higher-spin-matrix models to the extension of matrix models and HSGR which can be described by the system of constraints where δS MM rY s " 0, the rigid component of the equations of motion, are obtained from the variation of the action of the matrix model S MM rY s. We say that the higher spin matrix model (2.19) has the same type of S MM rY s.

Deformation quantization of type IIB strings
In the former section we omitted explicit reference to the ›-product. To transform our formulas into a ›-product form we shall assume that for any two functions F py; xq and Gpy; xq of commutative variables x and non-commutative variables y, which can be decomposed as F :" f pxqgpyq " gpyqf pxq, where f pxq and f 1 pxq are differential forms in M-space and gpyq and g 1 pyq are elements of the algebra UpYq, the product F py; xqGpy; xq has inserted a ›-product such that F py; xqGpy; xq :" pf pxqf 1 pxqqpgpyq›g 1 pyqq " pgpyq›g 1 pyqqpf pxqf 1 pxqq.
Let us define now the Groenewold-Moyal ›-product, for two arbitrary functions on the world-sheet f pyq and gpyq. The auxiliary variables ξ and ζ are also world-sheet type and the quantization (Planck) constant θ (with units of area), is the parameter of deformation from the classical juxtaposition product f pyqhpyq " hpyqf pyq. For example the ›-product of world-sheet coordinates is, where y α y β " y β y α is the classical part of the product and the quantum deformation is given by the product of θ and the epsilon tensor ǫ 12 "´ǫ 21 " 1. The ›-commutator is given by rf pyq, hpyqs › :" f pyq›hpyq´f pyq›hpyq, and it is simple to show that in the classical limit and that lim θÑ0 X›Y " XY " Y X.

Now let us deform the string actions (2.2) and (2.4)
. This is simple since these action principles can be expressed in terms of the Poisson brackets. The deformation-quantization version of the Nambu-Goto action is obtained substituting rX, Y s P by rX, Y s › {iθ, so that f " f . Thus in the classical limit, lim θÑ0 S θ NG rY, ψs " S NG rY, ψs.

Deforming the Schild action (2.2) results in
where we have used the values (2.5) and φ´1 is such that φ›φ´1 " 1 " φ´1›φ. In an operator representation, which can be achieved by means of a Wigner map, the equivalent of the action (3.20) is given by (2.6). In the classical limit, lim θÑ0 S θ Sch rY, ψ, φs " S Sch rY, ψ, φs.
The authors of [12,13] assume that the Poisson bracket produced in the classical limit of the respective matrix models is defined in a two dimensional phase space (the worldsheet). This condition may be relaxed since there is no data of the phase space dimension in the matrix model. In the deformation quantization approach we can declare the dimension of the classical phase space in the action itself, as an example of generalization of the model [12,13]. Now the tension T has been rescaled according to the choice of the dimensions of the phase space. Now in the classical limit what is produced is a theory of extended objects. In this case the ›-product will be given by f pyq›gpyq " ż d 2d ξ d 2d ζ pπθq 2d exp´´2 i θξ ζ¯f py`ξq gpy`ζq , (3.22) whereξζ " ξ α ζ α and we have used the symplectic matrix to rise or lower the phase space indexξ " ξ α according to the following conventions, Thus the ›-product of two vectors in phase space yields, The equations of motion obtained from (3.21) are, where tX, Y u › :" X›Y`Y ›X, is the ›-anti-commutator. The analogous to the IKKT model (2.1) is obtained when φ " 1 and d " 1 in (3.22). The equations of motion are obtained from the property for Grassmann parity |A| " 0, 1. dY I`r W, Y I s › " 0, dΨ`rW, Ψs › " 0, dΦ`rW, Φs › " 0, (4.25) rY J , tΦ´1, rY I , Y J s › u › s › " 0 , γ I rY I , Ψs › " 0 (4.26) I Ñ tµrns " rµ 1 µ 2¨¨¨µn s , µ 1 ă µ 2 ă¨¨¨ă µ n , µ k " 0, 1, ..., D´1u. (4.28) The kinetic (2.11) and rigid (2.12) equations are now (4.24)-(4.25) and (4.26)-(4.27) respectively. The rank of indices (4.28) will be explained below. In this specific case, the rigid equations are obtained from the variation of the action for type IIB matrix models (3.21). To the theory (4.24)-(4.27) we shall refer as to type IIB higher spin matrix model.
We have not declared what is the dimension of the spacetime D, neither the rank of the capital indices I, J or the algebras involved. The system above is formally integrable but the choices of the algebras involved and of the dimension of the space-time is an engineering problem, it depends of what we want to describe. In what follows we shall solve the system (4.24)-(4.27) employing similar techniques than in reference [33].
Like in Vasiliev's HSGR, let us consider the Heisenberg algebra, Y " try α , y β s › " iθC αβ ; α, β " 1, ..., 2du, whose universal enveloping algebra, i.e. the algebra of polynomials in Y, is given by UpYq " ty αpnq ; n " 0, ..., 8u, (4.29) where the symmetric products y αpnq :" y pα1 ›y α2 ›¨¨¨›y αnq " y α1 y α2¨¨¨yαn , (4.30) are by the properties of the ›-product equivalent to the classical monomial in the right hand side. In With respect to the ›-commutator product r¨,¨s › the second order polynomials y αp2q generate a representation of the spp2dq algebra, which for 2d " 2 rD{2s contains a representation of the anti de Sitter algebra sopD´1, 2q. The latter algebra can be used to make explicit the Lorentz covariance of the system (4.24)-(4.27). Indeed the variables y αpnq transform in the spin n{2 adjoint representation of the Lorentz algebra in D spacetime dimensions. From this observation it is clear that W contains at level n " 2 the AdS D Lorentz connection and for n ‰ 2 these are higher spin gravity gauge fields, which justify the "higher spin gravity" part of the title of this paper. This is, the AdS D space is a natural solution of the system of equations, if we put all the remaining fields to zero, for example. Now the spacetime dimension is specified to be D, according to the choice of the algebra (4.29).
Now we should clarify the meaning of the indices (4.28). The rigid equations (4.26)-(4.27) are written in a form that reminds us the type IIB matrix model field equations (2.7)-(3.23), but as we shall justify, the fields involved admit more general labels, i.e. with I labeling target space Lorentz multivectors instead of just vectors. As it was shown in [33], the phase space monomials y αp2q (α " 1, ..., 2 rD{2s ), in the classical level, parametrize an algebraic variety (say M) whose coordinates are conveniently labeled by Lorentz multivectors in D dimensions. M admits a covariant non-commutative deformation (quantization), i.e. introducing the ›-product (3.22) in the algebra of functions in the space M, which is Lorentz covariant. Surprisingly enough, the coordinates of the non-commutative version of M will satisfy the constraints (4.26)- (4.27). A posteriori this result justifies the use of multivector labels. Thus we can change the notation Y I by Y rµ1µ2¨¨¨µns ": Y µrns , µ k " 0, ..., D´1, for some values of n to be specified below. In order to observe this let us introduce the exterior algebra of real (Majorana) Dirac matrices, γ µrns :" γ rµ1¨¨¨γµns , µ " 0, ..., D´1, (4.35) constructed from the Clifford algebra tγ µ , γ ν u " 2 η µν , with diagpη µν q " p´1, 1,¨¨¨, 1q. Here n is in the set of integers 1, 2 Mod 4 and n ď D for D even or n ď rD{2s for D odd (see [34]), for which the γ µrns matrices are independent. The matrices (4.35) span a representation of the spp2 rD{2s q algebra, which acts on the space of 2d " 2 rD{2s -components spinors. The symmetric monomials y αp2q also span a representation of spp2 rD{2s q, with respect to the r¨,¨s › product, and indeed they can be given multivector labels using the bi-linear combinations provided by (4.35). This is are spp2 rD{2s q generators.
To make more explicit the correspondence (4.28) I Ø µrns let us split the interval I P t1, ..., dimpspp2 rD{2s qqu in subspaces where I rns are intervals of integer numbers of size`D n˘, in correspondence with the independent components of multivectors X µrns . To each element in the subset I rns we assign a single multivector index I P I rns Ñ µrns " rµ 1 µ 2¨¨¨µn s , µ 1 ă µ 2 ă¨¨¨ă µ n , µ " 0, 1,¨¨¨D´1.

Some solutions
where X µrns is given by (4.36), solve the systems of equations (4.39)-(4.41). It is trivial that Ψ " 0 solves (4.40). If Φ is constant then (4.39) is reduced to, That this equation is solved by (4.36) was shown in reference [33]. For that purpose we can use the identity, and that y α y α " 0.
As we observed, we can unify multivector type of matrix models and higher spin gravity.

Conclusions
We have presented a general class higher spin gravity models which incorporate matrix models in their own definition. The techniques are similar to standard Vasiliev's higher spin gravity, the difference is in the internal algebras and in part of the equations of motion which now are derived from matrix models. These theories put in closer contact higher spin gravity and string theory, as according to e.g. [10,12,13], matrix models are candidates for non-perturbative theories of strings and M-theory. We constructed and provided solutions for a particular model related to type IIB strings, and therefore called type IIB higher spin gravity. The matrix models constructed here extend those found in [33] in that now spacetime dimensions are added. An interesting aspect of these matrix models is that they incorporate coordinates of extended objects, as multivector coordinates of rank k are related to k-dimensional objects. For instance, part of the rigid equations of motion for (multivector) matrix model in 3`1 dimensions (without fermions and for a constant scalar field) can be reduced to the form, Some solutions of this system contain Plucker coordinates of planes through 3`1 dimensions as shown in [33]. An interesting problem to address is whether the classical limit of the respective matrix model (4.24)-(4.27), or (4.24), (4.25), (4.39)-(4.41), reproduces Polyakov's string theory in 3`1 dimensions with fine structure [35].
Note that in this paper we have not studied the details of the physical degrees of freedom contained in the constructed models. This can be done using similar methods than in matrix models and higher spin gravity, i.e. by means of perturbation theory around some solutions of the models, for instance the ones provided in section 4.1.
We would like to mention that multivector extensions of spacetime also appear in the context of E11-type string theories and supergravities [36][37][38] and in the study of higher spin fields dynamics in the approach [39][40][41]. These theories and ours suggest the existence of mutivector extended spacetimes which may be a necessary, or just convenient, for the formulation of quantum theories of gravity.
We expect that our proposal will be useful to find deeper connections between string theories, Mtheory and higher spin gravity. We encourage interested readers to construct and study in detail specific models falling in the category here presented.