Seed Conformal Blocks in 4D CFT

We compute in closed analytical form the minimal set of"seed"conformal blocks associated to the exchange of generic mixed symmetry spinor/tensor operators in an arbitrary representation (l,\bar l) of the Lorentz group in four dimensional conformal field theories. These blocks arise from 4-point functions involving two scalars, one (0,|l-\bar l|) and one (|l-\bar l|,0) spinors or tensors. We directly solve the set of Casimir equations, that can elegantly be written in a compact form for any (l,\bar l), by using an educated ansatz and reducing the problem to an algebraic linear system. Various details on the form of the ansatz have been deduced by using the so called shadow formalism. The complexity of the conformal blocks depends on the value of p=|l-\bar l | and grows with p, in analogy to what happens to scalar conformal blocks in d even space-time dimensions as d increases. These results open the way to bootstrap 4-point functions involving arbitrary spinor/tensor operators in four dimensional conformal field theories.


Introduction
There has recently been a revival of interest in the old idea of the bootstrap program [1,2] after ref. [3] observed that its applicability extends to Conformal Field Theories (CFTs) in more than two space-time dimensions (2D). Since ref. [3], several constraints have been imposed on CFTs data, namely spectrum of operators and Operator Product Expansion (OPE) coefficients, in CFTs in different dimensions, up to 6D (see e.g. ref. [4]). Imposing additional (mild and reasonable) assumptions, one can also compute CFT data of given CFTs, the 3D Ising model being probably the most striking example [5][6][7].
The bootstrap approach is a systematic way of imposing crossing symmetry in correlation functions. Analyzing 4-point functions is enough to get constraints on the CFT data. In order to be able to implement the bootstrap program, it is essential to be able to decompose the 4-point functions in terms of the individual contributions of the infinite number of primary operators (and all their descendants) that can be exchanged in any given kinematical channel, namely that appear in the OPE of the four external operators. For each primary operator, in particular, one has to resum the infinite series of associated descendant operators in what is called a Conformal Partial Wave (CPW). CPWs can be decomposed in terms of scalar functions known as Conformal Blocks (CBs). Before the advent of ref. [3], the only known CBs were those associated to symmetric traceless tensors exchanged in scalar 4-point functions in even number of dimensions [8,9], denoted for short scalar symmetric CBs in the following.
Not surprisingly, after ref. [3] significant progress has been made in computing CBs. Various techniques have been introduced to determine in an expanded form the scalar symmetric CBs in 3D [5,[10][11][12], where a general closed analytic expression has not been found so far. In particular, using the techniques of ref. [12] and the further developments in ref. [13], CBs associated to the exchange of fermion operators in 3D have recently been computed [14]. In ref. [15] it has been shown how to relate, in any number of dimensions, symmetric CBs in correlators of external traceless symmetric operators to the known scalar symmetric blocks. In ref. [16] the so called shadow formalism method [17][18][19][20], already used in ref. [8], has been further developed to compute any CB in any number of dimensions. Although very powerful, the shadow formalism leads to quite involved and not so enlightening expressions. Applications of this method for some specific correlators appeared in ref. [21]. Some other limits of the known CBs have been discussed in refs. [22,23], as well as their interpretation in terms of Witten diagrams in Anti de Sitter (AdS) space [24]. Despite significant progress, not much has been done in the analysis of CBs associated to mixed symmetry tensor (or fermion) operators, denoted simply mixed tensor CBs in the following. Such CBs are crucial to extend the bootstrap program to tensor correlators in CFTs in d > 3 space-time dimensions, where such operators can appear in the OPE between two external fields. 1 Mixed tensor CBs in 4D CFTs have recently been analyzed in ref. [25]. In particular it has been shown there how to relate, by means of differential operators, mixed tensor CBs appearing in an arbitrary spinor/tensor 4-point correlator (not necessarily traceless symmetric operators) to a basis of minimal mixed tensor CBs. These "seed" blocks arise from 4-point functions involving two scalars and two tensor fields in the (0, p) and (p, 0) representations of the Lorentz group, with p an arbitrary integer ((1, 0) is a fermion). Such 4-point functions are the simplest ones (i.e. with the least number of tensor structures) where ( + p, ) or ( , + p) mixed symmetry (bosonic or fermionic) tensors can be exchanged in some OPE limit, for any .
The aim of this paper is to compute the "seed" CBs identified in ref. [25]. We will be able to find in closed analytical form the set of seed CBs associated to the exchange of operators in the ( ,¯ ) representations of the Lorentz group. They are labelled by the positive integer p = | −¯ | and are thus infinite. We consider at the same time CBs associated to both bosonic (even p) and fermionic (odd p) mixed symmetry tensor operators. For each given p, one has to determine p + 1 CBs G (p) e , e = 0, 1, . . . , p, one for each tensor structure appearing in the corresponding CPW. Using the 6D embedding formalism in twistor space in index-free notation [16,26], we will be able to write in a compact form the system of Casimir equations satisfied by the p + 1 CBs. Solving the Casimir system is a hard task, that also requires the knowledge of some boundary conditions, like the asymptotic behaviour of G (p) e . We first attack the problem using the shadow formalism. With the use of some tricks, we find integral expressions of G (p) e for any p and , and explicit expressions for p = 1, 2 (and any ). The shadow formalism also allows us to get the asymptotic behaviour of G (p) e in the OPE limit u → 0, v → 1 for p = 1, 2 and any , and for any p and = 0, together with some other information on the structure of the blocks. Thanks to the knowledge acquired in this way, we will be able to go back to the Casimir system and solve it for any p and , using generalizations of the methods introduced in ref. [9] (and further refined in ref. [27]) to compute 6D symmetric CBs for scalar correlators. Like scalar blocks in higher even dimensions, the mixed tensor CBs are found using an ansatz given by a sum of hypergeometric functions with unknown coefficients c e m,n . In this way a system of p + 1 linear coupled differential equations of second order in two variables is reduced to an algebraic linear system for c e m,n . The set of non-trivial coefficients c e m,n , determined by solving the linear system, admits a useful geometric interpretation. They span a two-dimensional lattice in the (m, n) plane. For each CB labelled by e, the shape of the lattice is an octagon, with p and e dependent edges. For large p, the total number of coefficients c e m,n grows like p 3 and their explicit form becomes more and more complicated as p increases. We point out that a similar geometric interpretation applies also to the set of non-trivial coefficients x m,n entering the solution for the symmetric scalar blocks in even number of dimensions.
The structure of the paper is as follows. In section 2 we briefly review the results of ref. [25] and define the CPWs and the CBs of interest. In section 3, we derive the system of p + 1 Casimir equations satisfied by the CBs G (p) e , for any p. This is reported in eq. (3.16). In section 4 we compute the CPWs using the shadow formalism approach. In particular, in subsection 4.1 we derive compact integral expressions of the CBs for any p and , eq.(4.25). In subsection 4.2 we write a more explicit expression of the CBs for = 0 and in subsection 4.3 we find another integral expression for the CPW, eq.(4.40), more suitable to perform computations with = 0. The solution of the Casimir system of equations is described in section 5. In subsection 5.1 we derive, by extending the results found in section 4, the asymptotic behaviour of the CBs. We discuss the form of the ansatz in subsection 5.2, and finally we reduce the coupled differential equations to an algebraic system in subsection 5.3. The solution of the CBs is finally derived in subsection 5.4, eq.(5.36). In subsection 5.5 we draw an analogy between the mixed tensor blocks G (p) e and the symmetric scalar blocks in d even dimensions. We conclude in section 6. Various technical details, as well as the explicit form of the coefficient defining the fermionic CBs entering scalar-fermion correlators (p = 1) are reported in two appendices.
The explicit form of all the coefficients c e m,n entering the CBs (5.36) for p = 1, 2, 3, 4 can be downloaded from https://sites.google.com/site/dskarateev/downloads.

Deconstructing Conformal Partial Waves
In 4D CFTs, for a given 4-point function, CBs and CPWs are labelled by the quantum numbers of the exchanged primary operator and thus they depend on its scaling dimension ∆ and rep-resentation ( ,¯ ) of the 4D Lorentz group, with and¯ positive integers. Four-point functions involving only scalar fields are the best known. In any channel, the exchanged operators havē = , i.e. they are all and only traceless symmetric tensors. In this case CPW and CB are equivalent up to a kinematic factor and their analytic form has been derived in a remarkable compact form in refs. [8,9] for any ∆ and . Four-point functions involving tensor (or fermion) operators are considerably more complicated because different tensor structures arise and more operators can be exchanged. A generic fermion-tensor 4-point function can be parametrized as where I i are schematic Lorentz indices of the operators O i (x i ), is a kinematic factor, and T I 1 I 2 I 3 I 4 n (x i ) are kinematically determined tensor structures. The dynamical information of the 4-point function is encoded in the N 4 functions g n (u, v). As we mentioned, a bootstrap analysis requires to rewrite the 4-point function (2.1) in terms of the operators exchanged in any channel. In the s-channel , for instance, we have are the associated CPWs. The sum over the exchanged primary operators O r includes a sum over all possible representations ( ,¯ ) that can appear in the 4-point function and, for each representation, a sum over all the possible primaries, i.e. a sum over all possible scaling dimensions ∆ Or . It is useful to define δ = |¯ − | and rearrange the sum over ( ,¯ ) in a sum over, say, and δ. There is an important difference between these two sums. For any 4-point function, the sum over l extends up to infinity, while the sum over δ is always finite. More precisely, we have δ = 0, 2 , . . . , p − 2, p, O r bosonic δ = 1, 3 , . . . , p − 2, p, O r fermionic. (2.5) In both cases, the integer p is defined to be and is automatically an even or odd integer when O r is a boson or a fermion operator. There are several CPWs for each exchanged primary operator O r , depending on the number of allowed 3-point function structures. They admit a parametrization like the 4-point function itself, where g Or,n (u, v) are the CBs, scalar functions of u and v that depend on the dimensions and spins of the external and exchanged operators. Imposing crossing symmetry by requiring the equality of different channels is the essence of the bootstrap approach. In order to successfully bootstrap the correlator (2.1), it is necessary to know the explicit form of the CPWs (2.7), in particular the CBs g (i,j) Or,n (u, v). It has been shown in ref. [25] that the CPWs associated to an operator O ( , +p) (and similarly for its conjugate O ( +p, ) ) exchanged in the OPE channel (12)(34) of a 4-point function can be obtained from a single CPW W seed O ( , +p) as follows: where D i 12 and D i 34 are differential operators that depend on O 1,2 and O 3,4 , respectively. For even integer p = 2n, the seed CPWs are those associated to 4-point functions of two scalar fields with one (2n, 0) and one (0, 2n) bosonic operators, while for odd integer p = 2n + 1, they consist of 4-point functions of two scalar fields with one (2n + 1, 0) and one (0, 2n + 1) fermionic operators: 2 (2.10) In the above correlators, in the OPE channel (12)(34) primary operators O ( , +δ) and their conjugates O ( +δ, ) can be exchanged only with the values of δ indicated in eq. (2.5) and any .
There are several 4-point functions in which the operators O ( , +p) and O ( +p, ) are exchanged and in which the corresponding CPWs have a unique structure. Among these, the correlators (2.9) and (2.10) are the ones with the minimum number of tensor structures and hence the simplest. This is understood by noticing that for any value of δ (and not only for δ = p) the operators O ( , +δ) and their conjugates O ( +δ, ) appear in both the (12) and (34) OPE's with one tensor structure only, since there is only one tensor structure in the corresponding three-point functions: . Despite this simplification, the above computation is still technically challenging. A further great simplification occurs by using the 6D embedding formalism [28][29][30][31] in twistor space with index-free notation [16]. As we will see, among other things, this formalism spare us from explicitly writing tensor structures with open indices for the correlators (2.9) and (2.10). The 4D conformal group is isomorphic to the 6D Lorentz group SO(4, 2), so by embedding the 4D fields φ(x) into 6D counterparts Φ(X), the non-linear conformal transformations turn into linear 6D Lorentz transformations. 6D fields are defined on a 4 dimensional subspace: the projective (X ∼ λX) light-cone (X M X M = 0) (see e.g. refs. [32,33] for further details). Using the local isomorphism SO(4, 2) ∼ SU (2, 2), 4D Weyl spinors ψ α (x) can be embedded either into twistors Ψ a (X) subject to a transversality constraint [32] or to twistors Ψ b (X) subject to a gauge redundancy [16]. Following refs. [16,26], we adopt here the latter possibility. A general 4D primary field Oβ 1 ...β¯ α 1 ...α (x), with scaling dimension ∆ in the ( ,¯ ) representation is embedded in a 6D multitwistor field O a 1 ...a b 1 ...b¯ (X), homogeneous in X with degree τ = ∆ + ( +¯ )/2. We can saturate the indices of multi-twistor fields by multiplying them by auxiliary twistors S and S to get index-free scalar quantities: The gauge redundancy requires that effectively where X ab = X M Σ M ab , X ab = X M Σ ab M , Σ M and Σ M are the 6D chiral gamma matrices. The lightcone condition requires also XX = 0 (for all definitions and more details see ref. [26]). All tensor structures in twistor correlators can be written in terms of scalar functions of auxiliary twistors S's and S's. For n-point correlators one can find a basis of all possible linearly independent functions of (S 1 , . . . , S n , S 1 , . . . , S n ); for n = 4 such basis includes, among others, the following invariants (i = j = k = l): where X ij = X M i X jM . An independent basis for the p + 1 tensor structures appearing in the 6D uplift of the correlators (2.9) and (2.10) can be obtained from the invariants in eq.(2.15): e on ∆ and . The CBs depend also on the external operator dimensions, more precisely on a and b, defined as (2.18) For simplicity of notation, we no longer distinguish between even and odd values of p, since we can consider both cases simultaneously. It is then understood that in the corrrelator (2.16) F are 6D uplifts of 4D fermion fields for p odd. It is possible to get W seed (p) from W seed (p), or vice versa, using the results of ref. [25] and a parity transformation P. We have where is the CPW associated to the parity dual 4-point function We do not report here the explicit form of the differential operators ∇ ij , D i , d 3 andd 1 , as well as the action of parity on them and on the 6D SU (2, 2) invariants, that can be found in ref. [25]. In fact, we will not use eq.(2.19) to compute W seed (p), because we will find an easier way to directly compute both W seed (p) and W seed (p).
Instead of eq.(2.16), we could have considered the alternative 4-point function to calculate an analogue seed CPW W seed (p). Since eq.(2.22) is equal to eq.(2.16) under the permutation 3 ↔ 4, the CBs appearing in the decomposition of W seed (p) and W seed (p) are related as follows: The 4D CPWs W seed O ( , +p) and W seed O ( +p, ) are obtained by projecting to 4D their 6D counterparts W seed (p) and W seed (p). There is no need to explicitly perform such projection, because the 4D CBs are directly identified with their 6D counterparts. One has simply

The System of Casimir Equations
In this section we derive the system of second order Casimir equations for the seed conformal blocks defined in eq. (2.17). Before addressing the more complicated case of interest, let us recall how the Casimir equation for scalar correlators is derived. One starts by considering the 4-point function whereĈ is the quadratic Casimir operator. 3 Recasting the generators of the 4D conformal group in a 6D form asL M N , with M, N 6D indices, we havê The Casimir equation is derived by expressing eq.(3.1) in two different ways. On one hand, we can replace in eq.(3.2) the operatorL M N with its explicit action in terms of differential operators acting on the scalar fields inserted at the points x 1 and On the other hand, we might consider the (12) OPE. Scalar operators can only exchange symmetric traceless operators, so p = 0 in this case, and one has where T is a tensor structure factor whose explicit form will not be needed. In the latter view, we end up having the commutator ofĈ with O ( , ) which gives the Casimir eigenvalue is the value associated to an operator in the ( + p, ) or ( , + p) Lorentz representations. Using then eq.(2.4) one derives a differential equation for each CPW, for any fixed ∆ and .
The explicit form of the second order differential operator acting on the CPW or directly on the CB is best derived in the 4 + 2-dimensional embedding space. The CPW of scalar correlators is parametrized by a single conformal block G (0) 0 (z,z). When acting on scalar operators at x 1 and x 2 , the Lorentz generator can be written as Plugging eq.(3.6) in eq.(3.2), one finds after a bit of algebra the Casimir equation [9] ∆ (a,b;0) 2 where a and b are defined in eq.(2.18), u = zz and v = (1−z)(1−z). The second-order differential operator ∆ is defined as in terms of the second-order holomorphic operator The above derivation can be generalized for CPWs entering 4-point correlators of tensor fields. As we have seen in section 2, in the most general case the exchange of a given field O ( ,¯ ) is not parametrized by a single CPW, but by a set of CPWs W (i,j) , whose number depends on the number of tensor structures defining the three-point functions (12O) and (34O). In order to derive the second order differential equation satisfied by W (i,j) one has to properly identify the OPE coefficients λ i appearing in the generalization of eq.(3.3) with those in eq.(2.4). This is not needed for the seed correlators (2.16) since the CPW is unique, like in the scalar correlator. For each p, we have where C is the explicit differential form of the Casimir operator to be determined and E p is as in eq.(3.5). An identical equation is satisfied by W seed (p). Contrary to the scalar case, the single differential equation (3.10) for W seed (p) turns into a system of equations for the p + 1 CBs G We can get rid of all the twistor indices by defining the index-free Lorentz generators Given any 6D tensor O(X, S,S) , we have 13) whereL M N satisfy the Lorentz algebra The explicit form of the Casimir differential operator entering eq.(3.10) is obtained by plugging eq.(3.12) in eq.(3.2). The single equation (3.10) for the CPW turns into a system of secondorder coupled differential equations for the p + 1 conformal blocks G   e (G) is conceptually straightforward but technically involved. The main complication arises from the spin part of the Lorentz generator (3.12) that generates products of SU (2, 2) invariants not present in eq.(2.17). The new invariants are linearly dependent and must be eliminated using relations among them. See Appendix A of ref. [25] for a list of such relations. This is a lengthy step, that however can be automatized in a computer. When redundant structures have been eliminated, one is finally able to read from eq.(3.15) the Casimir system Cas (p) e (G). Despite the complicacy of the computation, the final system of p + 1 equations can be written into the following remarkably compact form: 17) and the coefficients E p are given in eq.(3.5). In eq.(3.16) it is understood that G e . Interestingly enough, only two differential operators enter into the Casimir system: the second-order operator (3.8) that already features p = 0, with coefficients a e , b e and c e given by a e ≡ a, b e ≡ b + (p − e), c e ≡ p − e , (3.18) and the new linear operator L(µ) given by Another remarkable property of the Casimir system (3.16) is that, for each given e and p, at most three conformal blocks mix with each other in a sort of "nearest-neighbour interaction": G e mixes only with G e+1 and G e−1 . The Casimir equations at the "boundaries" Cas Finding the solution of the system (3.16) is a complicated task, that we address in the next sections.

Shadow Formalism
Another method to obtain CBs in closed analytical form uses the so called shadow formalism. It was first introduced by Ferrara, Gatto, Grillo, and Parisi [17][18][19][20] and used in ref. [8] to get closed form expressions for the scalar CBs. In this section we apply the shadow formalism, using the recent formulation given in ref. [16], to get compact expressions for W seed (p) and W seed (p) in an integral form for any p and . 4 Using these expressions, we compute the CBs G Despite the power of the above technique, it is computationally challenging to go beyond the p = 2 case. Moreover, as we will see, we do not have any control on the final analytic form of CBs. In light of this, we will provide the full analytic solution for G (p) e and G (p) e , for any p, only in section 5, where we solve directly the set of Casimir differential equations by using an educated ansatz for the solution. The results obtained in this section are however of essential help to argue the proper ansatz. They will also allow us to get the correct physical asymptotic behaviour of G

CPW in Shadow Formalism
We start by briefly reviewing the shadow formalism along the lines of ref. [16], where the reader can find more details. The CPW associated to the exchange of a given operator O r with spin ( ,¯ ) in a correlator of four operators O n (X n ), n = 1, 2, 3, 4 (in embedding space and twistor language) is given by where ν is a normalization factor, the projector gluing two 3-point functions is given by and O r is the shadow operator Or(Y, YS,Ȳ S) .
The integral in eq.(4.1) would actually determine the CPW associated to the operator O r (X, S,S) plus its unwanted shadow counterpart, that corresponds to the exchange of a similar operator but with the scaling dimension ∆ → 4 − ∆. The two contributions can be distinguished by their different behaviour under the monodromy transformation X 12 → e 4πi X 12 . In particular, the physical CPW should transform with the phase e 2iπ(∆−∆ 1 −∆ 2 ) , independently of the Lorentz quantum numbers of the external and exchanged operators. This projection on the correct monodromy component explains the subscript M in the bar at the end of eq.(4.1). We use eq.(4.1) to get an integral form of W seed (p) and W seed (p) in eq.(2.17). The explicit expressions of the needed 3-point functions are given by is a kinematic factor and are SU (2, 2) invariants for three-point functions. The "shadow" 3-point function counterparts are given by Using the above relations, after a bit of algebra, one can write where and We will not need to determine the normalization factors ν andν in eqs. (4.8) and (4.9). Notice that the correct behaviour of the seed CPWs under X 12 → e 4πi X 12 is saturated by the factor X 12 multiplying the integrals in eqs. (4.8) and (4.9). Hence the latter should be projected to their trivial monodromy components M = 1, as indicated. Notice that eqs.(4.11) and (4.12) are related by a simple transformation: where P is the parity operator. We can recast the expression (4.11) in a compact and convenient form using some manipulations. We first define 3 variables Then we look for a relation expressing the generic N (p) in terms of the known N (0): where C p are Gegenbauer polynomials of rank p. Starting from eq.(4.11), after acting with the S and T derivatives, one gets where we have defined Ω ab = (X 2X1 X 0X3 X 4 ) ab . In order to relate N (p) above to N +p (0) in eq.(4.15), we look for an operator D satisfying We deduce that D should be bilinear inS 4 and S 2 and should commute with ( − → ∂SX 0 − → ∂T ). In addition to that, it should have the correct scaling in X's and should be gauge invariant, namely it should be well defined on the light-cone X 2 = 0 and preserve the conditions (2.14). It is not difficult to see that the choice D = D/(8X 01 X 04 ), where fulfills all the requirements. One has D(SΩT ) = (SS 2 )(S 4 X 3T ). Iterating it p times gives the desired relation: The operator D annihilates all the scalar products with the exception of X 12 , in which case we have DX 12 = I 2 , and we define The action on the s, t, and u variables is on Gegenbauer polynomials is and vanishes on J 42,01 and J 42,30 . Using recursively the identity for Gegenbauer polynomials we can write the following expression for N (p): where p w is the binomial coefficient and for compactness we have defined the dimensionful tensor structures Combining together eqs.(4.8), (4.9), (4.13), (4.14) and (4.24) we can finally write where ν andν are undetermined normalization factors.

Seed Conformal Blocks and Their Explicit Form for = 0
The computation of the CBs G (p) e and G (p) e starting form eq.(4.25) is a non-trivial task for generic and p, since we are not aware of a general formula for an integral that involves C p+1 −w (t) for p = 0. For any given , one can however expand the Gegenbauer polynomial, in which case the CBs G Recalling the definition of t in eq.(4.14), one realizes that the Gegenbauer polynomials in eq.(4.25), when expanded, do not give rise to intrinsically new integrals but just amounts to shifting the exponents in the denominator. The tensor structures in the numerators bring p open indices in the form X N 1 0 . . . X Np 0 , which can be removed by using eq.(3.21) in ref. [16]. In this way the problem is reduced to the computation of scalar integrals in 2h = 2(2 + p) effective dimensions, of the form: where A 01 +A 02 +A 03 +A 04 = 2h. The capital A 0i are used for the exponents in the denomentaor with all possible shifts introduced by the Gegenbaur polynomials. This integral is given by The derivative −∂/∂v in (z,z) coordinates equals In the case of = 0, all the above manipulations simplify drastically. The Gegenbauer polynomials C p+1 −w (t) vanishe for all the values w except for w = 0, leaving only one type of tensor structure: I p 1 for W seed (p) and I p 2 for W seed (p). This leads to a one-to-one correspondence between CBs and integrals: (4.32) The asymptotic behaviour of the CBs when z,z → 0 (u → 0, v → 1) for = 0 is easily obtained from eq.(4.31) by noticing that R (h) (z,z; A 02 , A 03 , A 04 ) is constant in this limit. Then we have lim z→0,z→0 By knowing that the CBs should be proportional to the factor in eq.(4.32), we can refine eq.(4.33) and write lim z→0,z→0 Notice that the behavior (4.34) and (4.35) of the CBs for z,z → 0 when = 0 is not guaranteed to be straightforwardly extended for any = 0. Indeed, we see from eq.(4.25) that for a given p, the generic CPW is obtained when ≥ p, in which case all terms in the sum over w are present.
All the values of < p should be treated separately.

Computing the Conformal Blocks for = 0
A useful expression of the CBs for generic values of can be obtained using eq.(4.19) and the known closed form of W seed (0). Recall that where a and b are as in eq.(2.18) for p = 0 and G (0) (z,z) are the known scalar CBs [8,9] expressed in terms of the function 5 k (a,b;c) ρ (z) ≡ z ρ 2 F 1 (a + ρ, b + ρ; c + 2ρ; z) . (4.38) 5 We adopt here the notation first used in ref. [3] for this function, but notice the slight difference in the definition: k there ρ = k here ρ/2 .
Comparing eq.(4.36) with eq.(4.25) for p = 0, one can extract the value of the shadow integral in closed form for generic spin [16]: where D = PD| 1↔3,2↔4 , as follows from eq.(4.13), The tensor integral is evaluated using SO(4, 2) Lorentz symmetry. One writes where n runs over all possible rank p traceless symmetric tensors τ n which can be constructed from X 1 , X 2 , X 3 , X 4 and η M N 's, with arbitrary scalar coefficients A n to be determined. Performing all possible contractions, which do not change the monodromy of the integrals, the A n coefficients can be solved as linear combinations of the scalar block integrals I defined in eq.(4.39), with shifted external dimensions. In this way, we have computed the CBs G p−e , for any p. More on this point in section 5. In all cases the CBs satisfy the Casimir system (3.16).
As mentioned at the end of subsection 4.2, the asymptotic behaviour of the CBs for z,z → 0 depends on whether ≥ p or not. For p = 1 we can expand the obtained solutions, which for ≥ 1 read as

Solving the System of Casimir Equations
The goal of this section is to find the explicit form of the conformal blocks G (p) e and G (p) e appearing in eq.(2.17) by solving the Casimir system (3.16). In doing it we adopt and expand the methods introduced by Dolan and Osborn in refs. [9,27] to obtain 6D scalar conformal blocks. We will mostly focus on the blocks G (p) e , since the same analysis will apply to G (p) e with a few modifications that we will point out.
Before jumping into details let us outline the main logical steps of our derivation. We first find, with the guidance of the results obtained in section 4, the behaviour of G (p) e and G (p) e in the limit z,z → 0 in which the Casimir system (3.16) can be easily solved. Using this information and eq.(4.32), we then write an educated ansatz for the form of the CBs. Using this ansatz, we reduce the problem of solving a system of linear partial differential equations of second order in two variables to a system of linear algebraic equations for the unknown coefficients entering the ansatz. Then we show that the non-zero coefficients in the ansatz admit a geometric interpretation. They form a two-dimensional lattice with an octagon shape structure. This interpretation allows us to precisely predict which coefficients enter in our ansatz for any value of p. Finally, we show that the linear algebraic system admits a recursive solution and we discuss the complexity of deriving full solutions for higher values of p.

Asymptotic Behaviour
Not all solutions of the Casimir system (3.16) give rise to sensible CBs. The physical CBs are obtained by demanding the correct boundary values for G where N e , λ (e) andλ (e) are parameters to be determined. For simplicity of notation we have omitted their p-dependence. The differential operators (3.8) and (3.19), when acting on eq.(5.1) give, at leading order in z andz, ∆ (ae,be;ce) → λ (e) (λ (e) − 1) + c e (λ (e) +λ (e) ) +λ (e) (λ (e) − 1) − λ (e) , (5.2) e with e = p. In the limit z,z → 0 it reads For generic values of , we have λ (e) =λ (e) . Hence we cannot have λ (p−1) + 1 < λ (p) in eq.(5.4), since this would imply that the last term dominates in the limit and N p−1 vanishes, in contradiction with the initial hypothesis (5.1). Let us first consider the case in which λ (p−1) + 1 > λ (p) , so that the terms in the second row of eq.(5.4), coming from G (p) p−1 , vanish. It is immediate to see that the only sensible solution for λ (p) andλ (p) that reproduce the known OPE limit for the p = 0 case is Notice that eq.(5.5) agrees with the asymptotic behaviour for the CBs G According to eq.(4.42), we expect λ (p−2) = λ (p−1) = λ (p) ,λ p−1 =λ (p) − 1,λ p−2 =λ (p) − 2 in eq.(5.6). In this case the last term is higher order in z and eq.(5.6) is satisfied by simply taking .
Notice that we have tacitly assumed above that λ (p) −λ (p) = − does not vanish, i.e. = 0. For = 0, more care is required and one should consider the first subleading term inz in the expansion (5.1). The above analysis can be iteratively repeated until the last equation Cas where λ (p) andλ (p) are as in eq.(5.5) and e = 0, . . . , p − 1. The asymptotic behaviour of the CBs in the OPE limit is given for any and p by lim z→0,z→0 We do not report the explicit form of the normalization factors N e , since they will be of no use in what follows. We still have to consider the case in which λ (p−1) + 1 = λ (p) in eq.(5.4). By looking at eq.(4.43), it is clear that this case corresponds to the asymptotic behaviour of the conjugate CBs G where The asymptotic behaviour of the conjugate CBs are given for any and p by lim z→0,z→0

The Ansatz
The key ingredient of the ansatz is the function k (5.14) Using eq.(5.14) one can define an eigenfunction of the operator ∆ (a,b;c) 0 as the product of two k's: These functions played an important role in ref. [9] for the derivation of an analytic closed expression of the scalar CBs in even space-time dimensions. In our case, the situation is much more complicated, because we have different blocks appearing in the Casimir equations. We notice, however, that the second order operator ∆ in each equation Cas e+1 are multiplied by first order operators only. Since, as we will shortly see, first order derivatives and factors of z andz acting on the functions F can always be expressed in terms of functions F with shifted parameters, a reasonable ansatz for the CBs is to take each G e proportional to a sum of functions of the kind F (ae,be;ce) ρ 1 , ρ 2 (z,z) for some ρ 1 and ρ 2 . Taking also into account eq.(4.32), found using the shadow formalism, the form of the ansatz for the blocks G where c e m,n are coefficients to be determined and the sum over the two integers m and n in eq.(5.17) is so far unspecified. Notice that all the functions F entering the sum over m and n have the same values of a e , b e and c e . Matching eq.(5.17) in the limit z,z → 0 with eq.(5.9) allows us to determine ρ 1 and ρ 2 , modulo a shift by an integer. We take 18) in which case the sum over n is bounded from below by n min = −p. At this value of n, we have m(n min ) = e − p. There is no need to discuss separately the behaviour of the blocks with ≤ p. Their form is still included in the ansatz (5.17) with the additional requirement that some coefficients c e m,n min should vanish. This condition is automatically satisfied in the final solution. In the next subsections we will discuss the precise range of the sum over m and n and explain how the coefficients c e m,n can be determined.

Reduction to a Linear System
The eigenfunctions F ± (a,b;c) ρ 1 , ρ 2 (z,z) have several properties that would allow us to find a solution to the system (3.16). In order to exploit such properties, we first have to express the system (3.16) for G  e (z,z) defined in eq.(5.17). We plug the ansatz (5.17) 6 Recall that the conformal blocks are even under z ↔z exchange, that leaves u and v unchanged. in eq.(3.16) and use the following relations to obtain the system of Casimir equations for g (p) e : We have split each Casimir equation in terms of three differential operators Cas 0 , Cas + , Cas − , that act on g e−1 , respectively. In order to avoid cluttering, we have omitted the obvious e and p dependences of such operators. Their explicit form is as follows: where A 0 , A − and A + are coefficients that in general depend on all the parameters involved: a, b, ∆, , e and p but not on z andz, namely they are just constants. For future purposes, in eqs.(5.24)-(5.26) we have only made explicit the dependence of A 0 , A − and A + on the integers m and n. The sum over (r, t) in each of the above terms runs over a given set of pairs of integers. We report in fig. 1 the values of (r, t) spanned in each of the three regions R 0 , R + and R − . We do not report the explicit and quite lengthy expression of the coefficients A 0 r,t , A + r,t and A − r,t , but we refer the reader again to appendix A where we provide all the necessary relations needed to derive them. Using eqs. (5.27) The functions F − appearing in eq.(5.27) are linearly independent among each other, since they all have a different asymptotic behaviour as z,z → 0. Hence the only way to satisfy eq.(5.27) is to demand that terms multiplying different F − vanish on their own: where m = m + r, n = n + t. The Casimir system is then reduced to the over-determined linear algebraic system of equations (5.28).

Solution of the System
In order to solve the system (5.28), we have to determine the range of values of (m, n) entering the ansatz (5.17), that also determines the size of the linear system. By rewriting the known p = 1 and p = 2 CBs found using the shadow formalism in the form of eq.(5.17), we have deduced the range in (m, n) of the coefficients c e m,n for any p (a posteriori proved using the results below). For each value of e, the non-trivial coefficients c e m,n span a two-dimensional lattice in the (m, n) plane. For each e, the shape of the lattice is an octagon, with p and e dependent edges. The position and shape of the generic octagon in the (m, n) plane is depicted in fig. 2. One has n min = − p, n max = e + p, m min = e − 2 p, m max = p . The system (5.28) is always over-determined, since it is spanned by the values (m , n ) whose range is bigger than the range of (m, n) ∈ Oct (p) e (spanning all the coefficients to be determined) due to the presence of (r, t) ∈ [−2, 2]. There are only N p − 1 linearly independent equations, because the system of Casimir equations can only determine conformal blocks up to an overall factor. The most important property of the system (5.28) is the following: while the number of equations grows with p, the total number of coefficients c e m,n entering any given equation in the system (5.28) does not. This is due to the "local nearest-neighbour" nature of the interaction between the blocks, for which at most three conformal blocks can enter the Casimir system This factor is generally non-vanishing, unless m = 0 and n = −p, in which case it vanishes for any ∆, and p. In this way eq.(5.33) selects c p 0,−p as the only non-vanishing coefficient at level n = −p for e = p. Notice that it is crucial that A 0 0,−2 (m, n)| e=p vanishes automatically for a given pair (m, n), otherwise either the whole set of equations would only admit the trivial solution c e m,n = 0, or the system would be infinite dimensional. A similar reasoning applies for the other three coefficients. One has in particular that are automatically satisfied because the three coefficients A 0 0,2 , A 0 2,0 and A 0 −2,0 vanish when evaluated for the specific values reported in eq.(5.34) for any ∆, and p.
The system (5.28) is efficiently solved by extracting a subset of N p − 1 linearly independent equations. This can be done by fixing the values (r, t) = (r * , t * ) entering the definitions of (m , n ). There are 4 very special subsets of the N p − 1 equations (corresponding to very specific values (r * , t * )) which allows us to determine the solution iteratively starting from eq.(5.28). They correspond to a solution where one of the four coefficients (5.32) is left undetermined, in other words (r * , t * ) can be set to be (0, −2), (0, 2), (2, 0) or (−2, 0). For instance, if we choose c 0 ≡ c p 0,−p as the undetermined coefficient, a recursion relation is found from eq.(5.28) by just singling out the term with t = −2 in A 0 and setting (r * , t * ) = (0, −2). Such a choice leads to m = m, n = n − 2, and one finally gets It is understood in eq. It is important to remind the reader that the CBs G  blocks, but no significant extra complications come from the external operators. An example of such phenomenon in a four fermion correlator is shown (though in a less significative way) in section 7.1 of ref. [25]. For any given 4-point function, after the use of the differential operators introduced in ref. [25], there is no need to compute the coefficients c e m,n for any a and b but only for the values of interest. This considerably simplifies the expression of c e m,n .

Analogy with Scalar Conformal Blocks in Even Dimensions
It is worth pointing out in more detail some similarities between the CBs G (p) e for mixed symmetry tensors computed above and the scalar conformal blocks G d in d > 2 even space-time dimensions (G 4 = G (0) 0 in our previous notation). The quadratic Casimir equation for scalar CBs in any number of dimensions is is the quadratic Casimir eigenvalue for traceless symmetric tensors. The explicit analytical form of scalar blocks in d = 2, 4, 6 dimensions has been found in refs. [8,9]. The same authors also found a relation between scalar blocks in any even space-time dimensionality, eq.(5.4) of ref. [9] (see also the more elegant eq.(4.36) of ref. [27]), that allows us to iteratively determine G d for . These five points form a slanted square in the (m, n) plane, centered at the origin. The explicit form of the coefficients x m,n is known, but it will not be needed in what follows. 8 It is natural to expect that eq.(5.41) should apply for any even d ≥ 4, with a number of nonvanishing coefficients that increases with d. 9 This is not difficult to prove. From the first relation in eq.(5.19) we can get the form of the Casimir equation for the function g d (z,z) defined in eq.(5.41), that can be written as Using the techniques explained in subsection 5.3 and the results of appendix A, it is now straightforward to identify which is the range of (m, n) of the non-vanishing coefficients x m,n for any d (see fig.4). 10 In d dimensions, the minimum and maximum values of m and n are given by The number N d of coefficients x m,n entering the ansatz (5.41) for scalar blocks in d even spacetime dimensions is easily computed by counting the number of lattice points enclosed in the slanted square. We have e , in analogy to the operator (4.35) of ref. [27] relating G d+2 to G d . It would be very useful to find, in this or some other way, a more compact expression for the blocks G On the contrary, no such simple guess seems to be possible for the power 2p + 1 entering G (p) e , given also the appearance of the operator L defined in eq.(5.19). As discussed, we have fixed the power 2p + 1 by means of the shadow formalism.

Summary and Conclusions
We have computed in this paper the seed CBs G e for any e, p, ∆, and external scaling dimensions, by solving the Casimir set of differential equations, that can be written in the compact form (3.16). The shadow formalism has been of fundamental assistance to deduce it and also as a useful cross check for the validity of the results. The final expression for the CBs is given in eq.(5.36), the most important formula in the paper. The CBs are expressed in terms of coefficients c e m,n , that can be determined recursively, e.g. by means of eq.(5.35). For each CB, the coefficients c e m,n span a 2D octagon-shape lattice in the (m, n) plane, with sizes that depend on p and e and increase as p increases. We have reported in Appendix B the explicit form of c e m,n for the simplest case p = 1. We have not reported the c e m,n for higher values of p, since their number and complexity grows with p. Their explicit form up to p = 4 can be downloaded from a website. The CBs up to p = 4 are enough to bootstrap many tensor correlators, including four conserved spin 1 currents. Aside from the obvious application in the numerical bootstrap, the knowledge of the CBs G (p) e should be useful in other contexts. Having analytical closed expressions for the blocks should be very useful to generalize the so called analytic bootstrap approach [36,37] to tensor correlators. It would also be interesting to explore holographic interpretations in AdS 5 of the CBs G (p) e and their possible utility in the formulation of higher spin theory in AdS 5 . It should be interesting, for bootstrap applications, to systematically study the correlation functions associated to the free theories of (p, 0) spinor/tensor fields, that admit a particle interpretation [38]. Indeed, aside from the interest per se, this study might give a useful analytically known benckmark point for future bootstrap analysis involving the CBs G (p) e found in this paper, very much like the role that free scalar theories play in actual analytical and numerical studies.
The somewhat surprisingly simple form of the Casimir system (3.16), where at most three blocks at a time can enter in a sort of local interaction, and the geometric interpretation of the coefficients c e m,n in terms of octagons, are perhaps an indication of a more fundamental symmetry principle. This should hopefully allows us to gain a better understanding of 4D CFTs or at least, less ambitiously, more compact expressions for the CBs G (p) e . Even in absence of a would be underlying symmetry, it is well possible that there is a better way to parametrize the blocks that we might have overlooked.
Acknowledgments A.C.E. thanks Emilio Trevisani and Mauro Valli for useful discussions. The work of M.S. was supported by the ERC Advanced Grant no. 267985 DaMESyFla. A.C.E. and M.S. gratefully acknowledge support from the Weizmann Institute of Science, and the Simons Center for Geometry and Physics, Stony Brook University at which some of the research for this paper was performed.

A Properties of the F Functions
In this Appendix we provide all the properties of the functions F (a,b;c) ρ 1 , ρ 2 needed for the system of Casimir equations and more specifically to derive eqs.(5.24)-(5.26). We will not consider the functions F ± (a,b;c) ρ 1 , ρ 2 here, since their properties can trivially be deduced from the ones below by demanding both sides to be symmetric/anti-symmetric under the exchange z ↔z.
The fundamental identities to be considered can be divided in two sets, depending on whether

B The Conformal Blocks for p = 1
We report in this appendix the full explicit solution for the two conformal blocks G