Chiral Algebras of Two-Dimensional SYK Models

We study chiral algebras in the $\bar{Q}$-cohomology of two dimensional SYK models with extended supersymmetry. In a special limit discovered in arXiv:1805.09325, we are able to construct explicitly a"vertical"single-particle higher-spin algebra that is bilinear in the fundamental fields. This algebra can be regarded as the counterpart, when going away from criticality, of the infrared emergent higher-spin symmetry of the $\mathcal{N}=(0,2)$ SYK model. Moreover, a second"horizontal"single-particle higher-spin algebra appears in this limit. Together with the vertical algebra they generate a stringy algebra with a"higher spin square"structure that is believed to appear in the tensionless limit of string theory. On the other hand, we do not find single-particle higher-spin algebra away from the special limit, which is consistent with the result in arXiv:1805.09325. Our analysis is carried out for each individual realization of the random couplings and for finite $N$ (and $M$), which in particular indicates that the conclusion in arXiv:1805.09325 is robust to $1/N$ corrections.

Besides the gravity sector, there is in addition a tower of operators [6][7][8][9] with finite anomalous dimensions. As suggested in [7], this tower of operators is analogous to the tower of operators in vector models and a further discussion on this relation can be found in [44]. Different bulk duals of the tower of operators are proposed in [11,[45][46][47][48].
Higher dimensional generalizations of the SYK model are interesting in making connections to other previously studied models and even to experimental realizations.
There have been proposals with a discrete spatial direction, see e.g. [49], or with an unconventional kinetic term [50], or with marginal irrelevant interactions [51].
In this paper we study a model, introduced in [1], that is defined in (continuous) 1+1 dimension with canonical kinetic terms and relevant SYK-like random coupling, which is a direct generalization of [52], and see also [41]. The model has an N = (0, 2) supersymmetry. In the infrared, this model is dominated by the set of melonic diagrams in the large-N limit and can be solved perturbatively. The N = (0, 2) supersymmetry is crucial for our discussion. The N = 2 supersymmetry in the right-moving sector makes the IR solution reliable; the absence of supersymmetry in the left-moving sector allows a one parameter family of such models. In certain limits of this family we observe emergent higher-spin symmetries in the infrared [1]. This provides an explicit illustration of a connection between SYK-like models and models with higher-spin symmetry: higher-spin theories can be thought as a subsector of some tensionless limit of string theory [53][54][55][56][57][58], while the SYK model should be holographically dual to some string theory with finite tension [7,59] 2 . In addition, it is explicitly observed that the chaotic behaviors disappear as the higher-spin symmetries emerge, which is consistent with the general picture. Furthermore, the pattern of how massive fields becomes conserved higher-spin currents agrees with the pattern of the breaking of the higher spin symmetries from a previously conformal perturbation computation [60], which provides another evidence to support our claim. The field content of this model, namely N chiral (bosonic) multiplets Φ a (a = 1, 2, · · · , N) and M Fermi multiplets Λ i (i = 1, 2, · · · , M), is identical to those in some two dimensional N = (0, 2) Landau-Ginzburg models on flat spacetime [61]. The only difference is the form of the interaction: the coupling of the SYK model is [44] random and the model does not have any global symmetry except for a U(1) charge symmetry.
In this paper we will discuss a property that is relevant to the N = (0, 2) SYK model, and probably also to the N = (0, 2) Landau-Ginzburg model, which is the chiral algebra of theQ cohomology. At the classical level the quasihomogeneity of the superpotential guarantees [62,63] the existence of the conformal stress tensor and implies that the classical action is invariant under the right-moving U(1) R symmetry. It is easy to verify that the potential of the SYK model is indeed quasi-homogeneous.
At the quantum level, we show that the SYK model satisfies a necessary condition for the existence of the stress tensor. This means there is always a Virasoro subalgebra in the chiral algebra of the 2d SYK model.
In this paper we further study properties of the higher-spin operators, which are those generators of the chiral algebra with spin greater than 2 (and their supersymmetric partners). We would like to understand what is the minimum subalgebra of the chiral algebra that includes at least one, in our case the one with the smallest spin 3 , higher-spin operator. As we will show in the paper, in the special limit q → N M , where q + 1 labels the rank of the interaction of the SYK model (2.1), one can identify two different higher-spin subalgebras that are isomorphic to the W ∞ [λ = 1] algebra studied in [65]. The generators of one of the algebras, which we call "vertical", can be expressed in terms of single-sum operators that are quadratics of the fundamental fields.
Here we denote by "single-sum" the operators involving one sum of the flavor indices of the chiral or the Fermi multiplet, which is an analogue of the single trace operator in matrix models. The generators of the other algebra can be expressed in terms of single-sum terms that are higher powers of the fundamental fields. The commutators of the two algebra generate a larger algebra that has the structure of a "higher spin square" [57,66].
One would also ask what is the relation between the chiral higher-spin subalgebra of theQ-cohomology and the chiral higher-spin algebra emerging in the special limit of the SYK model (2.1) in the infrared [1]. From our analysis in this paper, it is tempting to consider the vertical higher-spin algebra in the limit q → N M -whose generators are all quadratic in the fundamental fields -as the counterpart of the emergent higher-spin symmetry discussed in [1] when going away from the infrared SYK critical point.
However, once we are away from the limit q → N M the inclusion of a single higherspin operator leads to an infinite dimensional subalgebra that is larger than a conventional higher-spin algebra: the number of generators of this algebra at each single spin grows as the spin increases, which is in contrast with the conventional higher-spin algebra where there is one operator (supermultiplet) at each spin. This reflects the stringy nature of this algebra; there is not a simple way to separate out the leading Regge trajectory as opposed to the case in the 0+1 dimensional SYK model. This could also be related with the fact that these 2 dimensional SYK models are not maximally chaotic. Furthermore, the fact that there is no "single-particle" higher-spin subalgebra in the chiral algebra away from the q → N M limit can be understood holographically as there is not a tower of higher-spin fields in the bulk. Therefore this is consistent with the fact that the Lyapunov exponent of the early time out-of-time-order correlation function is not zero away from this special limit.

2D supersymmetric SYK model
We consider SYK models in continuous 1+1 dimensional spacetime that is of the class discussed in [52]. Our primary example of this class of model is the one with N = (0, 2) supersymmetry discussed in [1]. The model describes N chiral multiplets Φ a and M Fermi multiplets Λ i with a random coupling The model with M = N has an enhanced N = (2, 2) supersymmetry and the model reduces to the one discussed in [52], see also [41]. Each superfield contains one bosonic and one fermionic field where z is the holomorphic coordinate. The coupling J ia 1 ...aq is randomly chosen from a Gaussian distribution. It is relevant, with mass dimension one, so it dominates the physics in the infrared. The IR solution of this model in the limit is presented in [1].
As shown in [1], at generic value of µ, this model shares some common features of SYK-like models, such as being chaotic and the emergence of a conformal symmetry in the IR. An intriguing feature of this model is the emergence of higher spin operators in two different limits of the model. In each of the two limits one observes a tower of operators that become holomorphic whose left-moving conformal dimension vanish. There is also a tower of antiholomorphic operators in each of the limits as well. These operators close under the Operator Product Expansion (OPE) and generate a higher-spin symmetry algebra.
The appearance of the higher-spin operators and higher-spin symmetry is confirmed by the vanishing of the Lyapunov exponent in the two limits. We emphasize that such higher-spin symmetries appear only in the two special limits of the IR model, which mimic the limit where the string tension approaches to zero.

Chiral algebra in theQ cohomology
In this paper we discuss a slightly different type of algebra that is present in the model (2.1) for generic µ. This is the chiral algebra generated by the cohomology classes of one of the supercharges. Explicitly, the SYK model we are interested in have two supercharges We are interested in the cohomology of one supercharge, sayQ. Following (2.7), the elements of the cohomology can only havez dependence, up toQ exact terms. As a result, the elements of the cohomology generate a chiral algebra in the antiholomorphic sector [62,67], see also [68]. It is this chiral algebra that we would like to study in the 2 dimensional supersymmetric SYK models. 4 Such chiral algebras have been discussed for Landau-Ginzburg models, see e.g. [62,[68][69][70]. One crucial property of the chiral algebra, which follows from the fact that the algebra is defined in the cohomology ofQ, is that details of the interaction do not affect the form of the chiral algebra [62,69]. In particular, the chiral algebra admits free field representations in terms of the (β, γ)-ghost system for the chiral supermultiplet and the (b, c)-ghost system for the Fermi multiplet. This is also true for our model. A similar argument following [62,[69][70][71] indicates that the chiral algebra of our model can also be constructed in terms of the ghost fields. In the following sections, we explicitly construct some of the generators of the chiral algebra with higher spin. We choose to use the set of φ a ,∂φ a , λ i ,λ i , as the fundamental building blocks. 5 Furthermore, since the interaction term does not directly affect the form of the algebra, such algebras exist for each individual realization of the random coupling. Therefore our result does not rely on averaging over the random coupling, which makes our conclusion applicable to the full quantum mechanical model and avoids possible subtleties about the replica symmetries, see e.g. [72][73][74][75].
On the other hand, although the chiral algebra in the cohomology ofQ does not directly depends on the form of the interaction [69], the form of the interaction (2.1) does determine the global symmetries the model has, which imposes nontrivial conditions on the chiral algebras of the model. Firstly, our model has a U(1) R symmetry that transforms the various fields according to The fact that each term in the sum of the interaction is invariant under this symmetry meansα Since the interaction runs over all combination of the fermions, which means there is one equation (2.11) for each set of indices {i, a 1 , . . . , a q }. This large set of equations is only solved when α a = α ,α a =α , (2.12) 5 Our results can also be written in terms of the ghost fields via the following dictionary Following almost identical argument as in [62,63], this is the same condition to guarantee the existence of a left-moving stress tensor at the classical level.
In addition, there is another global U(1) L symmetry of the model (2.14) The invariance of the action under this symmetry requires the charges to satisfy These relations among the charges give another set of nontrivial constraints on the form of the generators of the chiral algebra.
Secondly, the interaction (2.1) breaks the complete permutation symmetry, and hence any non-abelian global symmetry, among the N chiral multiplets and the M Fermi multiplets. Therefore the operators in the chiral algebra do not need to respect such symmetries. This allows more general form of the chiral algebra, which is the crucial point so that a direct connection to the "higher spin square" is possible.

Enhance to N = (2, 2) supersymmetry
The model with M = N has an enhanced N = (2, 2) supersymmetry, and the model reduces to that discussed in [41,52].
In this case the U(1) R and U(1) L symmetry presented in the N = (0, 2) model In particular, since the supersymmetry is enlarged, the chiral and Fermi multiplets are combined into N = (2, 2) chiral supermultiplet, so their charges are related. In 17) and the charge of the N = (2, 2) supermultiplet under the U(1) v symmetry is simply As a consistency check, at M = N the condition (2.13) and (2.15) leads to The total R-charge under the U(1) v symmetry then becomes which is the same as the R-charge of the Landau-Ginzburg model with a Φ q+1 potential.
Apart from this, the rest computation of the chiral algebra at M = N is largely parallel to those of the N = (0, 2) model which we discuss in the next section.
3 Minimal subalgebras with higher-spin extension In this section we work out the first few operators of the chiral algebra and determine the OPE among them.
As discussed above, the higher-spin operators that generate the chiral algebra can be expressed in terms of the free fields (2.9) of the chiral multiplet and the Fermi multiplet with the OPEs 6 The (antiholomorphic) stress energy tensor of the model (2.1) is with the standard OPE where the central charge is given by With respect to this stress energy tensor the conformal dimensions for φ a (z),∂φ a (z), The α a and α a parameters take the values in (2.13) and (2.15). In particular, the central charge is Notice that our stress energy tensor takes a similar form as those of the Landau-Ginzburg model [63,68].
• Spin-1 There is a spin-1 operator corresponding to the U(1) L symmetry which is a Virasoro primary operator of spin 1 under the stress energy tensor if the following condition is satisfied This together with (2.15) fixesα parameter tõ In the following, we focus on the model with this value of the parameter. In particular, the central charge for this case becomes The level of the U(1) current algebra is determined by the OPE For our SYK model, the charges obey (2.13) and (2.15). This means 11) and the condition (3.7) reduces to The U(1) R current is in the same supermultiplet with the stress energy tensor and does not lead to new operator in the cohomological algebra. There is no other conserved spin-1 generator that leads to a symmetry of the model.

• Spin-3
Consider the following ansatz for a spin-3 operator where repeated indices are summed over. We look for a spin-3 Virasoro primary operator with This fixes some of the coefficients in the ansatz (3.13). Further requiring its OPE with the spin-1 operator to be regular, as well as its OPE with itself closes in the algebra, we obtain further conditions that fix all but three coefficients. One of the remaining coefficient corresponds to the normalization of the spin-3 operator and can be fixed by In the following we choose the canonical normalization, namely c W 3 = (c − 1)/3, after the spin-1 operator is factored out. One can then read off the spin-4 operator from the second order pole of the W 3 (z)W 3 (w) OPE. Further requiring the spin-4 operator to be Virasoro primary, together with the requirement of its OPE with the spin-1 is regular, fixes all but one parameter in the ansatz (3.13). After applying all these conditions, where a new spin-4 operator W 4 (w) appears. In the above expression we have also definedT that have regular OPE with the spin-1 operator J(z). The explicit form of the W 4 operator is not quite illuminating, so we omit it for simplicity. Notice that this spin-4 operator still depends on the remaining undetermined parameter. The above W 3 (z) W 3 (w) OPE agrees with the known result in the literature, say the one in [77] once we identify thec there as c − 1 in our expressions.
One might then propose that we can continue to construct operators with higher spins and to find a W ∞ algebra. However, we claim that the algebra generated in this way is not a conventional W ∞ algebra. The reason is the following. The conditions mentioned above are all that we could impose to get a closed algebra among the first few low spin operators. However, the fact that we cannot fix one parameter in the ansatz (3.13) is quite surprising. Indeed, the unfixed parameter in the spin-3 ansatz signals the existence of two different spin-3 primary operators. The parameter also appears in the spin-4 operator, which indicates that there are at least two spin-4 operators in the algebra. On the other hand, we know from the result in [78] that the family of bosonic W ∞ [λ] algebra is described by only two parameters, which relies on the assumption that the set of operators that generate the higher-spin algebra contains only one higher-spin operator for each given spin. Therefore the fact that there are more than one operators at each large enough spin is consistent with the fact that there are more parameters than the conventional W ∞ algebra. This implies that what we have found is not a higher-spin type algebra. In other words, the minimal subalgebra of this chiral algebra that contains at least one higher-spin primary operator is larger than the conventional higher-spin algebra.
We believe the above statement is correct. A simple test of this claim would be to compute higher order OPEs, such as the W 3 (z)W 4 (w) and the W 4 (z)W 4 (w) OPE, and test if the parameter does appear explicitly in these OPEs. Unfortunately, we do not have enough computational power to confirm this. We hope to come back to this computation in the future with an updated power of computation. Nevertheless, we can test our statement indirectly by considering a special model at M = N, where the model has an N = (2, 2) supersymmetry. In that case due to the larger symmetry, the result is more constrained and we can observe the existence of extra operators, comparing to the conventional SW ∞ algebra in e.g. [79], explicitly. Since the N = (2, 2) is a special case of the N = (0, 2) model, the result in the next section provides a test of the prediction we made on a special case. supersymmetry. This enhancement renders the cohomological chiral algebra supersymmetric. In the following we construct the first few higher-spin operators of this model explicitly.
Due to the presence of supersymmetry, the operators are organized into supermultiplets. The lowest supermultiplet contains a spin-1 current operator, two supercharges and the stress energy tensor. This multiplet of the N = (2, 2) Landau-Ginzburg model is worked out in e.g. [62]. The operators in our model are similar, the only difference is that the charges α a ,α i , p a andp i take the value in (2.17) and (2.19). Explicitly, the operators in this multiplet of our SYK model (2.1) are The central charge is given by The generators J(z), G ± (z), and T (z) satisfy the standard N = 2 superconformal algebra. Next we consider the higher-spin operators. The first higher-spin N = 2 multiplet consists of operators with spin (2, 5 2 , 5 2 , 3) . We only need to determine the new spin-2 operator and the other operators are its N = 2 supersymmetric descendants. The superconformal primary spin-2 operator reads W 2 (z) = (2 + q)(2q − 1)∂φ a∂φa (z) + (2 − q) φ a∂2φa (z) + 1 2 (1 − 2q)q λ i∂λi (z) We expect the inclusion of this multiplet leads to an infinite dimensional algebra [80]. We can proceed to compute the W 2 (z)W 2 (w) OPE where the central charge is (3.22) and The form of this OPE is identical to equation (3.30) of [81], so one is tempting to identify a quartic relation between the λ parameter of the SW ∞ [λ] algebra and the q, N parameter in our model But as we will argue in the following this mapping might not be very meaningful. The problem is the following. One can work out the other components of the multiplet, for example the two fermionic components in the multiplet are where repeated indices are all summed over. One can then compute the W + 5/2 (z)W + 5/2 (w) OPE. If the underlying algebra is indeed a SW ∞ [λ] algebra [79,80], one expects that the right-hand-side of this OPE should only contain normal ordered products of operators with smaller dimension since there is no generator in the SW ∞ [λ] algebra with two units of the spin-1 charge. Surprisingly, this is not the case for our algebra. One can compute the W + 5/2 (z)W + 5/2 (w) OPE explicitly and get where the new Virasoro primary spin-4 operator (3.31) carries two units of the spin-1 charges. One can further check that it is a descendant of a spin-7/2 supersymmetric primary operator with one unit of spin-1 charge It is easy to check that they belong to a multiplet with spin content ( 7 2 , 4, 4, 9 2 ). Therefore, we explicitly observe that the subalgebra we identified, which is the minimal one that contains a higher-spin operator, contains more operators than a conventional higher-spin SW ∞ algebra. In addition, based on the pattern in which unexpected multplets appear, the number of supersymmetric multiplets seems to grow as spin increases.

Anomaly
As discussed in section 2.2, at the classical level, the existence of a stress energy tensor in the algebra ofQ cohomology is guaranteed by the quasi-homogeneity condition of the potential Λ i J i (Φ a ). Now that we have constructed a few operators in the chiral algebra, we can test if the result is anomaly free at the quantum level. We do not know a set of sufficient conditions for this, but we have a necessary condition to check following [63].
We will take the example of the spin-2 Virasoro operator, and the computation for the other operators is similar. At the quantum level, we would like to check if [Q, T ] = 0 is consistent as an quantum operator identity. This amounts to check, for any operator O(x), the following identity where S µ is the supersymmetry current whose corresponding charge isQ and the appropriate ordering of the operators is implicitly adopted. Following [63], we can choose special operators O(x) where necessary condition for the quantum cohomology requirement could be obtained from (3.33). Explicitly, one rewrites (3.33) into the relation  (3.35) that are conjugates to the supercharges, namely Therefore theQ cohomology is the same as theD cohomology that is simpler to compute. Furthermore, as argued in [63], one can use the the free fields commutators, and the equation of motion 37) to compute the leading contribution to the OPEs.
The logic of the following computation is to choose some O(x) operator whose supersymmetry action [Q, O(x)} is known, then by computing the leading order term of the second term of (3.33) we get the value of the first term This gives Before doing the remaining OPE we first consider what would be the form of the result to the leading correction of the coupling in the potential. This is all we need for the purpose of our computation. Since the sum of the two terms of (3.34) must be zero, the leading term should be canceled with the leading terms of the first term The leading term of this expression only receives contributions from the free field part, so the only singular term in this expression must come from the OPE between λ i (x) and (the derivatives of) any linearλ i (w) term. The result will be simply a constant term due to (3.1). Therefore a necessary condition for the vanishing of the anomaly in the stress energy tensor is simply that the residual of the leading order pole of (3.40) is not proportional to the identity operator. It is easy to check that this is true for our potential J i (φ a ) since it is holomorphic in the chiral superfields while the relevant terms in T (w) contains both∂φ a and φ a . We can further consider the operator O(x) = φ a . Then it is straightforward to check that it does not lead to any new condition since the result {Q, φ a } is Grassmann odd and thus its further OPE with T (w) will not be proportional to the identity operator.
To sum up, we have shown that for the model (2.1) there could be an anomaly free stress energy tensor at the quantum level. Similar arguments apply to the other higherspin operators. Therefore we do not find extra consistency conditions at quantum level. 4 Higher spin square at the q → N M limit A property of the higher-spin algebra defined in the previous section is that it always involves multiple-sum terms that are not normal ordered products of low spin operators. The inclusion of these multiple-sum terms is required by the closure of the algebra, which, however, makes the holographic interpretation of the higher-spin operators less clear.
On the other hand, as shown in [1] explicitly, there is good evidence showing that at a special value, q → N M , one observes the emergence of a tower of conserved operators in the infrared SYK-like fix point of the model (2.1). They correspond to a tower of massless higher-spin fields in the bulk which resist the early time chaotic behavior and make the Lyapunov exponent vanish. Given the property that the chiral algebra is rigid along the RG flow, one naturally expects to see a higher-spin algebra of single-sum operators which are holographically dual to the tower of "single particle" higher-spin fields, to emerge in the chiral algebra. In addition, from the analysis of [1] one expects this to appear only at the above special limit.
So we search for higher-spin subalgebras of the chiral algebra that are generated by single-sum operators, which are analogues to the single-trace operators in matrix theories. It turns out that one can indeed find higher-spin "single-sum" subalgebras of the chiral algebra, and these subalgebras only close at the special value q → N M . Moreover, we find two, instead of one, higher-spin subalgebras. They do not commute with each other; commutators of operators from the two different subalgebra generate new operators. All the operators generated in such a manner can be organized into representations of the vertical higher-spin algebra. These properties indicate that at the special point q → N M there is a "higher-spin square" structure that was previously found in the symmetric orbifold theory [57,66]. As expected there, see also a related example [60], the higher-spin square is a stringy algebra that is closely related to the tensionless limit of the string theory. Given our expectation of the relation between the SYK-like models and some finite tension string theory, and our interpretation of the limit q → N M as a toy version of the tensionless limit, it is not surprising that we find such a "higher-spin square" structure in this limit.
In the following, we first discuss the two different higher-spin algebras. Then we discuss the higher-spin square structure they generate. At the end of this section, we comment on the relation between this algebra and the emergent higher-spin algebra discovered in [1].

The "vertical" higher-spin subalgebra
One of the higher-spin algebra at the limit q → N M consists of operators that are quadratic in the fundamental fields. We call the subalgebra generated by them the "vertical" higher-spin algebra in the notion of [57,66].
The stress energy tensor is the same as (3.2), which we recast here There is a primary spin-3 operator that is quadratic in the fundamental fields It has the following OPE where we have defined a primary spin-4 operator so that is has the normalization (4.7) Notice that we call the higher-spin algebra generated by these operators to be quadratic in the sense that the single particle operators, namely the terms with only one sum of the i or the a indices, are all quadratic in the fundamental fields. The last term, which involves the normal ordered product T T , in the W v 4 operators is added so that this operator is written in the primary basis; if we had chosen to work in the equally well-defined quasi-primary basis, this normal ordered product term can be dropped and all the terms in the operator become manifestly quadratic in the fundamental fields.
We have checked to a few more higher operators and we find only one operator at each spin. Therefore we conclude this is a conventional higher-spin type W ∞ algebra.

The "horizontal" higher-spin subalgebra
Since the algebra discussed in the last section is quadratic in terms of the fundamental fields, it is easy to check that all other operators in the chiral algebra form various representations of the vertical higher-spin algebra. In particular, single-sum operators with different number of fundamental fields are in different representations. In the following we will discuss an interesting observation that there is a higher-spin algebra structure on this set of representations. To see this we consider the higher-spin primary state of each representation of the vertical higher-spin algebra, namely the states with the lowest weight in each representation. It turns out that these lowest weight operators, which consist of single-sum terms of higher powers of the fundamental fields, generate another higher-spin algebra. We now explicitly construct the first few higher-spin operators in the horizontal algebra. The spin-2 operator is again (3.2). The primary spin-3 operator of this second higher-spin algebra is The Virasoro primary spin-4 operator reads where the repeated indices are summed over. The OPE between W h 3 (z) and itself turns out to be 3 (4.11) As the vertical algebra we checked a few more operators with higher spins and we

Higher spin square
The results in the previous subsections remind us about a similar structure, namely the "higher spin square", in 2d symmetric orbifold CFT [57,66,84]. In that context there are two different higher-spin symmetry algebras. One, which is referred to as the vertical higher-spin algebra, is generated by a tower of higher-spin generators that are quadratic in the fundamental fields. The other algebra, referred to as the horizontal higher-spin algebras, are generated by higher-spin operators that are higher powers of the fundamental fields that subject to a single sum of the repeated indices. By taking the commutators of these two algebras, a larger chiral algebra is generated and the new generators can all be organized into representation of the vertical higher-spin algebra.
In our model, we find a very similar structure at the q → N M limit. The algebra discussed in section 4.1 correspond to the generators of the vertical higher-spin algebra, and hence they are labeled by upper indices v. On the other hand, the operator discussed in section 4.2 is analogous to the horizontal higher-spin algebra, and hence they are labeled by upper indices h. The OPE of the operators in the two higher-spin algebra is non-trivial and one finds a much larger chiral algebra being generated. The latter is the realization of the stringy "higher spin square" in our model. operators in the higher spin square do close among themselves, they do not commute with the other operators of the chiral algebra that are not in the higher-spin square.
The spin-1 field is one example; it has nontrivial commutation relations with most of the operators in the higher-spin square and the commutators contains operators that are not in the higher spin square. Given that the central terms of the stress energy tensor and the spin-1 fields are both positive, and the IR fixed point of the model (2.1) is unitary. We believe there is no issue of non-unitarity; once the spin-1 fields are properly factored out, the central terms of all fields will be positive definite. We have seen partial evidence of this statement, but to verify this statement explicitly we have to properly factor out the spin-1 fields for all the operators in the higher-spin square. We defer this into a future project.
From this result, we observe that although there are no M/N dependence in the two higher-spin algebras discussed in section 4.1, 4.2, there are indeed explicitly N/M or q dependence in the higher-spin square. Since the two higher-spin algebras are fixed, both being the W ∞ [1] algebra, there is not likely further field redefinitions to remove the q dependence. Therefore q should be regarded as a new parameter that labels a family of different higher-spin square structure.

Relation with the emergent IR higher-spin symmetry
We have shown that there is not a single-sum W ∞ subalgebra in the chiral algebra of theQ-cohomology at any value of q = N M . On the other hand, in the limit q → N M , we found two different higher-spin subalgebra that further generate a higher-spin square subalgebra. The q → N M is the same limit where the IR higher-spin symmetry discussed in [1] emerges. Therefore it is natural to expect a close relation between the higherspin algebra in theQ-cohomology and the emergent IR higher-spin symmetry. But one would like to ask which higher-spin subalgebra, among those discussed in the section 4.1,4.2 and 4.3, is the "image" of this IR higher-spin algebra.
We believe that the vertical higher-spin algebra is the counterpart of the IR higherspin algebra in theQ-cohomology that can be extended away from the IR fixed point. The reason is the following. The emergent IR higher-spin symmetry is observed in the singlet channel of the 4-point function φ a φ aφb φ b and φ a φ aΛi Λ i . The antiholomorphic higher-spin operators running in this channel have the schematic form where the first two terms represent schematically the free field expressions of the higherspin operators, and the last term represents all the terms generated from the evolution by the Hamiltonian of the model. Due to the presence of the last term, we do not expect simple expressions of these operators in terms of the fundamental fields Φ a and Λ i . One the other hand, as explained originally in [62], which we briefly recast in section 2.2, the operators in theQ-cohomology and the algebraic relations among them continue to make sense even if the potential is tune to zero. In this imaginary process the operators (4.19) approaches their free field expressions, which makes it clear that their counterparts in theQ-cohomology should be those operators that are quadratic in the fundamental fields. 7 As a result, it is natural to consider the vertical higher-spin algebra found in section 4.1 to be related to the emergent IR higher-spin symmetry algebra discussed in [1].

Conclusion
In this paper we consider the chiral algebra of a class of 1+1 dimensional SYK models with N = (0, 2) and N = (2, 2) supersymmetry.
In the special limit q → N M we have constructed two different higher-spin algebras: one algebra is generated by bilinear "single-sum" operators, the other is generated by "single-sum" operators that are higher powers in the fundamental fields. The two algebras then generate a larger stringy algebra that has the structure of a "higher spin square". This gigantic stringy symmetry has been discussed previously [57,66,84]. It can be used as a guiding symmetry to organize the spectrum of certain string theories, and it is useful to clarify the relation between string theory, in some appropriate tensionless limit, and higher-spin theories. It is not surprising that a similar structure appears in the SYK model discussed in [1] and in this paper. The results in [1] indicates that as we tune a parameter, the model exhibits emergent higher-spin symmetry and mimics the transition from a tensile string theory to its tensionless limit. Therefore, the appearance of similar higher spin square structure becomes natural. As commented in section 4.4, it is natural to identify the vertical higher-spin subalgebra, which is generated by operators that are bilinear in the fundamental fields, with the emergent IR higher-spin algebra. It is thus very interesting to further clarify the explicit mapping between the two set of higher-spin-type symmetries.