Scalar-Fermion Analytic Bootstrap in 4D

In this work we discuss an analytic bootstrap approach [1,2] in the context of spinning 4D conformal blocks [3,4]. As an example we study the simplest spinning case, the scalar-fermion correlator $\langle\phi\psi\phi\bar\psi\rangle$. We find that to every pair of primary scalar $\phi$ and fermion $\psi$ correspond two infinite towers of fermionic large spin primary operators. We compute their twists and products of OPE coefficients using both s-t and u-t bootstrap equations to the leading and sub-leading orders. We find that the leading order is represented by the scalar-fermion generalized free theory and the sub-leading order is governed by the minimal twist bosonic (light scalars, currents and the energy-momentum tensor) and fermionic (light fermions and the suppersymmetric current) operators present in the spectrum.

The obtained CFT result can be interpreted in AdS as a system of scalar particles with large relative angular momentum . One can compute shifts of their binding energies due to gravity [82][83][84][85] and other interraction [86]. This matches nicely the sub-leading correction in (1.1) governed by the energy-momentum tensor on the CFT side.
Extending the analytic techniques to correlators with spin is crucial for studying generic CFT data. The works [87,88] addressed correlators with two scalars and two traceless symmetric operators 4 such as conserved currents J µ and the energy-momentum tensor T µν . This has allowed to prove [88] the conjectured conformal collider bounds [89]. The work [87] deals solely with 3D CFTs, whereas in [88] the authors work in d ≥ 3 dimensions. Analytic studies of fermions also got some attention recently [90,91].
Goal of the paper In this work we will discuss the analytic method of [1,2] in the context of generic spinning correlators in 4D [3,4,31,32]. 5 We will follow the notation of [32]. For performing computations in practice we will develop and attach to the paper a little Mathematica [93] code "analyticBootstrap4D.m", which should be used together with the "CFTs4D" package [32]. We mention the names of relevant functions throughout the text as [functionName].
As a demonstration in this paper we will study in detail the simplest spinning correlator in 4D, namely the scalar-fermion correlator 6 φ ψ φ ψ . (1.3) We will address both s-t and u-t channels and discuss how they complement each other. 7 To facilitate the discussion we will also study the spectrum and the OPE coefficients of the scaclar-fermion GFT. We will discuss how the GFT results are related to the large-spin analytic results. We will provide details of all the computations in the attached Mathematica notebook "Example.nb".
Outline of the paper In section 2 we discuss in details generic spinning 4-point functions.
We provide a recipe for studying them analytically and develop practical tools for performing computations in Mathematica. In section 3 we consider a demonstration of the strategy given in section 2 on the scalar-fermion 4-point function. A compact summary of the final result is given in section 3.2.3. In section 3.3 we address a scalar-fermion GFT and show how it is related to our analytic bootstrap results. We conclude in section 4. Details on the Ward identities for the scalar-fermion correlator are provided in appendix A.

Analytic Bootstrap
We start by discussing the main ingredients of general bootstrap equations in section 2.1 preparing the ground for section 2.2, where we provide an analytic bootstrap recipe in a schematic form.

Bootstrap Equations
We work in 4D Minkowski space where the conformal group is SO(2, 4) SU (2, 2). The local primary operators transform in a finite-dimensional representation of the U (1) × SU (2) × SU (2) sub-group and thus are labeled by the U (1) charge (the scaling dimension) ∆ and the pair of integers ( ,¯ ) which describes spin. 8 In what follows it will be more convenient to use instead of the scaling dimension another quantity called the twist defined as The inequalities in (2.1) are the unitary bounds, they are saturated by conserved operators [94,95]. In the ¯ = 0 case the conserved operators are free [96]. A generic local operator is denoted by where x is the position and s ands are auxiliary spinors (polarizations). The additional index k is used to label degenerate operators, which arise for example in the case of global symmetries or Dirac fermions. In the reminder of this paper we drop this label for simplicity. We use a convention for naming operators (2.2) such that Every operator (2.2) has its hermitian conjugate (2.4) 8 Traceless symmetric operators in 4D are those with =¯ .
A generic 4-point function of local operators has the following form when expanded in the s-channel where λ a are the OPE coefficients (3-point coupling constants) and W ab are the conformal partial waves (CPWs) which can be written as Any conformal block H ab I (z, z) can be expressed [3] in terms of the derivatives of simpler blocks called the seed blocks and dual seed blocks. Both seed and dual seed blocks have the following form [4] The coefficients c m,n are some constants [4,32] which have the following schematic form c e m,n = (−1) i p × rationalF unction m,n (a, b; , τ ). (2.9) An example of how the rationalF unction looks like can be found in appendix B in [4] for the p = 1 case. The functions F are defined as Here the parameters a and b depend on the twists of external operators, parameter c on the spin difference p defined in (2.8). The parameters ρ 1 and ρ 2 depend on twist and spin of the internal (exchanged) operator as where r and s are some integer or half-integer constants. See [4] for details. The k-functions are defined as For b = c − a the k-functions admit the following representation as a double-infinite sum (2.13) which simply follows from Taylor expansion around x = 0. The Pochhammer symbol and the digamma functions are denoted by respectively. Another important representation of (2.12) can be obtained in the large >> 1 limit and small x << 1, which reads as where K α is the modified Bessel function and the parameter ρ 1 is defined in (2.11). 9 This relation is a straightforward generalization of the results of appendix A.1.4 in [1]. For a given correlator one can write essentially 2 bootstrap equations: the s-t equation and the u-t equation

s-t bootstrap equation
We consider now the bootstrap equation (2.16). After decomposing it into a set of independent equations by stripping of independent tensor structures T I , all the equations get the following schematic form It can be studied in the z, z ∈ [0, 1] region. Here P O and P O represent products of the OPE coefficients We will study this equation in the z, z ∈ [1, +∞] region. We make a change of variables from (z,z) to (ω, ω) defined as Even though the hypergeometric series defining the k(ω) functions (2.12) is convergent in the ω ∈ [0, 1] region, the k(ω) functions is analytic in the whole complex plane modulo branch cuts [97], see figure 2.1. The following formula holds [98] k (a,b;c) The factor (−1) ±ρ appears due to the choice of analytic continuation ±iπ , namely one can choose to be right above or right below the left brunch cut, as indicated in figure 2.1 by arrows. We perform an analytic continuation of the functions k(ω) and k(ω) separately choosing opposite directions. It is required when moving to Euclidean space where ω = ω * . As a consequence we get where according to (2.11) ρ 1 − ρ 2 is an integer and one has After these manipulations instead of equation (2.20) we can write Structurally it has exactly the same form as (2.18) but with an extra (−1) . There is an equivalent and somewhat simpler way of dealing with the u-t bootstrap equations. One can consider a different correlation function (with positions of operators re-ordered) and to study an s-t channel (2.26) In this case the standard s-t techniques straightforwardly apply and no analytic continuation is needed.

The Analytic Bootstrap Recipe
In Minkowski space the variables z and z are independent quantities. The light-cone limit is defined as 10 We summarize here the recipe of [1] for studying analytically the bootstrap equations of the form (2.18) in the light-cone limit (2.27). We keep the discussion at the schematic level and pretend that the conformal blocks H are simply represented by a single F function (2.10) together with the (−1) i p factor coming from (2.9) 11 The left-hand side of (2.18) in the light-cone limit reads as 12 [KeepLeadingTerm] The first term comes from the identity operator that may or may not be present. This depends on the correlator under consideration. The second term describes the contribution of all the rest of operators in the theory. In the z → 0 limit only the lowest (minimal) twist operators contribute (light scalars or fermions with 1 < τ ≤ 2 and conserved bosonic or fermionic operators with τ = 2). We then recast the expression (2.29) into the final form by performing a series expansion around small 1 − z [kFunctionExpand]. 13 10 In the ω-variable we will use the light-cone limit ω 1 − ω 1. 11 We thus ignore the necessity of taking a linear combination in F with different coefficients cm,n and taking derivatives in z and z in general situation. We also drop the factor zz z−z 2 p+1 which in the limit (2.27) simply leads to (−z) 2 p+1 . 12 In writing (2.29) we assume ρ1 < ρ2. Depending on parameters r and s defined in (2.11) there might be a situation when ρ1 > ρ2. The second term in (2.29) should change accordingly. 13 The action of [kFunctionExpand] on k(z) is equivalent to (2.13) for x = 1 − z and b = c − a.
The right-hand side of (2.18) will be dominated (in the light-cone limit) by the large spin operators as was argued in [1,2]. 14 Using the definition (2.11) it is clear that the second term in (2.10) is negligible since it behaves as (1 − z) and should be dropped using [plugSeedBlocks[0], plugDualSeedBlocks[0]]. We approximate the sum over by an integral where the distribution density R τ ( ) indicates how operators with a given twist τ contribute to the sum. The right-hand side of (2.18) is then written as 15 [largeSpinExpand] In the large limit all the coefficients (2.9) simplify significantly but remain non-zero for generic parameters. They do not give any contribution apart from the (−1) factor. This is the reason we can safely leave them out of the discussion in (2.28).
In the large spin limit we expect the following behavior of the squared OPE coefficients multiplied by the density where the coefficients P 0 and P 1 contain the dependence on only via the factor (−1) which may or may not be present. In what follows we assume R τ ( ) = 1. 16 Presence of the 4 − factor in (2.32) guaranties sufficiently fast decay of the OPE coefficients required by the OPE convergence [99,100] and should cancel out the 4 factor in (2.31). The same type of ansatz was also used in [87,88]. Using (2.32) one can then evaluate the integral in (2.31) analytically term by term It is clear that, if the factor (−1) is present in the final expression, the integral would return zero, since contributions from even and odd spins are approximately the same and cancel each other out. The twist has the following behavior where τ 0 is the leading twits and γ is the coefficient of the sub-leading correction often called the anomalous dimension. Plugging (2.34) to (2.31) one should make an expansion around τ 0 for small γ D of the expression We then recast the expression (2.31) into the final form by performing a series expansion around small 1 − z [kFunctionExpand]. For example at the leading order in 1 − z one has Matching the left-and right-hand side After performing all prescribed manipulations one can equating (2.29) and (2.31) and determine the unknown coefficients in (2.32) and (2.34) by matching independently terms with equal powers where # 3 can only be 0 or 1. More precisely we start analyzing terms with # 3 = 0. Matching the power # 1 allows to determine the exponent B and matching the power # 2 allows to determine the leading twist τ 0 . We have a series expansion in 1 − z, to satisfy the equations order by order in 1 − z it is required to have an infinite tower of operators with the twists τ 0 which differ by 2n, where n is a non-negative integer. We can then determine the coefficient P 0 . By studying the terms with # 3 = 1 one gets access to the anomalous dimension. After extracting the coefficient γ and the exponent D, one can also compute P 1 .
In the case of bootstrap equations with no identity operator we do not get access to all the parameters, we will see an instance of this in section 3.2.2.

Structure of the Correlator
The correlator (3.1) has the following structure where g I (z,z) is an arbitrary function of conformally invariant variables and T I are the tensor structures defined as The objectsÎ 42 andÎ 42 31 represent the spin structure of the correlator, their explicit form is irrelevant for current discussion (see [32] for precise definitions).
We address now the s-, t-and u-channel OPE decompositions of the correlator (3.2). For brevity we will use the following names for bosonic (traceless symmetric) and fermionic operators consistent with (2.3) (

3.4)
It is convenient to write explicitly the contribution of operators and their hermitian conjugates to the sum (2.7). The decomposition (2.5) reads then The labels seed and dual seed indicate that the CPWs have the simplest possible structure among the spinning blocks in 4D, see (2.8). Here we work with p = 1 (dual) seeds. In the last equation the CPWs W a can be computed from the seed CPW with p = 0. We provide below all the relevant information on how to do it.
The products of OPE coefficients λ are generically denoted by P . The OPE coefficients λ are defined via the basis of tensor structures in a given 3-point function. We choose the following basis of structures for correlators involving Ψ and Ψ and for correlators involving O The object K 3 ,Î ij ,Ĵ i jk ,K jk i andK jk i encode spin and scaling structure of correlators. Their explicit form is irrelevant for current discussion (see [32] for precise definitions).
The four OPE coefficients defined via (3.8)-(3.11) are not independent of each other. Using complex conjugation of the defining correlators and permuting operators one gets the following relations 18 (3.14) Further constraints on the scaling dimensions and the OPE coefficients of external operators in (3.8)-(3.13) appear if operators Ψ and O are conserved. In case O is conserved, the correlators (3.12), (3.13) automatically satisfy the differential conservation equations. The only further constraint comes from the Ward identities discussed in appendix A which allow to relate the OPE coefficients to normalization constant of 2-point functions. We use the normalization convention of [32]. In the case of the energy-momentum tensor T we have In case Ψ is conserved, we deal with supersymmetric current. Conservation condition implies that in the correlators (3.8), (3.11) and in (3.9), (3.10) the twists of external operators should be related as τ φ = τ ψ and τ φ = τ ψ + 1 (3.16) 18 Since the products of OPE coefficients form positive (negative) definite quantities one can try to study the s-t bootstrap equations numerically using standard techniques. However since s-and t-channels are not unitary (or reflection positive) one does not expect to have any bounds. This is because (as we will see later) for such configuration the correlation function has an alternating factor (−1) which prevents one from constructing a positive definite functional.
respectively. This is expected in SUSY theories where scalars and fermions belong to the same multiplets. In principle we could have also used Ward identities to compute the OPE coefficients in (3.8)- (3.11) in the presence of SUSY current, but this is beyond the scope of our work. Finally we can construct tensor structures in (3.13) from the single scalar (seed) tensor structure by means of differential operators [26] where (•) stresses that we do not deal with the full correlator but only with its kinematic part (tensor structure), the differential basis of operators D and the matrix of coefficients M can be chosen as follows .

Analytic Results
We study the u-t and s-t bootstrap equation for the correlator (3.1). We perform their analytic analysis to the leading and sub-leading orders. The sub-leading corrections will depend on the OPE coefficients of the minimal (lowest) twist operators. Let us discuss in more details what minimal twist operator do we expect.
In the u-channel (3.7) we can exchange only traceless symmetric bosons, thus according to (2.1) the minimal twist operators can be light scalars, conserved = 1 currents and the energy-momentum tensor. We will denote twist and spin of these operator by τ Om and Om . In the s-channel (3.5) we can exchange only fermionic operators, according to (2.1) they could be light (1, 0) fermions or the super-symmetric current (2, 1). We will denote twist and spin of these operators by τ Ψm and Ψm .
We use the ansatz (2.34) and (2.32) for twists and products of OPE coefficients respectively. We report these expressions here again for convenience

Equations: u-t-channel
We equate (3.7) with (3.6) and perform a change of variables according to the prescription below (2.21). Schematically the bootstrap equation has the form (2.25) and can be studied in the light-cone limit Applying the procedure of section 2.2 we have found that there are two towers of large spin fermionic operators Ψ : X ( +1, ) , Y ( +1, ) (3.22) together with their hermitian conjugates X and Y . We determine their twists together with the products of OPE coefficients To the leading order twists read as where n is a non-negative integer spanning the tower of operators. Both anomalous dimensions are suppressed by the spin with the exponent and have the following coefficients 19 where we have defined We keep a = 1 and a = 2 to indicate explicitly that we deal with two OPE coefficients (not to be confused with orders in our large spin perturbation theory). The subscript b in the anomalous dimensions stays to indicate that it reproduces bosonic minimal-twist operators. The products of OPE coefficients have the following B exponents The exponents B XX and B Y Y are sub-leading to (3.29) and thus the products of OPE coefficients governed by this exponents are zero to the order we work in. In the leading order In the equations below one should make a replacement P a T → c −1 T P a T for the energy-momentum tensor, where cT > 0 is a real number which appears in the T T normalization. This happens because of the following. The normalization of 2-point functions is fixed only up to a positive real constant. In general we use this freedom to set the constant to 1. In the case of conserved operators we do not have such a freedom since the normalization of 2-point functions is completely fixed by the algebra they obey. To correct for this the above replacement is required.
where we have defined Finally the sub-leading order OPE coefficients are suppressed by large spin with the exponents (3.32) and the coefficients for n = 0 read as 20 n = 0 : We would like to note, that since our blocks are normalized with the factor (−1) according to (2.9), it was crucial to compensate this factor by another (−1) coming from the analytic continuation (2.24) to reproduce the identity operator in our bootstrap equations. Alternatively we could have worked with the s-t equations of another correlator as discussed around (2.26). In this case no analytic continuation is needed, however the products of OPE coefficients P will neatly organize themselves into positive (negative) definite quantities like in (3.14) times a factor (−1) . Which means that the final equations will never have a (−1) factor. This will allow to reproduce the identity operator correctly.

Equations: s-t-channel
The s-t bootstrap equation is formed by equating the expression (3.5) with (3.6). Schematically it has the form (2.18) and can be studied in the light-cone limit (2.27).
Applying the procedure of section 2.2, we have found that there are again two towers of large spin fermionic operators Ψ : X ( +1, ) , Y ( +1, ) (3.35) together with their hermitian conjugates X and Y . The operators (3.35) are in principle different from (3.22), we put a prime to indicate that. As before we determine their twists together with products of OPE coefficients In the leading order the twists read as where n is a non-negative integer spanning the tower of operators. 20 We have also found expression for n = 0, however their form is not presentable in the paper. We refer the reader to the accompanying code for details.
The coefficients P 0 X X and P 0 Y Y cannot contain a (−1) factor in order to reproduce correctly the absence of identity operator. We cannot however put further constraints on them, they are left undetermined. As a consequence the exponents B X X and B Y Y are also left undetermined.
The anomalous dimensions are suppressed by large with exponents D. Since the exponents B are unknown we can only compute the differences The coefficients of anomalous dimensions are given in terms of unknown P 0 as 21 where we have defined The sub-leading corrections to products of OPE coefficients are suppressed by the large spin with exponents C. Again, we can only determine the differences The coefficients P 1 are then given by n = 0 : where we have defined (3.44)

Summary
In previous sections we found that to satisfy the bootstrap equations in the u-t channel and in the s-t channel we are required to have operators (X, Y ) and (X , Y ) respectively. They look very similar: their leading twists Under the assumption that the identity operator in the u-t channel is fully reproduced only by the operators (X, Y ) we can argue as follows. Since the operators (X , Y ) cannot contribute to the u-channel one is forced to set P 0 X X = P 0 Y Y = 0, because these coefficients do not have a (−1) factor as was argued below (3.41) and would give a non-zero contribution otherwise. However this requirement is inconsistent with the coefficients of anomalous dimensions (3.39), (3.40) which would immediately blow up. Thus we conclude that (X, Y ) = (X , Y ).
We can then neatly combine the results of sections 3.2.1 and 3.2.2 as follows. The twists of large spin fermionic operators read as where the coefficients of anomalous dimensions γ b are given in (3.26), (3.27) and γ f are given in (3.39), (3.40). The OPE coefficients have the following form where P 0 XX and P 0 Y Y are given in (3.30), P 1 b XX and P 1 b Y Y are given in (3.33), P 1 f XX and P 1 f Y Y are given in (3.43) and P XX = P Y Y = 0 (3.49) to the same order of expansion. Looking at the final result (3.45)-(3.48) one notices that all the "fermionic" quantities (labeled by a superscript f ) contain a factor (−1) and thus their contribution to the ut bootstrap equation completely vanishes, while all the "bosonic" quantities (labeled by a superscript b) do not have a (−1) factor and thus their contribution completely vanishes in the s-t equation.

Generalized Free Theory
Here we consider a GFT with two primary fundamental fields: a real scalar φ and a fermion ψ. See footnote 3 for the definition of GFTs. All the double-twist composite primary operators Ψ ( ,¯ ) 22 are constructed from φ, ψ and their derivatives. We will not need their precise expressions and instead are interested only in their schematic form. For that consider a 3-point function φ(p 1 )ψ(p 2 )Ψ ( ,¯ ) (p 3 ). (3.50) It can be computed using Wick contractions and give a non vanishing result only if the operators Ψ ( ,¯ ) are of the following two types We choose the curly font for the operators (X , Y) in order to distinguish them from the operators (3.22) and (3.35). The twists of these operators are easy to read off from their schematic expressions In GFT one can also compute the 4-point function (3.2) using Wick contractions s-channel suites the best for determining the OPE coefficients of the operators (3.51) and (3.52). Equating the s-channel expansion (3.5) with (3.54) and splitting the obtained expression into independent tensor structures we get two equations of the schematic form where δ is the Kronecker symbol and f I ,n (z,z) represents conformal blocks together with the OPE coefficients. Expanding these equations around z = 0 andz = 0 and equating coefficients in front of the same powers of z andz on both sides, one gets a system of linear equations for the coefficients P X X and P YY . Only a finite number of coefficients P contribute at each order of expansion, thus, the system is exactly solvable. We found a solution for several values of and n and guessed the generic solution, which is , together with P GF T X X (n, ) = P GF T YY (n, ) = 0. (3.57) Notice that signs of the above expression neatly match the signs in (3.14) providing a nice consistency check.
t-channel is very similar to the s-channel. We can equate the t-channel expansion (3.6) with (3.54), split the resulting expression in two independent tensor structures and get two equations of the schematic form¨( for the OPE coefficients. These equations are apparently very different from (3.55) 23 but are in fact equivalent. Expanding (3.55) around z =z = 1 one finds a solution for the OPE coefficients which is identical to (3.56). This is another nice consistency check.
Connection to the Analytic Bootstrap Taking the leading behavior of (3.47), (3.48) and comparing it with (3.56) one can notice that the following holds This happens because the equations for the GFT OPE coefficients match the bootstrap equations in the light-cone limit. More precisely, take the t-channel GFT equation (3.58) and change variables according to (2.21) and the prescription below it. This expression will then be equivalent to (2.25) in the light-cone limit (3.21). Thus both equations should have the same asymptotic solution.

Conclusions
In this paper we discuss the analytic bootstrap approach [1,2] to general 4-point functions in 4D CFTs [3,4,31,32]. To simplify computations in practice we develop and attach to the paper a little Mathematica code "analyticBootstrap4D.m". We proceed by studying in detail the s-t and u-t bootstrap equations of the scalar-fermion correlator φ ψ φ ψ and show how the results from these two channels complement each other. We find that for every pair of scalar φ and fermion ψ operators in the spectrum correspond two infinite towers of large spin fermionic operators. We compute their twist and products of OPE coefficients in the leading and sub-leading order (3.45)-(3.48). The sub-leading corrections are governed by the minimal-twist operators in the CFT. They can be both bosonic (light scalars, conserved currents and the energy-momentum tensor) and fermionic (light (1, 0) and (0, 1) fermions and the supersymmetric current). We provide details of all the computations in the attached Mathematica file "Example.nb".
To complement the discussion we have also studied the scalar-fermion GFT. We derive closed-form expression for the GFT data (3.53) and (3.56). We show that in the large spin limit the GFT data perfectly matches the leading twist and OPE coefficients obtained using the analytic bootstrap.
In the same way of this paper one can study other spinning correlators. To complement numerical bootstrap studies in future it might be interesting to consider correlators of the form ψ ψψψ and JJJJ , where J is the conserved current. One can also study φJJ φ , where J is (2,1) fermion (spin 3/2 fermion). In case J is conserved, analogously to [87,88] one can access the supersymmetric conformal collider bounds.
The results of section 3 can be interpreted in the context of AdS. One can say that the double-twist operators X and Y on the CFT side are dual to a system of largely separated scalar and fermion particle orbiting each other with the large relative angular momentum . 24 The anomalous dimensions of operators X and Y would then correspond in AdS to energy shifts of the original system caused by interactions between orbiting particles. For our system we have two types of such interactions bosonic attractive/repulsive forces (gravity, Higgs-like force) or fermionic interactions which cause swapping of orbiting particles. It would be interesting to compute the energy shifts due to these interactions in AdS [82,86] 25 and compare them to anomalous dimensions found in this work. Here we only notice that the anomalous dimensions induced by the energy-momentum tensor are negative. This can be easily seen by plugging (3.15) into (3.26) and (3.27), indicating that gravity in AdS between scalar and fermion particle is attractive. Recent discovery of the Lorentizan inversion formula [74] for scalars brings analytic studies to a new level. Its generalization to spin is a matter of time. Our results might serve as a benchmark for checking future spinning inversion formula in 4D. 24 This system will also have two states ⊗ 1 2 = ( − 1 2 ) ⊕ ( + 1 2 ). Here we use half-integer values to label representations. It is different from the rest of the paper where we use only integer values. 25 In doing that techniques similar to [101] might be of help.
Combined with an index-free spinor notation  is the energy-momentum tensor, the charge D gives one constraint on the OPE coefficients 2λ 1 Oψψ − λ 2 Oψψ = 4i∆ ψ 3π 2 , (A.13) 27 We have defined the following short-hand notation for derivatives in auxiliary spinors The extra minus signs in the above expressions are needed to define the standard raising and lowering procedure in a consistent way. We also define coordinates in spinor notations as x αβ ≡ xµσ µ αβ andxα β ≡ x µσαβ µ . 28 Analogous expression can also be written for M µν and K µ .
where we have used [D, ψ(0, t)] = ∆ ψ ψ(0, t). One can verify that the result (3.15) obeys this constraint. To completely fix both OPE coefficients the charge D is obviously not enough and one should use other conserved charges, for example P µ . Unfortunately its action on the leading OPE term (A.12) vanishes and one should extend the OPE to sub-leading order. We find however that working with full correlation function at this point is more straightforward. We do not report the details, but only the final result (3.15).