1 Introduction

The large \(\mathcal{N}=4\) holography [1] connects the unitary Wolf space coset conformal field theory in two dimensions and the matrix extended higher spin theory on \(AdS_3\) space. The (large) \(\mathcal{N}=4\) supersymmetry plays an important role in this holography. One of the reasons why we want to understand the Wolf space coset construction is that this coset construction is a generalization of the free field construction. We obtain the latter by taking large level (the second order pole term of the OPE between spin-1 currents) k limit of the former. In other words, the large level k limit corresponds to the vanishing of ’t Hooft-like coupling constant [2]. In order to observe the behavior of finite ’t Hooft-like coupling constant, the Wolf space coset construction is necessary to describe its nontrivial structure fully. Then we obtain the finite (N, k) behavior in the Wolf space coset construction and this will provide some hints for the higher spin theory on \(AdS_3\) space at the quantum level. According to the results of [3], there exist different types of Wolf space cosets. One of them is given by orthogonal Wolf space coset we are interested in.Footnote 1 See also the relevant work in [7].

Contrary to the unitary Wolf space coset model [8,9,10,11], the orthogonal Wolf space coset model contains the lowest higher spin current of spin-2 rather than spin-1 [12]. This will make some calculations be rather involved. So far, the complete OPEs between the lowest 16 higher spin currents are not known. In [13], the lowest higher spin-2 current living in the lowest \(\mathcal{N}=4\) higher spin-2 multiplet (that is, there are 16 higher spin currents) in terms of orthogonal Wolf space coset fields for generic N was found. Then it is straightforward to obtain the remaining 15 higher spin currents from this higher spin-2 current by using the four supersymmetry generators of the large \(\mathcal{N}=4\) (non)linear superconformal algebra for fixed low values of N. For fixed \(N=4\), the three kinds of higher spin-3 currents were obtained from the OPEs between the higher spin currents of the above the lowest \(\mathcal{N}=4\) higher spin-2 multiplet. It was not clear how they appear in different \(\mathcal{N}=4\) multiplets at that time. We should look at the OPEs between the \(\mathcal{N}=4\) stress energy tensor and the possible \(\mathcal{N}=4\) multiplets by allowing the \(SO(\mathcal{N}=4)\) nonsinglet property to these \(\mathcal{N}=4\) multiplets.Footnote 2

In this paper, we reconsider the complete OPEs between the lowest 16 higher spin currents in the \(\mathcal{N}=4\) orthogonal Wolf space coset model. After determining the complete OPEs for the particular \(N=4\) where all the higher spin currents can be written in terms of orthogonal Wolf space coset fields, we discuss the \(N=5\) case by adding more fields. The main idea to this purpose is that by using the fundamental orthogonal Wolf space coset fields, we compute the various OPEs. When the new higher spin primary fields arise in the right hand side of the OPEs, then we should reorganize them under the \(SO(\mathcal{N}=4)\) symmetry and observe how they transform under the SO(4) symmetry.Footnote 3 For the SO(4) singlet \(\mathcal{N}=4\) multiplet in the unitary case, we can construct all the new higher spin currents living this multiplet once the lowest new higher spin current is determined with the help of the four supersymmetry generators.

However, if we have a single SO(4) nonsinglet \(\mathcal{N}=4\) multiplet which will appear in the right hand side of the our OPEs, then there are several elements on this multiplet. Each element transforms nontrivially under the 16 currents of the large \(\mathcal{N}=4\) linear superconformal algebra. The reason why we need to have this transformation is that we should calculate the OPEs between the 16 currents and the higher spin currents living in the SO(4) nonsinglet \(\mathcal{N}=4\) multiplet in order to use the Jacobi identities. More explicitly, the several lowest higher spin currents can be determined by the six spin-1 currents of the above 16 currents after computing the OPEs between them and reading off the first order poles. Once these lowest higher spin currents are obtained completely, then we can repeat the procedure for the singlet case because we can act the four supersymmetry generators on each lowest higher spin current.

The most difficult part of the present work is to write down all the possible orthogonal Wolf space composite fields appearing in the right hand sides of the OPEs in terms of the known (higher spin) currents. As the spin at the specific pole increases, the number of composite fields becomes large. When the new higher spin current appears, the situation is more involved. Even for \(N=5\), when the spin at the particular pole is large, then it is not obvious to observe how to express that pole in terms of known (higher spin) currents and a new primary higher spin current. Due to the many independent terms, it is not possible to solve the linear equations for the undetermined coefficients coming from the vanishing of the sum of the particular pole (written in terms of coset fields), the possible known composite (higher spin) current terms with arbitrary coefficients and a new higher spin current. We can calculate the OPE between the spin-\(\frac{1}{2}\) currents and the above particular pole term given in terms of coset fields. Then the first order pole of this OPE has spin which is less than the original spin by \(\frac{1}{2}\). This will make some computations easier. At the same time, we can calculate the OPE between the spin-\(\frac{1}{2}\) currents and the above sum of known composite terms and a new higher spin current. By comparing these two expressions, we can reduce the number of unknown coefficients and moreover, by using other conditions from other spin-\(1, \frac{3}{2}\) and 2 currents from the above \(\mathcal{N}=4\) primary, we can eventually determine all the coefficients if there is no new higher spin current. If not, the new higher spin current can be written in terms of the known composite terms as well as the extra terms which can be written in terms of coset fields.

Instead of considering all the possible \(136(=\sum _{i=1}^{16} i)\) OPEs (which arise from the OPEs between 16 higher spin currents), we focus on the 16 OPEs among them (the OPEs between 16 higher spin currents and the lowest higher spin-2 current or their reversed OPEs) because the remaining 120 OPEs can be extracted from the \(\mathcal{N}=4\) supersymmetry. That is, once the above 16 OPEs are obtained (this implies that we can write down the corresponding \(\mathcal{N}=4\) OPE explicitly by putting these five kinds of OPEs into the expansion of fermionic coordinates), then the \(\mathcal{N}=4\) superspace description allows us to write down the above 120 OPEs automatically by multiplying various super derivatives and putting the fermionic coordinates to zero in this OPE successively. Then after inserting the arbitrary coefficients which will depend on (N, k) explicitly (and possible other structure constants) in front of all the composite fields arising in the right hand side of 136 OPEs, the Jacobi identities can be used. Eventually, the complete OPEs can be determined and we will present them in a single OPE in \(\mathcal{N}=4\) superspace.

In Sect. 2, the \(\mathcal{N}=4\) orthogonal Wolf space coset model is reviewed.Footnote 4

In Sect. 3, based on the findings in [13] which is valid for \(N=4\), the new observations will be added.

In Sect. 4, based on the results of section 3, we consider \(N=5\) case. We will find various new higher spin currents (some of them are not present for \(N=4\) case). The 136 OPEs will be obtained eventually.

In Appendices, some of the detailed expressions described in previous sections are given.

The Thielemans package [29] is used.

An ancillary (mathematica) file \(\mathtt {ancillary.nb}\), where the complete OPEs with the explicit structure constants appearing in Appendices B and C are given, is included.

2 Review of \(\mathcal{N}=4\) orthogonal Wolf space coset model

We consider the Wolf space coset in the ‘supersymmetric’ version with groups \(G=SO(N+4)\) and \(H=SO(N) \times SO(4)\) as followsFootnote 5:

$$\begin{aligned} \hbox {Wolf}= \frac{G}{H} = \frac{SO(N+4)}{SO(N) \times SO(4)}. \end{aligned}$$
(2.1)

The group indices are denoted by

$$\begin{aligned}&G \quad \hbox {indices} : a, b, c, \ldots \nonumber \\&\quad =1,2, \ldots , \frac{1}{4} (N+4)\nonumber \\&\quad (N+3), 1^{*}, 2^{*}, \ldots , \left( \frac{1}{4} (N+4)(N+3) \right) ^{*}, \nonumber \\&\frac{G}{H} \quad \hbox {indices} : \bar{a},\bar{b},\bar{c},\cdots = 1,2, \ldots , 2N, 1^{*}, 2^{*}, \ldots , 2N^{*}.\nonumber \\ \end{aligned}$$
(2.2)

In the bosonic version, there exist 4N free fermions living in the extra SO(4N) group in the numerator of the coset at level 1.

The \(\mathcal{N}=1\) affine Kac–Moody algebra can be determined by the adjoint spin-1 current and the spin-\(\frac{1}{2}\) current of group \(G=SO(N+4)\). By adding the quadratic term in the fermions to the above spin-1 current, the operator product expansions between the ‘modified’ spin-1 current \(V^a(z)\) and the spin-\(\frac{1}{2}\) current \(Q^a(z)\) are described as

$$\begin{aligned} V^a(z) \, V^b(w)&= \frac{1}{(z-w)^2} \, k \, g^{ab} \nonumber \\&\quad -\frac{1}{(z-w)} \, f^{ab}_{\,\,\,\,\,\,c} \, V^c(w) +\cdots , \nonumber \\ Q^a(z) \, Q^b(w)&= -\frac{1}{(z-w)} \, (k+N+2) \, g^{ab} + \cdots . \end{aligned}$$
(2.3)

The level k is a positive integer. The metric can be obtained from \(g_{ab} = \frac{1}{2\, c_g} \, f_{ac}^{\,\,\,d} \, f_{bd}^{\,\,\,c}\) where \(c_g\) is the dual Coxeter number of the Lie algebra \(G=SO(N+4)\). That is, \(c_g=(N+2)\). The metric \(g_{ab}\) is given by the generators of \(SO(N+4)\) in the complex basis, \(g_{ab} = \frac{1}{2} \, \hbox {Tr}(T_a \, T_b)\). The commutation relation of generators is given by \([T_a, T_b]=f_{ab}^{\,\,\,c} \, T_c\).

For given \((N+4) \times (N+4)\) matrix, the above 4N coset indices (2.2) can be associated with the following locations with asterisk

$$\begin{aligned} \left( \begin{array}{rrrrr|rrrr} &{}&{}&{}&{}&{} {*} &{} {*} &{} {*} &{} {*} \\ &{}&{}&{}&{}&{} {*} &{} {*} &{} {*} &{} {*} \\ &{}&{} &{}&{}&{} \vdots &{} \vdots &{} \vdots &{} \vdots \\ &{}&{}&{}&{}&{} {*} &{} {*} &{} {*} &{} {*} \\ &{}&{}&{}&{}&{} {*} &{} {*} &{} {*} &{} {*} \\ \hline {*} &{} {*} &{} \cdots &{} {*} &{} {*} &{}&{}&{}&{} \\ {*} &{} {*} &{} \cdots &{} {*} &{} {*} &{}&{}&{}&{} \\ {*} &{} {*} &{} \cdots &{} {*} &{} {*} &{}&{}&{}&{} \\ {*} &{} {*} &{} \cdots &{} {*} &{} {*} &{}&{}&{}&{} \\ \end{array}\right) _{(N+4) \times (N+4)}. \end{aligned}$$
(2.4)

That is, the generators with 2N coset indices have two nonzero elements located at the above \(N \times 4\) and \(4 \times N\) off diagonal matrices in (2.4). The remaining 2N coset generators can be obtained from the above coset generators by transposing. Note that the size of two block diagonals at \(N=4\) is equal to each other.

2.1 The 11 currents of \(\mathcal{N}=4\) nonlinear superconformal algebra

The four supersymmetry currents of spin-\(\frac{3}{2}\), the six spin-1 currents of \(\hat{SU}(2)_k \times \hat{SU}(2)_N\), and the spin-2 stress energy tensor [30,31,32,33] can be described as follows:

$$\begin{aligned} \hat{G}^{0}(z)= & {} \frac{i}{(k+N+2)} \, Q_{\bar{a}} \, V^{\bar{a}}(z), \nonumber \\ \hat{G}^{i}(z)= & {} \frac{i}{(k+N+2)} \, h^{i}_{\bar{a} \bar{b}} \, Q^{\bar{a}} \, V^{\bar{b}}(z), \qquad i =1,2,3, \nonumber \\ \hat{A}_{i}(z)= & {} (-1)^{i+1} \frac{1}{4N} \, f^{\bar{a} \bar{b}}_{\,\,\,\,\,\, c} \, h^i_{\bar{a} \bar{b}} \, V^c(z), \nonumber \\ \hat{B}_{i}(z)= & {} -\frac{1}{4(k+N+2)} \, h^i_{\bar{a} \bar{b}} \, Q^{\bar{a}} \, Q^{\bar{b}}(z), \nonumber \\ \hat{T}(z)= & {} \frac{1}{2(k+N+2)^2} \left[ (k+N+2) \, V_{\bar{a}} \, V^{\bar{a}} +k \, Q_{\bar{a}} \, \partial \, Q^{\bar{a}}\right. \nonumber \\&\left. +f_{\bar{a} \bar{b} c} \, Q^{\bar{a}} \, Q^{\bar{b}} \, V^c \right] (z) \nonumber \\&- \frac{1}{(k+N+2)} \sum _{i=1}^3 ( (-1)^{i} \hat{A}_{i} + \hat{B}_{i} )^2 (z). \end{aligned}$$
(2.5)

Here the three almost complex structures are given by

$$\begin{aligned} h^1_{\bar{a} \bar{b}}= & {} \left( \begin{array}{cccc} 0 &{} 1 &{} 0 &{} 0 \\ -1 &{} 0 &{} 0 &{} 0 \\ 0 &{} 0 &{} 0 &{} 1 \\ 0 &{} 0 &{} -1 &{} 0 \\ \end{array} \right) , h^2_{\bar{a} \bar{b}} = \left( \begin{array}{cccc} 0 &{}i &{} 0 &{} 0 \\ -i &{} 0 &{} 0 &{} 0 \\ 0 &{} 0 &{} 0 &{} -i \\ 0 &{} 0 &{} i &{} 0 \\ \end{array} \right) , \nonumber \\ h^{3}_{\bar{a} \bar{b}}= & {} \left( \begin{array}{cccc} 0 &{} 0 &{} i &{} 0 \\ 0 &{} 0 &{} 0 &{} i \\ -i &{} 0 &{} 0 &{} 0 \\ 0 &{} -i &{} 0 &{} 0 \\ \end{array} \right) , \end{aligned}$$
(2.6)

where each entry is \(N \times N\) matrix. Note that we introduce \(h^{0}_{\bar{a} \bar{b}} \equiv g_{\bar{a} \bar{b}}\). The only coset indices associated with (2.4) appear in the fermionic fields.

We introduce the above 11 currents in different basis in order to describe them in the \(\mathcal{N}=4\) superspace as follows:

$$\begin{aligned} \hat{T}(z)\rightarrow & {} \hat{L}(z), \nonumber \\ \hat{G}^0(z)\rightarrow & {} \hat{G}_{ss}^2(z), \qquad \hat{ G}^1(z) \rightarrow \hat{G}_{ss}^3(z), \nonumber \\ \hat{G}^2(z)\rightarrow & {} -\hat{G}_{ss}^4(z), \qquad \hat{G}^3(z) \rightarrow \hat{G}_{ss}^1(z), \nonumber \\ \hat{A}^{\pm 1}(z)\rightarrow & {} \frac{i}{2} ( \hat{T}^{14} \mp T^{23})(z), \nonumber \\ \hat{A}^{\pm 2}(z)\rightarrow & {} \frac{i}{2} ( \hat{T}^{13} \pm \hat{T}^{24})(z), \,\, \hat{A}^{\pm 3}(z) \rightarrow \pm \frac{i}{2} ( \hat{T}^{12} \mp \hat{T}^{34})(z).\nonumber \\ \end{aligned}$$
(2.7)

We use \(\hat{A}_1(z) =-\hat{A}^{+1}(z)\), \( \hat{A}_2(z) =\hat{A}^{+2}(z)\), \(\hat{A}_3(z) =-\hat{A}^{+3}(z)\) and \(\hat{B}_i(z) = \hat{A}^{-i}(z)\). Then the 11 currents in this new basis satisfy Appendix (A.1) of [34] by using (2.3).

2.2 The 16 currents of \(\mathcal{N}=4\) linear superconformal algebra

In the linear version [3, 35,36,37,38] of \(\mathcal{N}=4\) superconformal algebra, there are also four spin-\(\frac{1}{2}\) currents and one spin-1 current

$$\begin{aligned} { F}_{11}(z)= & {} \frac{i }{\sqrt{2}} Q^{(2N+3)}(z), \qquad { F}_{22}(z)=-\frac{i }{\sqrt{2}} Q^{(2N+3)^*}(z), \nonumber \\ { F}_{12}(z)= & {} \frac{(1-i) }{2} Q^{(2N+2)^*}(z), \qquad \nonumber \\{ F}_{21}(z)= & {} \frac{(1+i) }{2} Q^{(2N+2)}(z), \nonumber \\ { U}(z)= & {} \frac{(1+i)}{2 \sqrt{2}} V^{(2N+2)} (z)+ \frac{(-1+i)}{2 \sqrt{2}} V^{(2N+2)^*}(z) \nonumber \\&+\frac{i}{(N+k+2)} Q^{(2N+1)} Q^{(2N+1)^*}(z) \nonumber \\&- \frac{i}{2(N+k+2)} \left( \sum _{a=1}^N Q^{a} Q^{a^*} - \sum _{a=N+1}^{2N} Q^{a} Q^{a^*} \right) (z).\nonumber \\ \end{aligned}$$
(2.8)

We use the following transformation in order to describe in the \(\mathcal{N}=4\) superspace

$$\begin{aligned}&\frac{i}{\sqrt{2}}(F_{12}+F_{21})(z) \rightarrow \Gamma ^{1}(z), \nonumber \\&\quad -\frac{1}{\sqrt{2}}(F_{12}+F_{21})(z) \rightarrow \Gamma ^{2}(z), \nonumber \\&-\frac{i}{\sqrt{2}}(F_{11}+F_{22})(z) \rightarrow \Gamma ^{3}(z), \nonumber \\&\quad -\frac{1}{\sqrt{2}}(F_{11}-F_{22})(z) \rightarrow \Gamma ^{4}(z), \qquad U(z) \rightarrow U(z). \nonumber \\ \end{aligned}$$
(2.9)

Furthermore, the six spin-1 currents of \(\hat{SU}(2)_{k+1} \times \hat{SU}(2)_{N+1}\), four supersymmetry currents of spin-\(\frac{3}{2}\) and the spin-2 stress energy tensor can be described as follows [30]:

$$\begin{aligned} T^{ij}(z)= & {} \hat{T}^{ij}(z) +\frac{2\,i}{(2+k+N)}\,\Gamma ^{i}\Gamma ^{j}(z), \nonumber \\ G^{i}(z)= & {} \hat{G}_{ss}^{i}(z)+\frac{2\,i}{(2+k+N)}\, U\Gamma ^{i}(z)\nonumber \\&+\varepsilon _{ijkl}\, \Bigg [\, \frac{4\,i}{3(2+k+N)^{2}}\,\Gamma ^{j}\Gamma ^{k}\Gamma ^{l}\nonumber \\&-\frac{1}{(2+k+N)}\,T^{jk}\Gamma ^{l} \,\Bigg ](z), \nonumber \\ L(z)= & {} \hat{L}(z)- \frac{1}{(2+k+N)}\Big (\,UU-\partial \Gamma ^{i} \Gamma ^{i}\,\Big )(z). \end{aligned}$$
(2.10)

Here we should use (2.5), (2.6), (2.7) and (2.8). The \(\mathcal{N}=4\) linear superconformal algebra can be realized from (2.9) and (2.10) in the orthogonal Wolf space coset model (2.1). The 16 currents are given by L(z), \(G^i(z)\), \(T^{ij}(z)\), \(\Gamma ^i(z)\) and U(z) which can be written in terms of \(Q^{\bar{a}}(z)\) and \(V^a(z)\).

2.3 The 16 lowest higher spin currents

The lowest higher spin-2 current in the \(\mathcal{N}=4\) orthogonal Wolf space is found in [13] and it is given by

$$\begin{aligned} \Phi _0^{(2)}(z)= & {} c_1 \, V_{\bar{a}} V^{\bar{a}}(z)+ c_2 \, \sum _{a':so(N)} V_{a'} V^{a'}(z)\nonumber \\&+ c_3 \, \sum _{a'':so(4)} V_{a''} V^{a''}(z)+ c_4 \, \sum _{i=1}^3 \hat{A}_i \hat{A}_i(z) \nonumber \\&+ c_5 \, \sum _{i=1}^3 \hat{B}_i \hat{B}_i(z) + c_6 \, Q_{\bar{a}} \partial Q^{\bar{a}}(z) \nonumber \\&+ c_7 \, \sum _{\mu =0}^3 h^{\mu }_{\bar{a} \bar{b}} h^{\mu }_{\bar{c} \bar{d}} f^{\bar{a} \bar{c}}_{\,\,\,\,\,\, e} Q^{\bar{b}} Q^{\bar{d}} V^e(z), \end{aligned}$$
(2.11)

where the coefficients areFootnote 6

$$\begin{aligned} c_1= & {} -\frac{(2 k^2 N+k^2+4 k N^2+6 k N+2 k+11 N^2-2 N-24)}{2 (k-1) N (k+N+2)^2}, \nonumber \\ c_2= & {} \frac{6 (2 k N+3 k+3 N+4)}{(k-1) N (k+N+2)^2}, \nonumber \\ c_3= & {} \frac{3 (k+N-2) (2 k N+3 k+3 N+4)}{2 (k-1) (k+2) (k+N+2)^2}, \nonumber \\ c_4= & {} \frac{2 (N+2) (2 k+N)}{(k+2) (k+N+2)^2}, \qquad c_5 = \frac{2 k (2 k+N)}{N (k+N+2)^2}, \nonumber \\ c_6= & {} \frac{k (N+2) (2 k+N)}{N (k+N+2)^3}, \qquad c_7 = \frac{(N+2) (2 k+N)}{4 N (k+N+2)^3}. \end{aligned}$$
(2.12)

The OPE between this higher spin-2 current and itself is described as

$$\begin{aligned} \Phi _{0}^{(2)}(z)\,\Phi _{0}^{(2)}(w)= & {} \frac{1}{(z-w)^{4}}\,\, c^{0} +\frac{1}{(z-w)^{2}}\, Q^{(2)}_{0}(w)\nonumber \\&+\frac{1}{(z-w)}\, \frac{1}{2}\,\partial Q^{(2)}_{0}(w) +\cdots , \end{aligned}$$
(2.13)

where the central term is

$$\begin{aligned} c_{0}^{0,4} = \frac{3k(2k+N)(2kN+3k+3N+4)(2k^2N+k^2+4kN^2+6kN+2k+11N^2-2N-24)}{(k-1)(k+2)N(k+N+2)^3}. \end{aligned}$$
(2.14)

The quasi primary field \(Q^{(2)}_{0}(w)\) depends on the higher spin-2 current itself \(\Phi ^{(2)}_{0}(w)\) and the spin-\(\frac{1}{2}, 1, 2\) currents of the \(\mathcal{N}=4\) linear superconformal algebra. The explicit form is given Appendix B. The 16 higher spin currents can be combined into one single \(\mathcal{N}=4\) super field as follows:

$$\begin{aligned}&\varvec{\Phi }^{(s=2)}(Z) \nonumber \\&\quad \equiv \Bigg ( \Phi _{0}^{(2)}(z), \Phi _{\frac{1}{2}}^{(2),i}(z), \Phi _{1}^{(2),ij}(z), \Phi _{\frac{3}{2}}^{(2),i}(z), \Phi _{2}^{(2)}(z) \Bigg ).\nonumber \\ \end{aligned}$$
(2.15)

In this paper, we construct 136 OPEs between the higher spin currents in (2.15) explicitly.Footnote 7

3 The OPE for \(N=4\)

In this section, we continue to calculate the OPEs between the lowest higher spin-2 multiplet for \(N=4\). Some of the OPEs were found in [13]. We would like to obtain the general structure of these OPEs which will give us some hints for the general N in next section.

3.1 The known facts

For fixed \(N=4\), it was straightforward to calculate the various higher spin currents in \(\mathcal{N}=4\) orthogonal Wolf space starting from the above higher spin-2 current in (2.11). One of the main results in [13] was to obtain the new three higher spin-3 currents which live in different higher spin multiplet. It was not clear how they appear in the right hand side of the whole 136 OPEs. They will turn out to be the lowest components of three SO(4) vector \(\mathcal{N}=4\) higher spin-3 multiplets in next subsection.

3.2 The complete OPEs in components and \(\mathcal{N}=4\) superspace

In order to observe the symmetry behind the presence of the new three higher spin-3 currents, we should go into the \(\mathcal{N}=4\) superspace. It is known that the above 16 currents can be combined into the following single \(\mathcal{N}=4\) super field [37]

$$\begin{aligned} \mathbf{J}(Z)= & {} - \Delta (z) + i \theta ^{j} \Gamma ^{j}(z)- \frac{i}{2}\theta ^{4-jk} T^{jk}(z)\nonumber \\&-\theta ^{4-j} (G^{j}-2 \alpha i \partial \Gamma ^{j})(z)+\theta ^{4-0} (2L-2\alpha \partial ^2 \Delta )(z) \nonumber \\\equiv & {} - \Delta (z) + i \theta ^{j} \Gamma ^{j}(z)- \frac{i}{2}\theta ^{4-jk} T^{jk}(z)\nonumber \\&-\theta ^{4-j} \, \tilde{G}^{j}(z)+ \theta ^{4-0} \, 2 \, \tilde{L}(z). \end{aligned}$$
(3.1)

Here we use the notation \(\theta ^{4-0}\) for the product of fermionic coordinates \(\theta ^{4-0} \equiv \theta ^1 \, \theta ^2 \, \theta ^3 \, \theta ^4\) and we have \( U(z) \equiv -\partial \, \Delta (z)\). The parameter \(\alpha \) appears in the above and is given by

$$\begin{aligned} { \alpha }\equiv & {} \frac{1}{2}\frac{(k^{+}-k^{-})}{(k^{+}+k^{-})}, \qquad k^+ \equiv k+1, \qquad k^- \equiv N+1. \end{aligned}$$
(3.2)

Then the explicit realization for these 16 currents described in previous section can be inserted into the above single \(\mathcal{N}=4\) SO(4) singlet super field.

It is known that the \(\mathcal{N}=4\) higher spin multiplet, which transforms nontrivially under the SO(4) [the index \(\alpha \) stands for this representation which is nothing to do with (3.2)], of (conformal) (super)spin s has the following OPE with the above \(\mathcal{N}=4\) stress energy tensor as follows [37]:

$$\begin{aligned} \mathbf{J }(Z_{1}) \, \varvec{\Phi }^{(s),\alpha }(Z_{2})= & {} \frac{\theta _{12}^{4-0}}{z_{12}^{2}} \, 2s\, \varvec{\Phi }^{(s),\alpha }(Z_{2})+ \frac{\theta _{12}^{4-i}}{z_{12}} \, D^{i} \varvec{\Phi }^{(s),\alpha }(Z_{2})\nonumber \\&+ \frac{\theta _{12}^{4-0}}{z_{12}} \, 2 \, \partial \varvec{\Phi }^{(s),\alpha }(Z_{2}) \nonumber \\&- \frac{i}{2} \,\frac{\theta _{12}^{4-ij}}{z_{12}}({ T}^{ij})^{\alpha \beta }\, \varvec{\Phi }^{(s),\beta }(Z_{2}) +\cdots .\nonumber \\ \end{aligned}$$
(3.3)

Note that for the SO(4) singlet higher spin multiplet the last term in (3.3) will disappear. We will see two kinds of \(T^{ij}\) which span the representation of the SO(4) Lie algebra in this OPE of this paper.

By using the 16 component fields for fixed \(\alpha \) as in

$$\begin{aligned}&\varvec{\Phi }^{(s),\alpha }(Z) \nonumber \\&\quad \equiv \Bigg ( \Phi _{0}^{(s),\alpha }(z), \Phi _{\frac{1}{2}}^{(s),i,\alpha }(z), \Phi _{1}^{(s),ij,\alpha }(z), \Phi _{\frac{3}{2}}^{(s),i,\alpha }(z), \Phi _{2}^{(s),\alpha }(z) \Bigg ), \!\!\!\nonumber \\ \end{aligned}$$
(3.4)

the various complicated component results of (3.3) are presented in Appendix A.

We can show that the three higher spin-3 currents (\(P^{(3)}(z)\), \(\tilde{Q}^{(3)}_{-}(z)\) and \(\tilde{R}^{(3)}_{+}(z)\)) found in [13] with proper change of basis can be written in terms of the three lowest components of the three SO(3) vector \(\varvec{\Phi }^{(s=3),\alpha }(Z)\) where \(\alpha =1,2,3\) and this higher spin-3 multiplet transforms as in (3.3) with

$$\begin{aligned} { T}^{12}= & {} \left( \begin{array}{cccc} 0 &{} i &{} 0 \\ -i &{} 0 &{} 0 \\ 0 &{} 0 &{} 0 \end{array}\right) = -T^{34}, \qquad { T}^{13}=\left( \begin{array}{cccc} 0 &{} 0 &{} i \\ 0 &{} 0 &{} 0 \\ -i &{} 0 &{} 0 \end{array}\right) = -T^{24}, \nonumber \\ { T}^{23}= & {} \left( \begin{array}{cccc} 0 &{} 0 &{} 0 \\ 0 &{} 0 &{} -i \\ 0 &{} i &{} 0 \end{array}\right) = -T^{14}. \end{aligned}$$
(3.5)

We can check that the above three SO(3) (inside of SO(4)) generators in (3.5) satisfy \([T^i,T^j]= \varepsilon ^{ijk} \, T^k\) where \(T^i \equiv -\frac{i}{2} \, \varepsilon ^{ijk}\, T^{jk}\).Footnote 8 Each higher spin-3 multiplet \(\varvec{\Phi }^{(s=3),\alpha }(Z)\) has 16 components of higher spin currents according to (3.4). Totally, we have 48 higher spin currents.

Furthermore, there exist the SO(4) singlet higher spin-4 multiplet \(\varvec{\Phi }^{(s=4)}(Z)\) and the four SO(4) vector higher spin-\(\frac{9}{2}\) multiplets \(\varvec{\Phi }^{(s=\frac{9}{2}),i}(Z)\) with \(i=1,2,3,4\). They have their component fields as in (3.4).

It turns out that there will be a problem to generalize the OPEs for general N without introducing the new primary fields which will be discussed in next section.

4 The OPE for \(N \ge 5\)

We would like to construct the OPEs between the 16 higher spin currents for generic N in component approach and in \(\mathcal{N}=4\) superspace.

4.1 What happens for \(N=5\)?

It is natural to take the OPEs found in previous section and introduce the arbitrary coefficients in front of the composite fields appearing in the right hand side of the OPEs. It is straightforward to apply the Jacobi identity in order to determine these coefficients completely. It turns out that there are no consistent solutions unless we introduce the new primary fields. Therefore, we focus on the case of \(N=5\) in order to understand the algebraic structures more clearly.

4.2 The new higher spin current of spin 2

Let us consider the OPE between the last component and the first component of (2.15). How we can obtain the last component in terms of orthogonal Wolf space coset fields from the first component? According to the OPEs in Appendix A with a singlet \(\alpha \), we obtain the second component using the OPE between \(G^i(z)\) and the first component \(\Phi _0^{(2)}(w)\). Then we can calculate the OPE between \(G^i(z)\) and the second component \(\Phi _{\frac{1}{2}}^{(2),j}(w)\) with \(i \ne j\). Then the third component can be determined. The fourth component can be obtained by the OPE between \(G^i(z)\) and the third component \(\Phi _{1}^{(2),jk}(w)\) with \(i=k\). Finally, the last component can be determined by the OPE between \(G^i(z)\) and the fourth component \(\Phi _{\frac{3}{2}}^{(2),j}(w)\) with \(i=j\).

The reason for describing this particular OPE rather than others is that the structure of the right hand side of this OPE will be simple because these two higher spin currents are SO(4) singlets rather than nonsinglets, although we have found the presence of this new higher spin-2 current in other OPEs.

Let us emphasize that the last component of (2.15) is not a quasi primary field. See also Appendix A. As in unitary case [8], we subtract the additional terms from the last component of (2.15) in order to make it to be primary field.Footnote 9 It turns out that we have the following OPE with implicit notation

$$\begin{aligned}&\Phi _{2}^{(2)}(z)\,\Phi _{0}^{(2)}(w)\nonumber \\&\quad = \frac{1}{(z-w)^{6}}\, 8 \, \alpha \, c_{0}^{0,4}\nonumber \\&\qquad +\frac{1}{(z-w)^{5}}\, Q^{(1)}_{2}(w) +\frac{1}{(z-w)^{4}}\,\Bigg [\, \frac{3}{2}\,\partial Q^{(1)}_{2} +Q^{(2)}_{2} \,\Bigg ](w) \nonumber \\&\qquad + \frac{1}{(z-w)^{3}}\,\Bigg [\, \partial ^{2} Q^{(1)}_{2} +\partial Q^{(2)}_{2} +Q^{(3)}_{2} \,\Bigg ](w) \nonumber \\&\qquad + \frac{1}{(z-w)^{2}}\,\Bigg [\, \frac{5}{12}\,\partial ^{3} Q^{(1)}_{2} +\frac{1}{2}\, \partial ^{2} Q^{(2)}_{2} +\frac{5}{6}\,\partial Q^{(3)}_{2} +Q^{(4)}_{2} \,\Bigg ](w) \nonumber \\&\qquad + \frac{1}{(z-w)}\,\Bigg [\, \frac{1}{8}\,\partial ^{4} Q^{(1)}_{2} +\frac{1}{6}\,\partial ^{3} Q^{(2)}_{2} +\frac{5}{14}\,\partial ^{2} Q^{(3)}_{2}\nonumber \\&\qquad +\frac{3}{4}\,\partial Q^{(4)}_{2} + Q^{(5)}_{2} \,\Bigg ](w) \nonumber \\&\qquad + p_{1}\,\sum \limits _{n=3}^{4} \frac{1}{(z-w)^{n}} \, \Big \{ \partial ^{2} \Phi _{0}^{(2)} \,\Phi _{0}^{(2)} \Big \}_{n}(w)\nonumber \\&\qquad + p_{2}\,\sum \limits _{m=2}^{4} \frac{1}{(z-w)^{m}}\, \Big \{ (L \Phi _{0}^{(2)})\,\Phi _{0}^{(2)} \Big \}_{m}(w) \nonumber \\&\qquad + \cdots . \end{aligned}$$
(4.1)

The central term is proportional to the previous central term in (2.14) together with (3.2).Footnote 10 Note that the above central term comes from these \(p_1\) and \(p_2\) terms.

We will use the following quasi primary fields with their spins, SO(4) indices i, j and the subscript indicating the number of fermionic coordinates

$$\begin{aligned}&Q_0^{(2)}(z) \, ; \, Q_{\frac{1}{2}}^{\left( \frac{3}{2}\right) ,i}(z), \, Q_{\frac{1}{2}}^{\left( \frac{5}{2}\right) ,i}(z), \nonumber \\&Q_{\frac{1}{2}}^{\left( \frac{7}{2}\right) ,i}(z) \, ; \, Q_{1}^{(1),ij}(z), \, Q_{1}^{(2),ij}(z), \, Q_{1}^{(3),ij}(z), \, Q_{1}^{(4),ij}(z) \, ; \nonumber \\&Q_{\frac{3}{2}}^{\left( \frac{1}{2}\right) ,ij}(z), \, Q_{\frac{3}{2}}^{\left( \frac{3}{2}\right) ,ij}(z), \, Q_{\frac{3}{2}}^{\left( \frac{5}{2}\right) ,ij}(z), \nonumber \\&Q_{\frac{3}{2}}^{\left( \frac{7}{2}\right) ,ij}(z), \, Q_{\frac{3}{2}}^{\left( \frac{9}{2}\right) ,ij}(z) \, ; \, Q_2^{(1)}(z), \nonumber \\&Q_2^{(2)}(z), \, Q_2^{(3)}(z), \, Q_2^{(4)}(z), \, Q_2^{(5)}(z). \end{aligned}$$
(4.2)

Note that the spin is given by the number inside the bracket and we do not add the subscript for the spin, contrary to the notation of (3.4).Footnote 11

First of all, the fifth order pole has the spin-1 current U(w) of the \(\mathcal{N}=4\) linear superconformal algebra from Appendix B. The next fourth order pole contains the descendant field \(\partial U(w)\) with the known coefficient and other terms.

We observe that there exists a new primary higher spin field of spin-2 denoted by \(X_0^{(2)}(w)\) which cannot be written in terms of the known composite fields of the currents and higher spin currents as in Appendix B. That is,

$$\begin{aligned} Q_2^{(2)}(w)= & {} w_{1,2} \, \Phi _0^{(2)}(w)+ w_{2,2} \, X_0^{(2)}(w) + \cdots , \end{aligned}$$
(4.3)

where other remaining terms are given in Appendix B. See also Appendix E for explicit form for the \(X_0^{(2)}(w)\) for \(N=5\). By considering the condition that the fourth order pole of the OPE between \(\Phi _0^{(2)}(z)\) and \(X_0^{(2)}(w)\) should vanish, we can determine the structure constant appearing the \(\Phi _0^{(2)}(w)\) of (4.3). For \(N=5\), we obtain this particular structure constant as

$$\begin{aligned}&w_{1,2}\equiv C_{(4)(2)}^{(2)}\Bigg |_{N=5} \nonumber \\&\ = -\frac{36}{5 (k-1) (k+2) (k+7)^2 (28 k+61) (11 k^2+132 k+241)} \nonumber \\&\qquad \times (7227 k^7+201718 k^6+2017067 k^5\nonumber \\&\qquad +8606534 k^4+ 13128257 k^3 \nonumber \\&\qquad - 11460814 k^2-54096247 k-42478238). \end{aligned}$$
(4.4)

At the moment, it is rather difficult to determine the N generalization of (4.4) because although we can expect the N dependence for the denominator of (4.4) by increasing the N values, the numerical values appearing in front of kth power in the numerator are functions of N. Even if we can try to calculate (4.4) for seven (which is the maximum power of k) N values where \(N=5,8,9,12,13,16\) and 17, it will take too much time to extract the higher spin-4 \(\Phi _2^{(2)}(z)\). In this paper, the above structure constant for generic N is not determined.

The next third order pole can be expressed in terms of the descendant fields and other known composite fields where there are two higher spin dependent terms \(U \,\Phi _0^{(2)}(w)\) and \(\Gamma ^i \, \Phi _{\frac{1}{2}}^{(2), i}(w)\). The four component fields of spin-\(\frac{5}{2}\) in (2.15) arise in this pole.

The second and first order poles appearing in the third and fourth lines of (4.1) will be described later subsection. We will observe that there will be a primary higher spin-4 current.

We can easily see that the singular terms appearing in the last two terms of (4.1) can be rewritten as follows.Footnote 12

Moreover, for the second term, by introducing

$$\begin{aligned} \Big \{ (L \Phi _{0}^{(2)})\,\Phi _{0}^{(2)} \Big \}_{n+2} \equiv E^{(4-n)}_{2}, \quad n=0,1,2, \end{aligned}$$
(4.5)

the following relations for three in (4.5) can be obtained from the OPE \((L \Phi _{0}^{(2)})(z) \,\Phi _{0}^{(2)}(w)\), where the previous relation (2.13) is used, in terms of \(Q_0^{(2)}(w)\), L(w) and \(\Phi _0^{(2)}(w)\),

$$\begin{aligned} E^{(2)}_{2}(w)= & {} (\, 4 Q^{(2)}_{0}+c^{0,4}_{0}L \,)(w), \nonumber \\ E^{(3)}_{2}(w)= & {} \left( \, \frac{5}{2} \, \partial \, Q^{(2)}_{0}+c^{0,4}_{0} \, \partial \, L \,\right) (w), \nonumber \\ E^{(4)}_{2}(w)= & {} \left( \, \frac{1}{2}\,\partial ^{2} Q^{(2)}_{0} + L Q^{(2)}_{0} +\frac{1}{2}\,\partial ^{2} L +2\,\Phi ^{(2)}_{0}\Phi ^{(2)}_{0} \,\right) (w), \nonumber \\ \end{aligned}$$
(4.6)

where \(c^{0,4}_{0}\) is the central term of the OPE between the lowest higher spin-2 current in (2.14).

Therefore, we can present the above OPE (4.1), together with (4.6), in complete form as follows:

$$\begin{aligned}&\Phi _{2}^{(2)}(z)\,\Phi _{0}^{(2)}(w) \nonumber \\&\quad = \frac{1}{(z-w)^{6}}\, 8 \, \alpha \, c_{0}^{0,4} +\frac{1}{(z-w)^{5}}\, Q^{(1)}_{2}(w) \nonumber \\&\qquad + \frac{1}{(z-w)^{4}}\,\Bigg [\, \frac{3}{2}\,\partial Q^{(1)}_{2} +Q^{(2)}_{2}-6\,p_{1}\,Q^{(2)}_{0} -p_2 \, E^{(2)}_{2} \,\Bigg ](w) \nonumber \\&\qquad + \frac{1}{(z-w)^{3}}\,\Bigg [\, \partial ^{2} Q^{(1)}_{2} +\partial Q^{(2)}_{2} +Q^{(3)}_{2}\nonumber \\&\qquad -p_1 \, \partial \, Q_0^{(2)} - p_2 \, E^{(3)}_{2} \,\Bigg ](w) \nonumber \\&\qquad + \frac{1}{(z-w)^{2}}\,\Bigg [\, \frac{5}{12}\,\partial ^{3} Q^{(1)}_{2} +\frac{1}{2}\, \partial ^{2} Q^{(2)}_{2}\nonumber \\&\qquad +\frac{5}{6}\,\partial Q^{(3)}_{2} +Q^{(4)}_{2} -p_2 \, E^{(4)}_{2} \,\Bigg ](w) \nonumber \\&\qquad + \frac{1}{(z-w)}\,\Bigg [\, \frac{1}{8}\,\partial ^{4} Q^{(1)}_{2} +\frac{1}{6}\,\partial ^{3} Q^{(2)}_{2} +\frac{5}{14}\,\partial ^{2} Q^{(3)}_{2}\nonumber \\&\qquad +\frac{3}{4}\,\partial Q^{(4)}_{2} + Q^{(5)}_{2} \,\Bigg ](w) \nonumber \\&\qquad + \cdots , \end{aligned}$$
(4.7)

where the relations (3.2), (2.14), (4.6) and (4.13) are used. All the singlet quasi primary fields appearing in this OPE are obtained and we present the partial expressions given in Appendix B.

4.3 The new higher spin currents of spin 3

Let us describe the OPE between the higher spin-3 (primary) current transforming as the SO(4) adjoint representation and the higher spin-2 current. We observe the following OPE

$$\begin{aligned}&\Phi _{1}^{(2),ij}(z)\,\Phi _{0}^{(2)}(w)\nonumber \\&\quad = \frac{1}{(z-w)^{4}}\, Q^{(1),ij}_{1}(w) +\frac{1}{(z-w)^{3}}\,\Bigg [\, \partial Q^{(1),ij}_{1}+ Q^{(2),ij}_{1} \,\Bigg ](w) \nonumber \\&\qquad + \frac{1}{(z-w)^{2}}\,\Bigg [\, \frac{1}{2}\,\partial ^{2} Q^{(1),ij}_{1} +\frac{3}{4}\,\partial Q^{(2),ij}_{1} + Q^{(3),ij}_{1} \,\Bigg ](w) \nonumber \\&\qquad + \frac{1}{(z-w)}\,\Bigg [\, \frac{1}{6}\,\partial ^{3} Q^{(1),ij}_{1} +\frac{3}{10}\,\partial ^{2} Q^{(2),ij}_{1}\nonumber \\&\qquad +\frac{2}{3}\,\partial Q^{(3),ij}_{1} + Q^{(4),ij}_{1}\, \Bigg ](w) \nonumber \\&\qquad - \sum \limits _{n=1}^{4}\frac{1}{(z-w)^{n}} \, (i\leftrightarrow j) +\cdots . \end{aligned}$$
(4.8)

Note that the higher spin-3 currents \(\Phi _{1}^{(2),ij}(z)\) are antisymmetric under the interchange of the index i and the index j. The last line of (4.8) implies that we should take the three lines with \( i \leftrightarrow j\) with minus sign. The quasi primary fields appearing in the fourth and third order poles, which are written in terms of the known composite fields are given in Appendix B.

After subtracting the descendant fields, the second order pole contains the six SO(4) adjoint higher spin-3 currents \(\Phi _0^{(3),\alpha }(w)\) with adjoint \(\alpha \) (that is, there will be 96 higher spin currents in these six higher spin multiplets \(\varvec{\Phi }^{(3),\alpha }(Z_2)\)).Footnote 13

By introducing the following six generators \(M^{\alpha }\) which are \(4 \times 4\) matrices with \(\alpha =1, 2, \ldots , 6\) of SO(4)

$$\begin{aligned}&{ M^{1}} \equiv L_{1}= \left( \begin{array}{cccc} 0 &{} 0 &{} 0 &{} 0\\ 0 &{} 0 &{} -i &{} 0\\ 0 &{} i &{} 0 &{} 0\\ 0 &{} 0 &{} 0 &{} 0 \end{array}\right) , { M^{2}} \equiv L_{2}= \left( \begin{array}{cccc} 0 &{} 0 &{} i &{} 0\\ 0 &{} 0 &{} 0 &{} 0\\ -i &{} 0 &{} 0 &{} 0\\ 0 &{} 0 &{} 0 &{} 0 \end{array}\right) ,\nonumber \\&{ M^{3}} \equiv L_{3}= \left( \begin{array}{cccc} 0 &{} -i &{} 0 &{} 0\\ i &{} 0 &{} 0 &{} 0\\ 0 &{} 0 &{} 0 &{} 0\\ 0 &{} 0 &{} 0 &{} 0 \end{array}\right) , \nonumber \\&{ M^{4}} \equiv K_{1}= \left( \begin{array}{cccc} 0 &{} 0 &{} 0 &{} -i\\ 0 &{} 0 &{} 0 &{} 0\\ 0 &{} 0 &{} 0 &{} 0\\ i &{} 0 &{} 0 &{} 0 \end{array}\right) ,\quad { M^{5}} \equiv K_{2} =\left( \begin{array}{cccc} 0 &{} 0 &{} 0 &{} 0\\ 0 &{} 0 &{} 0 &{} -i\\ 0 &{} 0 &{} 0 &{} 0\\ 0 &{} i &{} 0 &{} 0 \end{array}\right) ,\nonumber \\&{ M^{6}} \equiv K_{3} =\left( \begin{array}{cccc} 0 &{} 0 &{} 0 &{} 0\\ 0 &{} 0 &{} 0 &{} 0\\ 0 &{} 0 &{} 0 &{} -i\\ 0 &{} 0 &{} i &{} 0 \end{array}\right) , \nonumber \\ \end{aligned}$$
(4.9)

we can write down

$$\begin{aligned} Q_1^{(3),ij}(w) = w_{0,3} \, (M^{\alpha })^{ij} \, \Phi _0^{(3),\alpha }(w) + \cdots , \end{aligned}$$
(4.10)

as the one in Appendix B.Footnote 14 There are commutation relations \([L_i, L_j] = i \, \varepsilon _{ijk} \, L_k\), \([K_i, K_j] = i \, \varepsilon _{ijk} \, L_k\), and \([L_i, K_j] = i \, \varepsilon _{ijk} \, K_k\) [40]. The relative coefficients of (4.10) can be determined by using the defining nontrivial OPEs of \(L(z) \, \Phi _0^{(3),\alpha }(w)\), \(G^i(z) \, \Phi _0^{(3),\alpha }(w)\) and \(T^{ij}(z) \, \Phi _0^{(3),\alpha }(w)\) in Appendix A.

By recalling that from the relation (3.3) or Appendix A, the OPE \(T^{ij}(z) \, \Phi _0^{(3),\alpha }(w)\) between the spin-1 currents of the \(\mathcal{N}=4\) linear superconformal algebra and the lowest higher spin-3 currents contains the nontrivial singular terms \((T^{ij})^{\alpha \beta } \, \Phi _0^{(3),\beta }(w)\) where the generators \(T^{ij}\) in the SO(4) adjoint representation are given by

$$\begin{aligned} { T}^{12}= & {} \left( \begin{array}{cccccc} 0 &{} -i &{} 0 &{} 0 &{} 0 &{} 0 \\ i &{} 0 &{} 0 &{} 0 &{} 0 &{} 0 \\ 0 &{} 0 &{} 0 &{} 0 &{} 0 &{} 0 \\ 0 &{} 0 &{} 0 &{} 0 &{} -i &{} 0 \\ 0 &{} 0 &{} 0 &{} i &{} 0 &{} 0 \\ 0 &{} 0 &{} 0 &{} 0 &{} 0 &{} 0 \end{array}\right) , \qquad { T}^{13}=\left( \begin{array}{cccccc} 0 &{} 0 &{} -i &{} 0 &{} 0 &{} 0 \\ 0 &{} 0 &{} 0 &{} 0 &{} 0 &{} 0 \\ i &{} 0 &{} 0 &{} 0 &{} 0 &{} 0 \\ 0 &{} 0 &{} 0 &{} 0 &{} 0 &{} -i \\ 0 &{} 0 &{} 0 &{} 0 &{} 0 &{} 0 \\ 0 &{} 0 &{} 0 &{} i &{} 0 &{} 0 \end{array}\right) , \nonumber \\ { T}^{14}= & {} \left( \begin{array}{cccccc} 0 &{} 0 &{} 0 &{} 0 &{} 0 &{} 0 \\ 0 &{} 0 &{} 0 &{} 0 &{} 0 &{} -i \\ 0 &{} 0 &{} 0 &{} 0 &{} i &{} 0 \\ 0 &{} 0 &{} 0 &{} 0 &{} 0 &{} 0 \\ 0 &{} 0 &{} -i &{} 0 &{} 0 &{} 0 \\ 0 &{} i &{} 0 &{} 0 &{} 0 &{} 0 \end{array}\right) , \qquad { T}^{23}=\left( \begin{array}{cccccc} 0 &{} 0 &{} 0 &{} 0 &{} 0 &{} 0 \\ 0 &{} 0 &{} -i &{} 0 &{} 0 &{} 0 \\ 0 &{} i &{} 0 &{} 0 &{} 0 &{} 0 \\ 0 &{} 0 &{} 0 &{} 0 &{} 0 &{} 0 \\ 0 &{} 0 &{} 0 &{} 0 &{} 0 &{} -i \\ 0 &{} 0 &{} 0 &{} 0 &{} i &{} 0 \end{array}\right) , \nonumber \\ { T}^{24}= & {} \left( \begin{array}{cccccc} 0 &{} 0 &{} 0 &{} 0 &{} 0 &{} i \\ 0 &{} 0 &{} 0 &{} 0 &{} 0 &{} 0 \\ 0 &{} 0 &{} 0 &{} -i &{} 0 &{} 0 \\ 0 &{} 0 &{} i &{} 0 &{} 0 &{} 0 \\ 0 &{} 0 &{} 0 &{} 0 &{} 0 &{} 0 \\ -i &{} 0 &{} 0 &{} 0 &{} 0 &{} 0 \end{array}\right) , \qquad { T}^{34}=\left( \begin{array}{cccccc} 0 &{} 0 &{} 0 &{} 0 &{} -i &{} 0 \\ 0 &{} 0 &{} 0 &{} i &{} 0 &{} 0 \\ 0 &{} 0 &{} 0 &{} 0 &{} 0 &{} 0 \\ 0 &{} -i &{} 0 &{} 0 &{} 0 &{} 0 \\ i &{} 0 &{} 0 &{} 0 &{} 0 &{} 0 \\ 0 &{} 0 &{} 0 &{} 0 &{} 0 &{} 0 \end{array}\right) ,\nonumber \\ \end{aligned}$$
(4.11)

where they satisfy \([T^{ij}, T^{kl}]= i (\delta ^{ik} \, T^{jl}-\delta ^{il} \, T^{jk}- \delta ^{jk} \, T^{il} + \delta ^{jl} \, T^{ik})\).

Sometimes we need to have the (quasi)primary fields in order to express the OPE with known coefficients appearing in the descendant fields. Then we should introduce each primary field at the third, fourth and fifth components of (3.4). The point here is to consider the possible terms with correct spins, SO(4) vector index and SO(4) adjoint index. It turns out that we obtain the following decomposition where the primary field has a tilde

$$\begin{aligned}&\varvec{\Phi }^{(s),\alpha }(Z)\nonumber \\&\quad = \Phi _{0}^{(s),\alpha }(z)+\theta ^{i}\, \Phi _{\frac{1}{2}}^{(s),i,\alpha }(z) \nonumber \\&\qquad +\frac{1}{2} \, \theta ^{4-ij} \,\Phi _{1}^{(s),ij,\alpha }(z) + \theta ^{4-i}\, \Phi _{\frac{3}{2}}^{(s),i,\alpha }(z) + \theta ^{4-0}\, \Phi _{2}^{(s),\alpha }(z) \nonumber \\&\quad \equiv \Phi _{0}^{(s),\alpha }(z)+\theta ^{i}\, \Phi _{\frac{1}{2}}^{(s),i,\alpha }(z)\nonumber \\&\qquad +\frac{1}{2}\theta ^{4-ij}\Bigg [\widetilde{\Phi }_{1}^{(s),ij,\alpha } +\frac{i}{ s}\, ({ T}_R^{ij})^{\alpha \beta }\, \partial \Phi _{0}^{(s),\beta } \Bigg ](z) \nonumber \\&\qquad + \theta ^{4-i}\,\Bigg [ \,\widetilde{\Phi }_{\frac{3}{2}}^{(s),i,\alpha }+ \frac{2\,\alpha }{(2s+1)}\,\partial \Phi _{\frac{1}{2}}^{(s),i,\alpha }\nonumber \\&\qquad -\frac{2i}{(2\, s+1)} \, ({ T}_R^{ij})^{\alpha \beta }\, \partial \Phi _{\frac{1}{2}}^{(s),j,\beta }\, \Bigg ](z) \nonumber \\&\qquad + \theta ^{4-0}\,\Bigg [\, \widetilde{\Phi }_{2}^{(s),\alpha }- p_{1}\,\partial ^{2}\Phi _{0}^{(s),\alpha }-p_{2}\,L\,\Phi _{0}^{(s),\alpha }\nonumber \\&\qquad +\frac{i}{2(2\, s+1)}\,({ T}_L^{ij})^{\alpha \beta } \, \partial \Phi _{1}^{(s),ij,\beta } \nonumber \\&\qquad + ({ T}_L^{ij} { T}_R^{ij})^{\alpha \beta }\,( p_{3}\,\partial ^{2}\Phi _{0}^{(s),\beta } +p_{4}\,L\,\Phi _{0}^{(s),\beta }) \, \Bigg ](z), \end{aligned}$$
(4.12)

where the coefficients in the last component depend on N, k and s (together with (3.2)) and are given by

$$\begin{aligned}&p_{1} = -\frac{2\alpha \,(3+3k+3N+3kN+26s+13ks+13Ns)}{(3+3k+3N+3kN-4s+ks+Ns+6kNs+16s^{2}+8ks^{2}+8Ns^{2})}, \nonumber \\&p_{2} = \frac{12(k-N)s(1+s)}{(3+3k+3N+3kN-4s+ks+Ns+6kNs+16s^{2}+8ks^{2}+8Ns^{2})}, \nonumber \\&p_{3}= \frac{(-15 - 6 k - 6 N + 3 k N + 8 s + 4 k s + 4 N s)}{2 (1 + s) (3 + 3 k + 3 N + 3 k N - 4 s + k s + N s + 6 k N s + 16 s^2 + 8 k s^2 + 8 N s^2)}, \nonumber \\&p_{4}= \frac{3 (2 + k + N)}{(3 + 3 k + 3 N + 3 k N - 4 s + k s + N s + 6 k N s + 16 s^2 + 8 k s^2 + 8 N s^2)}. \end{aligned}$$
(4.13)

Note that the values \(p_1\) and \(p_2\) also appear in the corresponding higher spin-s multiplet in the unitary Wolf space coset model [8]. The expression (4.12) holds for SO(4) nonsinglet higher spin multiplets in the unitary Wolf space coset. Moreover, the following quantities with (3.2) are introduced

$$\begin{aligned} \widetilde{{ T}}^{ij}\equiv & {} \frac{1}{2}\,\varepsilon _{ijkl}\,{ T}^{kl}, \quad { T}^{ij}_{L} \equiv \frac{1}{2} \, { T}^{ij}+{ \alpha }\,{ \widetilde{T}}^{ij}, \nonumber \\ { T}^{ij}_{R}\equiv & {} { \alpha }\,{ T}^{ij}+\frac{1}{2}\,{ \widetilde{T}}^{ij}. \end{aligned}$$
(4.14)

For the SO(4) vector representation \(\alpha \), we will see similar construction in next subsection.

Then we can check that the SO(4) adjoint higher spin multiplet (4.12) satisfies the relation (3.3) with (4.11). Its component relations are given in Appendix A.

The first order pole of (4.8) contains other higher spin currents in various way. For example, the other components of the \(\mathcal{N}=4\) multiplets \(\varvec{\Phi }^{(3),\alpha }(Z)\) can arise.

4.4 The new higher spin currents of spin \(\frac{7}{2}\)

We consider the OPE between the four higher spin-\(\frac{5}{2}\) currents and the higher spin-2 current. It turns out that we obtain

$$\begin{aligned} \Phi _{\frac{1}{2}}^{(2),i}(z)\,\Phi _{0}^{(2)}(w)= & {} \frac{1}{(z-w)^{3}}\, Q^{\left( \frac{3}{2}\right) ,i}_{\frac{1}{2}}(w)\nonumber \\&+\frac{1}{(z-w)^{2}}\,\Bigg [\, \frac{2}{3}\,\partial Q^{\left( \frac{3}{2}\right) ,i}_{\frac{1}{2}} + Q^{\left( \frac{5}{2}\right) ,i}_{\frac{1}{2}} \,\Bigg ](w) \nonumber \\&+ \frac{1}{(z-w)}\,\Bigg [\, \frac{1}{4}\,\partial ^{2} Q^{\left( \frac{3}{2}\right) ,i}_{\frac{1}{2}}\nonumber \\&+ \frac{3}{5}\,\partial Q^{\left( \frac{5}{2}\right) ,i}_{\frac{1}{2}} + Q^{\left( \frac{7}{2}\right) ,i}_{\frac{1}{2}}\, \Bigg ](w)+\cdots .\!\!\!\nonumber \\ \end{aligned}$$
(4.15)

The quasi primary fields appearing in the above poles of (4.15) are given in Appendix B. In particular, the quasi primary field \( Q^{\left( \frac{5}{2}\right) ,i}_{\frac{1}{2}}(w)\) contains \( \Phi ^{(2),i}_{\frac{1}{2}}(w)\) which is the second component of the lowest higher spin-2 multiplet in (2.15). In the first order pole, there exist four new primary fields \( \Phi ^{\left( \frac{7}{2}\right) ,i}_{0}(w)\)Footnote 15 as well as the composite fields containing the SO(4) adjoint \(\Phi _0^{(3),\alpha }(w)\) (and \(\Phi _{\frac{1}{2}}^{(3),i, \alpha }(w)\)), the SO(4) singlet \(\Phi _0^{(2)}(w)\) (\(\Phi _{\frac{1}{2}}^{(2),i}(w)\), \(\Phi _{1}^{(2),ij}(w)\) and \(\Phi _{\frac{3}{2}}^{(2),i}(w)\)) and the other SO(4) singlet \(X_0^{(2)}(w)\) (\(X_{\frac{1}{2}}^{(2),i}(w)\), \(X_{1}^{(2),ij}(w)\) and \(X_{\frac{3}{2}}^{(2),i}(w)\)). In other words, we have

$$\begin{aligned} Q_{\frac{1}{2}}^{\left( \frac{7}{2}\right) ,i}(w) = w_{0,\frac{7}{2}} \, \Phi _0^{\left( \frac{7}{2}\right) ,i}(w) + \cdots . \end{aligned}$$
(4.16)

The abbreviated part is given in Appendix B.

The \(\mathcal{N}=4\) four SO(4) vector higher spin-\(\frac{7}{2}\) multiplets transform under the stress energy tensor as follows [37]. The OPE looks like (3.3) with \(s =\frac{7}{2}\):

$$\begin{aligned} \mathbf{J }(Z_{1}) \, \varvec{\Phi }^{(s),\mu }(Z_{2})= & {} \frac{\theta _{12}^{4}}{z_{12}^{2}} \, 2s\, \varvec{\Phi }^{(s),\mu }(Z_{2})+ \frac{\theta _{12}^{4-i}}{z_{12}} \, D^{i} \varvec{\Phi }^{(s),\mu }(Z_{2})\nonumber \\&+ \frac{\theta _{12}^{4}}{z_{12}} \, 2 \, \partial \varvec{\Phi }^{(s),\mu }(Z_{2}) \nonumber \\&- \frac{i}{2} \,\frac{\theta _{12}^{4-ij}}{z_{12}}({ T}^{ij})^{\mu \nu }\, \varvec{\Phi }^{(s),\nu }(Z_{2}) +\cdots , \nonumber \\ \end{aligned}$$
(4.17)

where the \(T^{ij}\) matrix is the generator of the SO(4) vector representation (4.9)

$$\begin{aligned} { T}^{12}= & {} \left( \begin{array}{cccc} 0 &{} -i &{} 0 &{} 0\\ i &{} 0 &{} 0 &{} 0\\ 0 &{} 0 &{} 0 &{} 0\\ 0 &{} 0 &{} 0 &{} 0 \end{array}\right) , \, { T}^{13}=\left( \begin{array}{cccc} 0 &{} 0 &{} -i &{} 0\\ 0 &{} 0 &{} 0 &{} 0\\ i &{} 0 &{} 0 &{} 0\\ 0 &{} 0 &{} 0 &{} 0 \end{array}\right) , \, \nonumber \\ { T}^{14}= & {} \left( \begin{array}{cccc} 0 &{} 0 &{} 0 &{} -i\\ 0 &{} 0 &{} 0 &{} 0\\ 0 &{} 0 &{} 0 &{} 0\\ i &{} 0 &{} 0 &{} 0 \end{array}\right) , { T}^{23} = \left( \begin{array}{cccc} 0 &{} 0 &{} 0 &{} 0\\ 0 &{} 0 &{} -i &{} 0\\ 0 &{} i &{} 0 &{} 0\\ 0 &{} 0 &{} 0 &{} 0 \end{array}\right) , \, \nonumber \\ { T}^{24}= & {} \left( \begin{array}{cccc} 0 &{} 0 &{} 0 &{} 0\\ 0 &{} 0 &{} 0 &{} -i\\ 0 &{} 0 &{} 0 &{} 0\\ 0 &{} i &{} 0 &{} 0 \end{array}\right) , \, { T}^{34}=\left( \begin{array}{cccc} 0 &{} 0 &{} 0 &{} 0\\ 0 &{} 0 &{} 0 &{} 0\\ 0 &{} 0 &{} 0 &{} -i\\ 0 &{} 0 &{} i &{} 0 \end{array}\right) , \end{aligned}$$
(4.18)

where the commutators \([T^{ij}, T^{kl}]\) satisfy the previous relations described in (4.11). The corresponding component OPEs of (4.17) can be obtained from Appendix A by considering the \(T^{ij}\) matrix as the ones in (4.18). Moreover, we should use the corresponding primary fields (if we need them) according to (4.12) by substituting the \(4\times 4\) matrices in (4.18).

4.5 The higher spin current of spin 4

In the second order pole of (4.1), the quasi primary field \(Q_2^{(4)}(w)\) contains the SO(4) singlet higher spin-4 current \(\Phi _0^{(4)}(w)\). The other part of the quasi primary field \(Q_2^{(4)}(w)\) is given in Appendix B. Note that the composite field \(X_0^{(2)} \, X_0^{(2)}(w)\) (as well as other dependent terms) is absorbed in \(\Phi _0^{(4)}(w)\) such that \(\Phi _0^{(4)}(w)\) should transform as the primary field under the stress energy tensor. If the composite field \(X_0^{(2)} \, X_0^{(2)}(w)\) is not included in the SO(4) singlet higher spin-4 current, then we should calculate the OPE \(X_0^{(2)}(z) \, X_0^{(2)}(w)\) (and their \(\mathcal{N}=4\) version) in order to use its normal ordered product for the Jacobi identity. We observe that all the components of the \(\mathcal{N}=4\) multiplets \(\varvec{\Phi }^{(2)}(Z)\) and \(\mathbf{X}^{(2)}(Z)\) appear in the quasi primary field \(Q_2^{(4)}(w)\).

Moreover, once we combine the above higher spin-4 current with \(\Phi _0^{(2)} \, \Phi _0^{(2)}(w)\) term as well as other terms as in

$$\begin{aligned} \tilde{\Phi }_0^{(4)}(w)= & {} w_{11,4}\,\Phi ^{(2)}_0\Phi ^{(2)}_0(w) +(M^{\alpha })^{ij}(\, c_{1}\,\Phi _{1}^{(3),ij,\alpha } \nonumber \\&+\cdots +\varepsilon ^{ijkl}c_{8}\,\Gamma ^{k}\Gamma ^{l}\Phi _{0}^{(3),\alpha }\,)(w) + \delta _{\mu }^{i}\, c_{9}\,\Phi _{\frac{1}{2}}^{\left( \frac{7}{2}\right) ,i,\mu }(w) \nonumber \\&+ ( c_{10}\,\Phi _{2}^{(2)} +c_{11}\,\tilde{L}\Phi _{0}^{(2)}\nonumber \\&+\cdots +c_{32}\,\varepsilon ^{ijkl}\, T^{ij}T^{kl}\Phi _{0}^{(2)})(w)\nonumber \\&+c_{33}\, X_{2}^{(2)}(w) +c_{34}\,\tilde{L} X_{0}^{(2)}(w) \nonumber \\&+ \cdots + c_{55}\,\varepsilon ^{ijkl}\, T^{ij}T^{kl}X_{0}^{(2)}(w)\nonumber \\&+c_{56}\,\tilde{L}\tilde{L}(w) +c_{57}\,\tilde{L}UU(w) \nonumber \\&+ \cdots +c_{155}\, \varepsilon ^{ijkl}\,\Gamma ^{i} \Gamma ^{j}\partial \Gamma ^{k}\partial \Gamma ^{l}(w), \end{aligned}$$
(4.19)

then the structure constant \(w_{11,4} \equiv C_{(4)(2)}^{(4)}\) does not appear in the remaining OPEs. The coefficients \(c_{10}, \ldots , c_{32}\) in (4.19) also depend on the structure constant \(C_{(4)(2)}^{(2)}\) appeared in the Sect. 4.2. We have checked that this feature arises also in the unitary case [8].

The final first order pole of (4.1) can be obtained and it turns out that there is no new primary field. All the terms after subtracting the descendant fields can be written in terms of the known composite fields (including the higher spin-\(\frac{9}{2}\) currents which will be described in next subsection). We expect that the \(\mathcal{N}=4\) higher spin-5 multiplets (we have not found in this paper) will appear by considering the other OPEs between the \(\mathcal{N}=4\) multiplets we have found in this paper.Footnote 16

4.6 The higher spin currents of spin \(\frac{9}{2}\)

We consider the OPE between the four higher spin-\(\frac{7}{2}\) currents and the higher spin-2 current. It turns out that we obtain

$$\begin{aligned}&\Phi _{\frac{3}{2}}^{(2),i}(z)\,\Phi _{0}^{(2)}(w)\nonumber \\&\quad = \frac{1}{(z-w)^{5}}\, Q^{\left( \frac{1}{2}\right) ,i}_{\frac{3}{2}}(w)\nonumber \\&\qquad +\frac{1}{(z-w)^{4}}\,\Bigg [\, 2\,\partial Q^{\left( \frac{1}{2}\right) ,i}_{\frac{3}{2}}+Q^{\left( \frac{3}{2}\right) ,i}_{\frac{3}{2}} \,\Bigg ](w) \nonumber \\&\qquad +\frac{1}{(z-w)^{3}}\,\Bigg [\, \frac{3}{2}\,\partial ^{2} Q^{\left( \frac{1}{2}\right) ,i}_{\frac{3}{2}} +\partial Q^{\left( \frac{3}{2}\right) ,i}_{\frac{3}{2}} +Q^{\left( \frac{5}{2}\right) ,i}_{\frac{3}{2}} \,\Bigg ](w) \nonumber \\&\qquad + \frac{1}{(z-w)^{2}}\,\Bigg [\, \frac{2}{3}\,\partial ^{3} Q^{\left( \frac{1}{2}\right) ,i}_{\frac{3}{2}} +\frac{1}{2}\,\partial ^{2} Q^{\left( \frac{3}{2}\right) ,i}_{\frac{3}{2}} \nonumber \\&\qquad +\frac{4}{5}\,\partial Q^{\left( \frac{5}{2}\right) ,i}_{\frac{3}{2}} + Q^{\left( \frac{7}{2}\right) ,i}_{\frac{3}{2}} \,\Bigg ](w) \nonumber \\&\qquad + \frac{1}{(z-w)}\,\Bigg [\, \frac{5}{24}\,\partial ^{4} Q^{\left( \frac{1}{2}\right) ,i}_{\frac{3}{2}} +\frac{1}{6}\,\partial ^{3} Q^{\left( \frac{3}{2}\right) ,i}_{\frac{3}{2}} +\frac{1}{3}\,\partial ^{2} Q^{\left( \frac{5}{2}\right) ,i}_{\frac{3}{2}} \nonumber \\&\qquad +\frac{5}{7}\,\partial Q^{\left( \frac{7}{2}\right) ,i}_{\frac{3}{2}} +Q^{\left( \frac{9}{2}\right) ,i}_{\frac{3}{2}} \, \Bigg ](w) \nonumber \\&\qquad + \frac{2\,\alpha }{5}\,\sum \limits _{n=2}^{4} \frac{1}{(z-w)^{n}} \, \Big \{ \partial \Phi _{\frac{1}{2}}^{(2),i} \,\Phi _{0}^{(2)} \Big \}_{n}(w) +\cdots . \end{aligned}$$
(4.20)

Again, the OPE (4.20) consists of two parts. The first four lines comes from the corresponding four primary higher spin-\(\frac{7}{2}\) currents \(\tilde{\Phi }_{\frac{3}{2}}^{(2),i}(z)\) and the last line comes from the additional term described in (4.12).Footnote 17

By using the relation (by using the notation of [39])

$$\begin{aligned} \Big \{ \partial \Phi _{\frac{1}{2}}^{(2),i} \,\Phi _{0}^{(2)} \Big \}_{n+1}(w) = -n\, \Big \{ \Phi _{\frac{1}{2}}^{(2),i} \,\Phi _{0}^{(2)} \Big \}_{n}(w), \quad n=1,2,3, \end{aligned}$$

where the previous relation (4.15) can be used, we can present the above OPE as follows:

$$\begin{aligned}&\Phi _{\frac{3}{2}}^{(2),i}(z)\,\Phi _{0}^{(2)}(w)\nonumber \\&\quad = \frac{1}{(z-w)^{5}}\, Q^{\left( \frac{1}{2}\right) ,i}_{\frac{3}{2}}(w)\nonumber \\&\qquad +\frac{1}{(z-w)^{4}}\,\Bigg [\, 2\,\partial Q^{\left( \frac{1}{2}\right) ,i}_{\frac{3}{2}}+Q^{\left( \frac{3}{2}\right) ,i}_{\frac{3}{2}} + R^{\left( \frac{3}{2}\right) ,i}_{\frac{3}{2}} \,\Bigg ](w) \nonumber \\&\qquad +\frac{1}{(z-w)^{3}}\,\Bigg [\, \frac{3}{2}\,\partial ^{2} Q^{\left( \frac{1}{2}\right) ,i}_{\frac{3}{2}} +\partial Q^{\left( \frac{3}{2}\right) ,i}_{\frac{3}{2}} +Q^{\left( \frac{5}{2}\right) ,i}_{\frac{3}{2}} +R^{\left( \frac{5}{2}\right) ,i}_{\frac{3}{2}} \,\Bigg ](w) \nonumber \\&\qquad + \frac{1}{(z-w)^{2}}\,\Bigg [\, \frac{2}{3}\,\partial ^{3} Q^{\left( \frac{1}{2}\right) ,i}_{\frac{3}{2}} +\frac{1}{2}\,\partial ^{2} Q^{\left( \frac{3}{2}\right) ,i}_{\frac{3}{2}}\nonumber \\&\qquad +\frac{4}{5}\,\partial Q^{\left( \frac{5}{2}\right) ,i}_{\frac{3}{2}} + Q^{\left( \frac{7}{2}\right) ,i}_{\frac{3}{2}} +R^{\left( \frac{7}{2}\right) ,i}_{\frac{3}{2}} \,\Bigg ](w) \nonumber \\&\qquad + \frac{1}{(z-w)}\,\Bigg [\, \frac{5}{24}\,\partial ^{4} Q^{\left( \frac{1}{2}\right) ,i}_{\frac{3}{2}} +\frac{1}{6}\,\partial ^{3} Q^{\left( \frac{3}{2}\right) ,i}_{\frac{3}{2}} +\frac{1}{3}\,\partial ^{2} Q^{\left( \frac{5}{2}\right) ,i}_{\frac{3}{2}}\nonumber \\&\qquad +\frac{5}{7}\,\partial Q^{\left( \frac{7}{2}\right) ,i}_{\frac{3}{2}} +Q^{\left( \frac{9}{2}\right) ,i}_{\frac{3}{2}} \, \Bigg ](w) \nonumber \\&\qquad + \cdots , \end{aligned}$$
(4.21)

where we introduce

$$\begin{aligned} R_{\frac{3}{2}}^{\left( \frac{7}{2}\right) ,i}(w)\equiv & {} -\frac{2 \, \alpha }{5}\, \left( \frac{1}{4} \, \partial ^2 \, Q_{\frac{1}{2}}^{\left( \frac{3}{2}\right) ,i} +\frac{3}{5} \, \partial \, Q_{\frac{1}{2}}^{\left( \frac{5}{2}\right) ,i} + Q_{\frac{1}{2}}^{\left( \frac{7}{2}\right) ,i}\right) (w), \nonumber \\ R_{\frac{3}{2}}^{\left( \frac{5}{2}\right) ,i}(w)\equiv & {} -\frac{4 \, \alpha }{5} \, \left( \frac{2}{3} \, \partial \, Q_{\frac{1}{2}}^{\left( \frac{3}{2}\right) ,i} + Q_{\frac{1}{2}}^{\left( \frac{5}{2}\right) ,i}\right) (w), \nonumber \\ R_{\frac{3}{2}}^{\left( \frac{3}{2}\right) ,i}(w)\equiv & {} -\frac{6 \, \alpha }{5} \, Q_{\frac{1}{2}}^{\left( \frac{3}{2}\right) ,i}(w), \end{aligned}$$
(4.22)

from the last line of (4.20). We can see that the fifth and fourth order terms of (4.21) consist of the composite fields from the 16 currents as in Appendix B. The third and second order poles of (4.21) have the higher spin currents (found before) dependent terms. The first order pole has the following form

$$\begin{aligned} Q_{\frac{3}{2}}^{\left( \frac{9}{2}\right) ,i}(w)= & {} w_{0,\frac{9}{2}} \, \Phi _0^{\left( \frac{9}{2}\right) ,i}(w) + \cdots , \end{aligned}$$
(4.23)

where the other remaining terms are given in Appendix B. Each component of the higher spin-\(\frac{9}{2}\) current can be obtained by following the procedure described in the footnote 14. We see that all kinds of higher spin currents appear in this quasi primary field. We can use (4.17) together with the above \(4\times 4\) matrices in (4.18) in order to obtain the OPEs between this higher spin multiplet and the \(\mathcal{N}=4\) stress energy tensor. Its component results are given in Appendix A.

4.7 The fundamental 16 OPEs

Therefore, the fundamental 16 OPEs (five kinds of OPEs) are given by (2.13), (4.7), (4.8), (4.15) and (4.21). The structure constants are written in terms of N, k and \(C_{(4)(2)}^{(2)}\) which will be given in the \(\mathtt{ancillary.nb}\). These will determine the remaining 120 OPEs by using the \(\mathcal{N}=4\) supersymmetry soon.

We summarize the higher spin currents appearing in various quasi primary fields (4.2) of these OPEs in Table 1. Some of the higher spin currents are not present in this Table 1 and they will arise in the remaining 120 OPEs. For simplicity, we do not include the dependence of 16 currents of the large \(\mathcal{N}=4\) linear superconformal algebra.

Table 1 The structure constant \(w_{1,2} \equiv C_{(4)(2)}^{(2)}\) appears in front of the above quasi primaries appearing in the OPEs except \(Q_{\frac{1}{2}}^{\left( \frac{5}{2}\right) ,i}(z)\) and \(Q_2^{(3)}(z)\), which do not have the components of \(\mathcal{N}=4\) multiplet \(\mathbf{X}^{(2)}(Z)\). The quasi primaries, which are not in this list, have the composite fields consisting of 16 currents. Although the \(X_0^{(2)}(z)\) dependence in \(Q_2^{(3)}(z)\) does not appear, its dependence appears in \(\partial \, Q_2^{(2)}(z)\) in (4.7). Similarly, the \(\Phi _0^{(4)}(z)\) dependence in \(Q_2^{(5)}(z)\) does not appear and its dependence appears via \(\partial \, Q_2^{(4)}(z)\) in (4.7)

4.8 One single \(\mathcal{N}=4\) super OPE

From the fundamental 16 OPEs (2.13), (4.15), (4.8), (4.21) and (4.7) (that is, five different kinds of OPEs), we can generalize them in \(\mathcal{N}=4\) superspace by taking [8]

$$\begin{aligned} U(w)\rightarrow & {} \partial \, \mathbf{J}(Z_2), \nonumber \\ \Gamma ^i(w)\rightarrow & {} -i \, D^i \mathbf{J}(Z_2) \equiv -i \, \mathbf{J}^i(Z_2), \nonumber \\ T^{ij}(w)\rightarrow & {} -\frac{i}{2!} \, \varepsilon ^{ijkl} \, D^k D^l \mathbf{J}(Z_2) \equiv -\frac{i}{2!} \, \varepsilon ^{ijkl} \, \mathbf{J}^{kl}(Z_2), \nonumber \\ \tilde{G}^i(w)\rightarrow & {} \frac{1}{3!} \, \varepsilon ^{ijkl} \, D^j D^k D^l \mathbf{J}(Z_2) \equiv \frac{1}{3!} \,\varepsilon ^{ijkl} \, \mathbf{J}^{jkl}(Z_2), \nonumber \\ \tilde{L}(w)\rightarrow & {} \frac{1}{2 \cdot 4!} \, \varepsilon ^{ijkl} \, D^{i} D^j D^k D^l \mathbf{J}(Z_2) \nonumber \\&\equiv \frac{1}{2 \cdot 4!} \,\varepsilon ^{ijkl} \, \mathbf{J}^{ijkl}(Z_2), \nonumber \\ \Phi _{0}^{(s),\alpha }(w)\rightarrow & {} \varvec{\Phi }^{(s),\alpha }(Z_2), \nonumber \\ \Phi _{\frac{1}{2}}^{(s),i,\alpha }(w)\rightarrow & {} D^i \varvec{\Phi }^{(s),\alpha }(Z_2), \nonumber \\ \Phi _{1}^{(s),ij,\alpha }(w)\rightarrow & {} -\frac{1}{2!} \, \varepsilon ^{ijkl} \, D^k D^l \varvec{\Phi }^{(s),\alpha }(Z_2), \nonumber \\ \Phi _{\frac{3}{2}}^{(s),i,\alpha }(w)\rightarrow & {} -\frac{1}{3!} \, \varepsilon ^{ijkl} \, D^j D^k D^l \varvec{\Phi }^{(s),\alpha }(Z_2), \nonumber \\ \Phi _{2}^{(s),\alpha }(w)\rightarrow & {} \frac{1}{ 4!} \, \varepsilon ^{ijkl} \, D^{i} D^j D^k D^l \varvec{\Phi }^{(s),\alpha }(Z_2), \qquad \nonumber \\&\alpha = \hbox {singlet, }\,\,\hbox { adjoint, }\,\,\hbox { vector}, \end{aligned}$$
(4.24)

where \(\tilde{G}^i(w)\) and \(\tilde{L}(w)\) are given in (3.1) and putting the relevant fermionic coordinates. For the singlet, adjoint and vector representations, we substitute the corresponding indices into the \(\alpha \). In doing this, there are additional terms arising from the summation over the same indices. We present the quasi primary fields in \(\mathcal{N}=4\) superspace in Appendix C. The total number of terms in Appendix B and Appendix C is little different from each other.

Then the single \(\mathcal{N}=4\) super OPE between the SO(4) singlet higher spin-2 multiplet can be summarized by (after rearranging (D.1))

$$\begin{aligned}&\varvec{\Phi }^{(2)}(Z_{1})\,\varvec{\Phi }^{(2)}(Z_{2}) = \frac{1}{z_{12}^{4}}\, c^{0,4}_{0} +\frac{\theta _{12}^{4-0}}{z_{12}^{6}}\, 8 \, \alpha \, c_{0}^{0,4} +\frac{\theta _{12}^{4-i}}{z_{12}^{5}}\, \mathbf{Q}^{\left( \frac{1}{2}\right) ,i}_{\frac{3}{2}}(Z_{2}) \nonumber \\&\quad +\frac{\theta _{12}^{4-0}}{z_{12}^{5}}\, \mathbf{Q}^{(1)}_{2}(Z_{2}) + \frac{\theta _{12}^{4-ij}}{z_{12}^{4}}\, \mathbf{Q}^{(1),ij}_{1}(Z_{2})\nonumber \\&\quad +\frac{\theta _{12}^{4-i}}{z_{12}^{4}}\,\Bigg [\, 2\,\partial \mathbf{Q}^{\left( \frac{1}{2}\right) ,i}_{\frac{3}{2}}+\mathbf{Q}^{\left( \frac{3}{2}\right) ,i}_{\frac{3}{2}} +\mathbf{R}^{\left( \frac{3}{2}\right) ,i}_{\frac{3}{2}} \,\Bigg ](Z_{2}) \nonumber \\&\quad +\frac{\theta _{12}^{4-0}}{z_{12}^{4}}\,\Bigg [\, \frac{3}{2}\,\partial \mathbf{Q}^{(1)}_{2} +\mathbf{Q}^{(2)}_{2} +\mathbf{R}^{(2)}_{2}\,\Bigg ](Z_{2})\nonumber \\&\quad + \frac{\theta _{12}^{i}}{z_{12}^{3}}\, \mathbf{Q}^{\left( \frac{3}{2}\right) ,i}_{\frac{1}{2}}(Z_{2}) + \frac{\theta _{12}^{4-ij}}{z_{12}^{3}}\,\Bigg [\, \partial \mathbf{Q}^{(1),ij}_{1} +\mathbf{Q}^{(2),ij}_{1} \,\Bigg ](Z_{2}) \nonumber \\&\quad +\frac{\theta _{12}^{4-i}}{z_{12}^{3}}\,\Bigg [\, \frac{3}{2}\,\partial ^{2} \mathbf{Q}^{\left( \frac{1}{2}\right) ,i}_{\frac{3}{2}} +\partial \mathbf{Q}^{\left( \frac{3}{2}\right) ,i}_{\frac{3}{2}} +\mathbf{Q}^{\left( \frac{5}{2}\right) ,i}_{\frac{3}{2}} +\mathbf{R}^{\left( \frac{5}{2}\right) ,i}_{\frac{3}{2}} \,\Bigg ](Z_{2}) \nonumber \\&\quad +\frac{\theta _{12}^{4-0}}{z_{12}^{3}}\,\Bigg [\, \partial ^{2} \mathbf{Q}^{(1)}_{2} +\partial \mathbf{Q}^{(2)}_{2} +\mathbf{Q}^{(3)}_{2} +\mathbf{R}^{(3)}_{2}\,\Bigg ](Z_{2}) +\frac{1}{z_{12}^{2}}\, \mathbf{Q}^{(2)}_{0}(Z_{2}) \nonumber \\&\quad + \frac{\theta _{12}^{i}}{z_{12}^{2}}\,\Bigg [\, \frac{2}{3}\,\partial \mathbf{Q}^{\left( \frac{3}{2}\right) ,i}_{\frac{1}{2}} + \mathbf{Q}^{\left( \frac{5}{2}\right) ,i}_{\frac{1}{2}} \,\Bigg ](Z_{2})\nonumber \\&\quad + \frac{\theta _{12}^{4-ij}}{z_{12}^{2}}\,\Bigg [\, \frac{1}{2}\,\partial ^{2} \mathbf{Q}^{(1),ij}_{1} +\frac{3}{4}\,\partial \mathbf{Q}^{(2),ij}_{1} +\mathbf{Q}^{(3),ij}_{1} \,\Bigg ](Z_{2}) \nonumber \\&\quad + \frac{\theta _{12}^{4-i}}{z_{12}^{2}}\,\Bigg [\, \frac{2}{3}\,\partial ^{3} \mathbf{Q}^{\left( \frac{1}{2}\right) ,i}_{\frac{3}{2}} +\frac{1}{2}\,\partial ^{2} \mathbf{Q}^{\left( \frac{3}{2}\right) ,i}_{\frac{3}{2}}\nonumber \\&\quad +\frac{4}{5}\,\partial \mathbf{Q}^{\left( \frac{5}{2}\right) ,i}_{\frac{3}{2}} + \mathbf{Q}^{\left( \frac{7}{2}\right) ,i}_{\frac{3}{2}} +\mathbf{R}^{\left( \frac{7}{2}\right) ,i}_{\frac{3}{2}} \,\Bigg ](Z_{2}) \nonumber \\&\quad +\frac{\theta _{12}^{4-0}}{z_{12}^{2}}\,\Bigg [\, \frac{5}{12}\,\partial ^{3} \mathbf{Q}^{(1)}_{2} +\frac{1}{2}\, \partial \mathbf{Q}^{(2)}_{2} +\frac{5}{6}\,\mathbf{Q}^{(3)}_{2} +\mathbf{Q}^{(4)}_{2} +\mathbf{R}^{(4)}_{2} \,\Bigg ](Z_{2}) \nonumber \\&\quad +\frac{1}{z_{12}}\, \frac{1}{2}\partial \mathbf{Q}^{(2)}_{0}(Z_{2}) + \frac{\theta _{12}^{i}}{z_{12}}\,\Bigg [\, \frac{1}{4}\,\partial ^{2} \mathbf{Q}^{\left( \frac{3}{2}\right) ,i}_{\frac{1}{2}} + \frac{3}{5}\,\partial \mathbf{Q}^{\left( \frac{5}{2}\right) ,i}_{\frac{1}{2}} + \mathbf{Q}^{\left( \frac{7}{2}\right) ,i}_{\frac{1}{2}} \,\Bigg ](Z_{2}) \nonumber \\&\quad + \frac{\theta _{12}^{4-ij}}{z_{12}}\,\Bigg [\, \frac{1}{6}\,\partial ^{3} \mathbf{Q}^{(1),ij}_{1} +\frac{3}{10}\,\partial ^{2} \mathbf{Q}^{(2),ij}_{1} +\frac{2}{3}\,\partial \mathbf{Q}^{(3),ij}_{1} +\mathbf{Q}^{(4),ij}_{1} \,\Bigg ](Z_{2}) \nonumber \\&\quad +\frac{\theta _{12}^{4-i}}{z_{12}}\,\Bigg [\, \frac{5}{24}\,\partial ^{4} \mathbf{Q}^{\left( \frac{1}{2}\right) ,i}_{\frac{3}{2}} +\frac{1}{6}\,\partial ^{3} \mathbf{Q}^{\left( \frac{3}{2}\right) ,i}_{\frac{3}{2}} +\frac{1}{3}\,\partial ^{2} \mathbf{Q}^{\left( \frac{5}{2}\right) ,i}_{\frac{3}{2}}\nonumber \\&\quad +\frac{5}{7}\,\partial \mathbf{Q}^{\left( \frac{7}{2}\right) ,i}_{\frac{3}{2}} + \mathbf{Q}^{\left( \frac{9}{2}\right) ,i}_{\frac{3}{2}} \,\Bigg ](Z_{2}) \nonumber \\&\quad +\frac{\theta _{12}^{4-0}}{z_{12}}\,\Bigg [\, \frac{1}{8}\,\partial ^{4} \mathbf{Q}^{(1)}_{2} +\frac{1}{6}\,\partial ^{3} \mathbf{Q}^{(2)}_{2} +\frac{5}{14}\,\partial ^{2} \mathbf{Q}^{(3)}_{2} \nonumber \\&\quad +\frac{3}{4}\,\partial \mathbf{Q}^{(4)}_{2} + \mathbf{Q}^{(5)}_{2}\, \Bigg ](Z_{2}) +\cdots . \end{aligned}$$
(4.25)

The central term is presented in (2.14). The maximum number of super spin is given by 5 and the corresponding (composite) higher spin currents appear in the last line of (4.25). Due to the space of the paper, we cannot write down all the operators in the right hand side of (4.25). The partial expressions of quasi primary super fields corresponding to the component fields in (4.2) are given in Appendix C (together with \(\mathtt{ancillary.nb}\)).

Here we introduce the following quantities (\(\mathcal{N}=4\) expressions of (4.22), (4.6) and the one in the footnote 12)Footnote 18 with (3.2)

$$\begin{aligned} \mathbf{R}_{\frac{3}{2}}^{\left( \frac{7}{2}\right) ,i}(Z_2)\equiv & {} -\frac{2 \, \alpha }{5}\, (\frac{1}{4} \, \partial ^2 \, \mathbf{Q}_{\frac{1}{2}}^{\left( \frac{3}{2}\right) ,i} +\frac{3}{5} \, \partial \, \mathbf{Q}_{\frac{1}{2}}^{\left( \frac{5}{2}\right) ,i} + \mathbf{Q}_{\frac{1}{2}}^{\left( \frac{7}{2}\right) ,i})(Z_2),\\ \mathbf{R}_{\frac{3}{2}}^{\left( \frac{5}{2}\right) ,i}(Z_2)\equiv & {} -\frac{4 \, \alpha }{5} \, (\frac{2}{3} \, \partial \, \mathbf{Q}_{\frac{1}{2}}^{\left( \frac{3}{2}\right) ,i} + \mathbf{Q}_{\frac{1}{2}}^{\left( \frac{5}{2}\right) ,i})(Z_2),\\ \mathbf{R}_{\frac{3}{2}}^{\left( \frac{3}{2}\right) ,i}(Z_2)\equiv & {} -\frac{6 \, \alpha }{5} \, \mathbf{Q}_{\frac{1}{2}}^{\left( \frac{3}{2}\right) ,i}(Z_2),\\ \mathbf{R}^{(4-n)}_{2}(Z_2)\equiv & {} - p_2 \, \mathbf{E}^{(4-n)}_{2}(Z_2)\\&-p_{1}\,n(n+1)\,\mathbf{Q}^{(4-n)}_{0}(Z_2),\qquad n =0,1,2. \end{aligned}$$

Note that \(\mathbf{R}_{\frac{3}{2}}^{\left( \frac{7}{2}\right) ,i}(Z_2)\), \(\mathbf{R}_{\frac{3}{2}}^{\left( \frac{5}{2}\right) ,i}(Z_2)\), \(\mathbf{R}_2^{(3)}(Z_2)\) and \(\mathbf{R}_2^{(4)}(Z_2)\) are not quasi primary. From (4.3), (4.10), (4.16), (4.19) and (4.23), the new primary higher spin \({ \mathcal N}=4\) multiplets in (4.25) arise in the following quasi primary ones

$$\begin{aligned} \mathbf{Q}^{(2)}_{2}(Z_2)= & {} c_{1}^{2,4}\,\varvec{\Phi }^{(2)}(Z_2) +c_{2}^{2,4}\,\mathbf{X}^{(2)}(Z_2) + \cdots , \nonumber \\ \mathbf{Q}^{(3),ij}_{1}(Z_2)= & {} c_{0}^{1,2}\,(M^{\alpha })^{ij}\varvec{\Phi }^{(3),\alpha }(Z_2) + \cdots , \nonumber \\ \mathbf{Q}^{\left( \frac{7}{2}\right) ,i}_{\frac{1}{2}}(Z_2)= & {} c_{0}^{\frac{1}{2},1}\,\varvec{\Phi }^{\left( \frac{7}{2}\right) ,i}(Z_2) +\cdots ,\nonumber \\ \mathbf{Q}^{(4)}_{2}(Z_2)= & {} c_{0}^{2,2}\,\varvec{\Phi }^{(4)}(Z_2) + \cdots , \nonumber \\ \mathbf{Q}^{\left( \frac{9}{2}\right) ,i}_{\frac{3}{2}}(Z_2)= & {} c_{0}^{\frac{3}{2},1}\,\varvec{\Phi }^{\left( \frac{9}{2}\right) ,i}(Z_2) +\cdots , \end{aligned}$$
(4.26)

where some of the abbreviated parts are given in Appendix C. Its component relations can be found in Appendix B. We can easily figure out the SO(4) indices and it is rather nontrivial to observe the second quasi primary field of (4.26) where SO(4) adjoint index is contracted with the one in the matrix M (4.9).

Schematically, we can present the above OPE as follows:

$$\begin{aligned} {[}\varvec{\Phi }^{(2)} \cdot \varvec{\Phi }^{(2)}]= & {} [\mathbf{I}] + \theta ^{4-0}\, ([\varvec{\Phi }^{(2)}] \nonumber \\&+ [\mathbf{X}^{(2)}]) + \theta ^{4-ij}\, (M^{\alpha })^{ij} \, [\varvec{\Phi }^{(3),\alpha }]\nonumber \\&+ \theta ^{i}\, [\varvec{\Phi }^{\left( \frac{7}{2}\right) ,i}]+ \theta ^{4-0}\, [{\tilde{\varvec{\Phi }}}^{(4)}] \nonumber \\&+ \theta ^{4-i}\, [\varvec{\Phi }^{\left( \frac{9}{2}\right) ,i}], \end{aligned}$$
(4.27)

where \([\mathbf{I}]\) stands for the large \(\mathcal{N}=4\) linear superconformal family of identity operator and the various composite fields consisting of \(\mathcal{N}=4\) multiplet \(\mathbf{J}\) up to the super spin 5 can appear. We also insert the SO(4) vector indices i, j and the SO(4) matrix \(M^{\alpha }\) in (4.9). Note that in the component approach described in previous subsections, the SO(4) indices are present in the two higher spin currents of the left hand side of the given OPE. In the \(\mathcal{N}=4\) superspace description (4.27), all the SO(4) indices are contracted with the ones in the fermionic coordinates. In order to obtain the component results, we can act various super derivatives both sides of the \(\mathcal{N}=4\) OPE to restore the SO(4) indices. Here the \(\mathcal{N}=4\) multiplet \({\tilde{\varvec{\Phi }}}^{(4)}(Z)\) has its lowest component given in (4.19) and contains the quadratic \(\mathcal{N}=4\) multiplet \(\varvec{\Phi }^{(2)}(Z)\). Compared with the \(\mathcal{N}=3\) example [14], there are more higher spin multiplets contracted with fermionic coordinates.

4.9 The 136 OPEs between the 16 lowest higher spin currents for generic N

We can calculate the remaining \(136-16=120\) OPEs from (4.25) by taking the super derivatives \(D^i_1\) or \(D^j_2\) both sides of (4.25) and putting \(\theta _1^k=0=\theta ^l_2\). By introducing the various coefficients in front of composite fields appearing in 136 OPEs (the number of coefficients is 2000 or so and coefficients are denoted by \(w_{1,s}, \ldots , w_{2043,s}\) in Appendix B) and using the Jacobi identities, we obtain the 136 OPEs between the 16 lowest higher spin currents (2.15) for generic N. Note that this number of coefficients is huge compared to the unitary case in [8]. Furthermore, after substituting these coefficients into (4.25) back, we obtain the final single OPE with fixed coefficients which depend on N and k (as well as \(C_{(4)(2)}^{(2)}\)).Footnote 19 Because there are new 16 \(\mathcal{N}=4\) multiplets (\(\mathbf{X}^{(2)}(Z)\), \(\varvec{\Phi }^{(3),\alpha }(Z)\), \(\varvec{\Phi }^{\left( \frac{7}{2}\right) ,\mu }(Z)\), \(\varvec{\Phi }^{(4)}(Z)\) and \(\varvec{\Phi }^{\left( \frac{9}{2}\right) ,\mu }(Z)\)) in the right hand side of the 136 OPEs, we have \(16 \times 16=256\) higher spin currents (in components) totally. Under the 16 currents of the large \(\mathcal{N}=4\) superconformal algebra, they transform nontrivially as in Appendix A. It is an open problem to determine the OPEs between these 256 higher spin currents (or the OPEs between the lowest 16 higher spin currents and those 256 higher spin currents) systematically.

5 Conclusions and outlook

We have described one single \(\mathcal{N}=4\) super OPE (4.25) between the lowest higher spin-2 multiplet in the \(\mathcal{N}=4\) superspace. As in the abstract, there exist several \(\mathcal{N}=4\) higher spin multiplets in the right hand side of this OPE.

There are open problems we can consider in the future as follows:

\(\bullet \) Higher spin algebra in the bulk theory

In [42], the free field construction at \(\lambda =0\) by using the bosons and fermions is presented. Maybe at this particular \(\lambda =0\) case, the full higher spin algebra can be described. In other words, the commutators and anticommutators for the higher spin currents (including the 16 currents) can be determined with complete structure constants. The final goal is to obtain the higher spin algebra at finite \(\lambda \) which will provide the corresponding algebra in the dual conformal field theory at the classical level. Contrary to the unitary case, the orthogonal case needs to obtain the appropriate truncation on the matrix elements observed in [42].

\(\bullet \) Three-point functions

One way to check the dual relation between the orthogonal Wolf space coset model and the higher spin theory on \(AdS_3\) space is to compute the three-point functions of the two scalars and the higher spin currents. According to the results of this paper, there are many higher spin currents from the single OPE (4.25). It is an open problem to obtain the remaining 15 higher spin currents in terms of the orthogonal Wolf space coset fields explicitly and to calculate the eigenvalues of the zero modes, by following the procedure studied in [43,44,45,46,47]. Although this will be rather involved, once we obtain them, then it is straightforward to compute the three-point functions at finite N and k.

\(\bullet \) \(\mathcal{N}=2\) superspace description

In principle, we can rewrite the above \(\mathcal{N}=4\) superspace OPE in terms of various 10 OPEs in \(\mathcal{N}=2\) superspace. One merit for this description is that contrary to the \(\mathcal{N}=4\) superspace OPE we have described so far, the \(\mathcal{N}=2\) superspace description enables us to write down in terms of quasi (super) primary fields completely. In doing this, it is rather nontrivial to obtain the correct component fields for the SO(4) nonsinglet \(\mathcal{N}=4\) multiplets. The relevant work in this direction appeared in [14].

\(\bullet \) The large k limit

We can examine the behavior of large k limit (for example, see the work of [48]) from the results we have obtained in this paper. We take the large k limit in the structure constants appearing in the right hand sides of the OPEs we have found. We can read off the leading behavior of k of the right hand sides. Even the \(N=5\) results are enough to analyze this large k limit. We expect to observe the realization of vanishing of ’t Hooft-like coupling constant \(\lambda =\frac{(N+1)}{(N+k+2)} \rightarrow 0\) for fixed N [2]. It is also interesting to observe whether there is an extension of the small \(\mathcal{N}=4\) superconformal algebra along the line of [2] or not.

\(\bullet \) The nonlinear version

It is an open problem to obtain the above single OPE (4.25) in the context of nonlinear version which is an extension of the large \(\mathcal{N}=4\) nonlinear algebra. Due to the fact that we do not know its \(\mathcal{N}=4\) superspace version, we need to present the whole 136 OPEs. Although the lowest 16 higher spin currents are primary under the corresponding stress energy tensor, they do transform nontrivially with respect to other 10 currents. In principle, because we do have the complete OPE in the linear version, it is straightforward to obtain them in the nonlinear version although the careful analysis should be done. In this paper, we applied some computations in this nonlinear basis, although we did not present them explicitly (some OPEs are rather complicated).