Luis F Alday

Gregory P Korchemsky

Revisiting instanton corrections to the Konishi multiplet

We revisit the calculation of instanton effects in correlation functions in N = 4 SYM involving the Konishi operator and operators of twist two. Previous studies revealed that the scaling dimensions and the OPE coefficients of these operators do not receive instanton corrections in the semiclassical approximation. We go beyond this approximation and demonstrate that, while operators belonging to the same N = 4 supermultiplet ought to have the same conformal data, the evaluation of quantum instanton corrections for one operator can be mapped into a semiclassical computation for another operator in the same supermultiplet. This observation allows us to compute explicitly the leading instanton correction to the scaling dimension of operators in the Konishi supermultiplet as well as to their structure constants in the OPE of two half-BPS scalar operators. We then use these results, together with crossing symmetry, to determine instanton corrections to scaling dimensions of twist-four operators with large spin.

Introduction

Recently impressive progress has been achieved in understanding the properties of fourdimensional maximally supersymmetric N = 4 Yang-Mills theory in the planar limit, see [START_REF] Beisert | Review of AdS/CFT Integrability: An Overview[END_REF]. Thanks to integrability of the theory in this limit, it becomes possible to compute various quantities for an arbitrary 't Hooft coupling constant. At weak coupling, the resulting expressions agree with the results of explicit perturbative calculation whereas at strong coupling they match the predictions coming from the AdS/CFT correspondence. Much less is known however about the properties of N = 4 SYM beyond the planar limit and in particular about nonperturbative effects induced by instanton corrections.

The motivation for studying instanton corrections is multifold. Firstly, N = 4 SYM possesses S-duality [START_REF] Montonen | Magnetic Monopoles as Gauge Particles?[END_REF][START_REF] Sen | Dyon -monopole bound states, selfdual harmonic forms on the multimonopole moduli space, and SL(2,Z) invariance in string theory[END_REF][START_REF] Witten | Some comments on string dynamics[END_REF], namely, invariance under modular SL(2, Z) transformations acting on the complexified coupling constant

τ = θ 2π + 4πi g 2 . (1.1)
Instantons are expected to play a crucial role in restoring the invariance of the spectrum of scaling dimensions under the S-duality group. Secondly, previous studies revealed a remarkable similarity between instanton corrections to correlation functions in N = 4 SYM at weak coupling and the dual supergravity amplitudes induced by D-instantons in type IIB string theory [START_REF] Bianchi | Instantons in supersymmetric Yang-Mills and D instantons in IIB superstring theory[END_REF][START_REF] Dorey | Yang-Mills instantons in the large N limit and the AdS / CFT correspondence[END_REF][START_REF] Dorey | Multi-instanton calculus and the AdS / CFT correspondence in N=4 superconformal field theory[END_REF]. This suggests that the AdS/CFT correspondence can be tested beyond the planar limit to include instanton effects. Finally, the crossing symmetry of correlation functions leads to nontrivial constraints for the conformal data of the theory. They have been used in [START_REF] Beem | The N = 4 Superconformal Bootstrap[END_REF] to derive bounds for the scaling dimensions of leading twist operators of various spins. These bounds are expected to be saturated at fixed points of the S-duality group [START_REF] Beem | Resummation and S-duality in N=4 SYM[END_REF][START_REF] Alday | Modular interpolating functions for N=4 SYM[END_REF], were instanton contributions cannot be neglected.

In this paper, we revisit the calculation of instanton corrections to various correlation functions in N = 4 SYM at weak coupling. To compute such corrections we follow the standard approach (see reviews [START_REF] Dorey | The Calculus of many instantons[END_REF][START_REF] Belitsky | Yang-Mills and D instantons[END_REF][START_REF] Bianchi | Instantons and Supersymmetry[END_REF]). Namely, we decompose all fields into the sum of classical instanton solutions and fluctuations and then integrate out the latter. In the semiclassical approximation the quantum fluctuations can be neglected and the correlation functions can be reduced to finite-dimensional integrals over the collective coordinates of instantons

O(1) . . . O(n) inst = dµ phys e -S inst O(1) . . . O(n) , (1.2)
where all fields on the right-hand side are replaced by their expressions on the instanton background.

In the simplest case of the SU (2) gauge group, the one-instanton solution depends on bosonic collective coordinates ρ and x 0 defining the size of the instanton and its location as well as on 16 fermionic coordinates ξ A α and ηA α (with A = 1, . . . , 4 and α, α = 1, 2), reflecting the invariance of the equations of motion under N = 4 superconformal transformations. The corresponding integration measure over the collective coordinates for the one-instanton sector in the SU (2) gauge group is [START_REF] Bianchi | Instantons in supersymmetric Yang-Mills and D instantons in IIB superstring theory[END_REF] dµ phys e -S inst = g 8 2 34 π 10 e 2πiτ d 4 x 0 dρ ρ 5 d 8 ξ d 8 η .

(1.3)

For the correlation function (1.2) to be different from zero, the product of operators O(1)...O(n) should absorb all 16 fermion modes. This property can be used to show the vanishing of leading instanton corrections to various correlation functions. The relations (1.2) and (1.3) can be generalized to the SU (N ) gauge group for the one-instanton solution [START_REF] Dorey | Yang-Mills instantons in the large N limit and the AdS / CFT correspondence[END_REF][START_REF] Dorey | The Calculus of many instantons[END_REF] and for multiinstanton solutions at large N , see [START_REF] Dorey | Multi-instanton calculus and the AdS / CFT correspondence in N=4 superconformal field theory[END_REF].

Applying this approach we can systematically take into account the instanton effects to correlation functions and, then, extract the corresponding corrections to the conformal data of the theory. In particular, the instanton corrections to the scaling dimensions of operators have the following general form [START_REF] Bianchi | Instantons and Supersymmetry[END_REF]:

1 γ inst = n≥1 e 2πinτ + e -2πinτ ∞ k=0 γ n,k (N )g 2k , (1.4)
where the two terms inside the brackets, e 2πinτ and e -2πinτ , describe the leading contribution of n instantons and n anti-instantons, respectively, and perturbative fluctuations produce subleading corrections suppressed by powers of the coupling constant g 2 . Previous studies of four-point and two-point correlation functions revealed [START_REF] Bianchi | On the logarithmic behavior in N=4 SYM theory[END_REF][START_REF] Arutyunov | Perturbative and instanton corrections to the OPE of CPOs in N=4 SYM(4)[END_REF][START_REF] Bianchi | Properties of the Konishi multiplet in N=4 SYM theory[END_REF][START_REF] Kovacs | On instanton contributions to anomalous dimensions in N=4 supersymmetric Yang-Mills theory[END_REF] that for many operators in N = 4 SYM, including those with lower bare dimension (Konishi operator) and twist-two operators, the leading instanton corrections vanish, γ n,0 (N ) = 0. Going beyond the semiclassical approximation, one can envisage two possible scenarios: (i) the instanton corrections vanish to all orders in the coupling constant, γ n,k = 0, due to some symmetry, or (ii) the instanton corrections do not vanish, namely γ n,k = 0 for k ≥ K, but they are suppressed by a power of the coupling constant g 2K . The first scenario seems to be incompatible with the expected S-duality of N = 4 SYM, whereas to test the second scenario would require to take into account quantum corrections making the calculation of instanton effects extremely complicated. This explains, in part, why little progress has been made in improving the existing results over the last decade.

In this paper we demonstrate, for the first time, that the scaling dimensions of the Konishi operator as well as other members of the corresponding N = 4 supermultiplet receive instanton corrections at order O(g 4 ), that is γ n,0 = γ n,1 = 0 but γ n,2 = 0. We identify the leading nonvanishing correction and compute the corresponding coefficient γ n,2 . We also evaluate the three-point correlation function of the Konishi operator and two half-BPS scalar operators and show that it receives a nonvanishing instanton correction at order O(g 2 ). While operators of the same supermultiplet ought to have the same anomalous dimension (and related OPE coefficients with two half-BPS operators), we observe that quantum instanton computations for some operators map to semiclassical instanton computations for others! This allows us to make progress.

Using these results, we obtain the instanton contribution to the asymptotic behavior of the four-point correlation function of half-BPS operators in the light-cone limit and, then, employ crossing symmetry to compute the instanton corrections to twist-four operators with high spin.

The paper is organized as follows. In Sect. 2 we review the conventional instanton calculus in N = 4 SYM and apply it to compute instanton effects to various correlation functions involving the Konishi operator. We show that, due to a different coupling constant dependence of the leading instanton corrections to two-and four-point correlation functions, it is to possible to compute O(g 4 ) corrections to the scaling dimension of the Konishi operator. Furthermore, we consider anomalous dimensions as well as OPE coefficients (with two half-BPS operators) corresponding to general twist-two operators. We show that at order O(g 2 ) only the OPE coefficient corresponding to the Konishi supermultiplet gets instanton constributions. In Sect. 3 we use this information in order to infer the instanton contribution to the light-cone asymptotics of correlation function of four half-BPS operators. This information, together with crossing symmetry, is then used to compute the instanton corrections to twist-four operators with large spin. Concluding remarks are presented in Sect. 4. Useful definitions are included in two appendices.

Instanton corrections to correlation functions

In this Section, we evaluate instanton corrections to two-and three-point correlation functions of various operators in N = 4 SYM theory in the semiclassical approximation. As was explained in the previous section, the calculation amounts to evaluating the product of operators in the background of instantons and integrating the resulting expression over the collective coordinates.

Instanton in N = 4 SYM

The Lagrangian of N = 4 super Yang-Mills theory describes a gauge field A µ , (anti)gaugino fieds λ A α and λ α A , as well as scalars φ AB satisfying the reality condition φAB =

1 2 ǫ ABCD φ CD L = 1 g 2 tr - 1 16 F 2 αβ - 1 16 F 2 α β - 1 4 D α α φ AB D α α φAB -2i λ αA D αβ λ A β + √ 2λ αA [ φAB , λ B α ] - √ 2 λ αA [φ AB , λ α B ] + 1 8 [φ AB , φ CD ][ φAB , φCD ] + i θ 8π 2 tr 1 16 F 2 αβ - 1 16 F 2 α β . (2.1)
Here we used spinor notations (see Appendix A for our conventionts) and denoted by F αβ and F α β the (anti) self-dual part of gauge field strength tensor

F µν = -i[D µ , D ν ] with covariant derivative D µ = ∂ µ + i[A µ , ] and D α α = D µ (σ µ ) α α.
All fields take value in the SU (N ) algebra, e.g. A µ = A a µ T a , with the generators satisfying [T a , T b ] = if abc T c and normalized as tr(T a T b ) = 1 2 δ ab . By definition, the instanton is a solution to the classical equations of motion. Due to our choice of normalisations in the Lagrangian (2.1), all elementary fields in the instanton background, A µ , λ A α , λ α A and φ AB , are independent of the coupling constant. Their explicit expressions for the SU (N ) gauge group are rather complicated and only few terms in their expansion in powers of 8N fermionic collective modes are currently known [START_REF] Dorey | The Calculus of many instantons[END_REF][START_REF] Belitsky | Yang-Mills and D instantons[END_REF][START_REF] Bianchi | Instantons and Supersymmetry[END_REF]. Significant simplification occurs however for the SU (2) gauge group. In this case, the general oneinstanton solution can be obtained by applying N = 4 superconformal transformations to a special solution to the equations motion

A (0) µ = 2 η a µν (x -x 0 ) ν T a (x -x 0 ) 2 + ρ 2 , φ AB,(0) = λ A,(0) α = λ α,(0) A = 0 , (2.2)
where

A (0)
µ is the well-known one-instanton solution in pure Yang-Mills theory. It depends on the collective coordinates ρ and x 0 defining the size and the position of the instanton, respectively. Here η a µν are the 't Hooft symbols and the SU (2) generators are related to Pauli matrices T a = σ a /2.

The field configuration (2.2) is annihilated by half of the superconformal generators Q αA and S A α . Applying the remaining superconformal transformations exp

ξ A α Q α A + η αA S αA to (2.
2), we obtain a solution to the full classical equations of motion that depends on 16 fermionic collective coordinates, ξ A α and η αA . By virtue of the SU (4) R-symmetry, the resulting expression for the instanton configuration takes the following form

A µ = A (0) µ + A (4) µ + • • • + A (16)
µ , λ αA = λ αA,(1) + λ αA, (5) (13) , (14) ,

+ • • • + λ αA,
φ AB = φ AB,(2) + φ AB,(6) + • • • + φ AB,
λ αA = λ(3) αA + λ(7) αA + • • • + λ(15) αA , (2.3)
where

A (n)
µ denotes the component of the gauge field that is homogenous in ξ A α and η αA of degree n and similar for other fields.

The leading terms of the expansions (2.3) have been worked out in Ref. [START_REF] Belitsky | Yang-Mills and D instantons[END_REF]. For our purposes we will also need subleading terms. Their direct calculation is more involved, e.g. finding φ AB, (6) amounts to applying Q-and S-transformations to (2.2) subsequently six times. There is, however, a shortcut that simplifies the calculation significantly. Namely, the subleading corrections depend on the instanton field A (0)

α α = iA (0)
µ (σ µ ) α α and fermion collective coordinates, ξ A α and η αA . It turns out that the requirement for the fields (2.3) to have the correct properties with respect to conformal symmetry, R-symmetry and gauge transformations, fixes their general form up to a few constants. The latter can be determined by requiring the fields (2.3) to satisfy the classical equations of motion derived from (2.1). Going through the calculation we have found the expressions for the subleading corrections to gauge field A (4)

α α = iA (4)
µ (σ µ ) α α and scalar φ (6),AB . Together with the leading correction φ (2),AB they are given by A (4)

α α = - 1 12 ǫ ABCD ζ A α ζ βB (ζD β αF ζ) CD - 1 2 ǫ ABCD ζ A α ηB α (ζF ζ) CD , φ (2),AB = 1 √ 2 (ζF ζ) AB , φ (6),AB = - 1 20 √ 2 ǫ CDEF (ζ 2 ) AC (ζ 2 ) BD (ζF 2 ζ) EF , (2.4)
where we have introduced a short-hand notation for a particular linear x-dependent combination of fermionic modes

ζ A α (x) = ξ A α + x α α η αA , (2.5)
and for various Lorentz contractions of fermion modes

(ζ 2 ) AB = (ζ 2 ) BA = ζ βA ǫ βγ ζ γB , (ζF ζ) AB = -(ζF ζ) BA = ζ αA F αβ ζ βB , (ζF 2 ζ) AB = -(ζF 2 ζ) BA = ζ αA F αβ ǫ βδ F δγ ζ γB . (2.6)
Here F αβ = F βα is the strength tensor for the SU (2) instanton2

(F αβ ) i j = 8ρ 2 [(x -x 0 ) 2 + ρ 2 ] 2 (ǫ iα δ j β + ǫ iβ δ j α ) , (2.7)
where i and j are the SU (2) indices.

A peculiar feature of φ (6),AB in (2.4) is that it depends on the fermionic modes only through the variable ζ defined in (2.5). This is not the case however for the gauge field A (4) α α. The difference is due to different transformation properties of fields under conformal transformations. By virtue of superconformal invariance, the action of N = 4 SYM evaluated on the instanton configuration (2.3) does not depend on the fermionic modes ξ and η and is given by

S inst = d 4 x L(x) = 8π 2 g 2 -iθ = -2πiτ , (2.8)
where τ is the complex coupling constant (1.1).

Normalization of operators

Later in this section we shall compute the leading instanton corrections to scaling dimensions and OPE coefficients of various composite gauge invariant operators built from scalar fields φ AB and φAB = 1 2 ǫ ABCD φ CD . Before doing this, we have to carefully examine the normalization of the operators. The reason for this is that, due to our definition of the Lagrangian (2.1), the free scalar propagator depends on the coupling constant

φa AB (1)φ b,CD (2) = g 2 δ ab δ C A δ D B -δ D A δ C B D(x 12 ) , φ a,AB (1)φ b,CD (2) = g 2 δ ab ǫ ABCD D(x 12 ) , (2.9)
where we used a shorthand notation for

x 12 = x 1 -x 2 and D(x) = 1/(4π 2 x 2 )
. Taking this into account, we define the simplest scalar operators of bare dimension 2

O 20 ′ (x, Y ) = 1 g 2 Y AB Y CD tr(φ AB φ CD ) , K(x) = 1 g 2 tr( φAB φ AB ) , (2.10)
where Y AB is an antisymmetric SU (4) tensor satisfying ǫ ABCD Y AB Y CD = 0. It was introduced to project the product of two scalar fields onto the irreducible SU (4) representation 20 ′ . The half-BPS operator O 20 ′ is annihilated by half of the N = 4 supersymmetries and its scaling dimension is protected from quantum corrections. The Konishi operator K is the simplest unprotected operator.

It is also convenient to introduce the following operator of bare dimension 4

K ′ (x) = 1 g 4 tr([Z, X][Z, X]) , (2.11)
where we used the standard notation for complex scalar fields Z = φ 14 and X = φ 24 . The operator K ′ is a supersymmetric descendant of the Konishi operator,

K ′ ∼ δ 2 Qδ 2
Q K, and, as a consequence, its anomalous dimension and OPE coefficient with two half-BPS operators coincide with those of the Konishi operator.

The definition of the operators in (2.10) and (2.11) involves inverse powers of the coupling constant, one per each scalar field. They were introduced in order to ensure that the correlation functions scale as O(g 0 ) in the Born approximation. Indeed

K(1)K(2) Born = 12(N 2 -1)D 2 (x 12 ) , K ′ (1) K′ (2) Born = 3 4 N 2 (N 2 -1)D 4 (x 12 ) , O 20 ′ (1)O 20 ′ (2) Born = 1 2 (y 2 12 ) 2 (N 2 -1)D 2 (x 12 ) , (2.12)
where y 2 12 ≡ ǫ ABCD Y 1,AB Y 2,CD . The presence of different powers of the coupling constant in the definition of the operators (2.10) and (2.11) leads to important consequences for instanton corrections to the correlation functions. Since the fields in the instanton background (2.3) do not depend on the coupling constant, the dependence on g 2 of the correlation function in the semiclassical approximation comes solely from the SU (2) integration measure over the moduli of instantons (1.3) and the from powers of 1/g 2 accompanying each operator. In this way, we find

K(1)K(2) inst = O(g 4 e 2πiτ ) , K ′ (1) K′ (2) inst = O(e 2πiτ ) , (2.13)
where the extra factor of g 4 in the first relation arises due to a different power of 1/g 2 in the definitions (2.10) and (2.11).

In general, two-point correlation functions develop logarithmic singularities and generate corrections to anomalous dimensions of operators. Then, assuming that the correlation functions (2.13) are different from zero, we would deduce that instanton corrections to anomalous dimensions of operators K and K ′ should have different dependence on the coupling constant:

γ K ′ = O(e 2πiτ ) whereas γ K = O(g 4 e 2πiτ
). However, this contradicts the fact that the two operators belong to the same supermultiplet and, therefore, their anomalous dimensions have to coincide. In other words, the superconformal symmetry dictates that γ K ′ has to vanish in the semiclassical approximation. 3 To get a nonzero result for γ K ′ we have to go beyond this approximation and take into account quantum fluctuation around instanton onfigurations. The calculation of quantum corrections to K ′ (1) K′ (2) inst is way more complicated but the resulting expression for γ K ′ ought to match γ K = O(g 4 e 2πiτ ) obtained in the semiclassical approximation. We shall use this observation to compute the leading instanton correction to the scaling dimension of the Konishi operator in Sect. 2.4.

To observe another interesting feature of the Konishi operator, we examine the coupling dependence of the following correlation functions

O 20 ′ (1)O 20 ′ (2)K(3) inst = O(g 2 e 2πiτ ) , O 20 ′ (1)O 20 ′ (2)O 20 ′ (3)O 20 ′ (4) inst = O(e 2πiτ ) , (2.14)
where each operator brings in the factor of 1/g 2 multiplied by the factor of g

Instanton profile of operators

As the first step, we evaluate the operators (2.10) and (2.11) in the instanton background for the SU (2) gauge group. We start with the Konishi operator (2.10) and replace scalar fields by their expressions (2.3). This leads to

K = 1 g 2 tr(φ AB φAB ) = 1 g 2 K (4) + K (8) + K (12) + K (16) , (2.15)
where K (n) denotes the contribution containing n fermion modes. Notice that K (n) is independent of the coupling constant. Since the scalar field has at least two fermion modes, the expansion starts with K (4) . By virtue of the SU (4) symmetry, the number of fermion modes in the subsequent terms of the expansion increases by four units. The last term of the expansion contains the maximal number of fermion modes. We find, in agreement with [START_REF] Bianchi | Properties of the Konishi multiplet in N=4 SYM theory[END_REF], that the first term on the right-hand side of (2.15) vanishes

K (4) (x) = tr [φ (2),AB φ(2) AB ] ∼ ǫ ABCD (ζ 2 ) AB (ζ 2 ) CD = 0 , (2.16)
where ζ = ξ + xη is a linear combination of fermion modes defined in (2.5). Here in the second relation we took into account that (ζ 2 ) AB = ζ αA ζ B α is symmetric with respect to SU (4) indices. As a result, the expression for the Konishi operator contains eight fermion modes at least and the leading term is given by

K (8) (x) = 2 tr [φ (2),AB φ(6) AB ] = -3 2 × 2 15 × ρ 6 [ρ 2 + (x -x 0 ) 2 ] 6 [ζ(x)] 8 , (2.17)
where ζ 8 = A,α ζ αA is the product of eight fermion modes. Expressions for higher components of (2.15) are more complicated but we do not need them for our purposes.

Let us consider the half-BPS operator O 20 ′ (x, Y ) defined in (2.10). Since this operator is annihilated by half of the N = 4 supercharges, it depends on four ξ and four η fermion modes. As a consequence, its expansion in powers of Grassmann variables is shorter as compared with (2.15)

O 20 ′ (x, Y ) = 1 g 2 Y AB Y CD tr(φ AB φ CD ) = 1 g 2 O (4) 20 ′ + O (8) 20 ′ , (2.18)
where the two terms on the right-hand side involve four and eight fermion modes respectively, and are given by

O ( 4
)
20 ′ = Y AB Y CD tr(φ (2),AB φ (2),CD ) , O (8)
20 ′ = 2Y AB Y CD tr(φ (2),AB φ (6),CD ) . (2.19)
As was already mentioned, the scalar fields φ (2),AB and φ (6),AB depend on fermion modes only through their linear combination (2.5) and the same is obviously true for the components (2.19). This property alone implies that O (8)

20 ′ has to vanish. Indeed, if O (8)
20 ′ were different from zero, it would be proportional to ζ 8 = (ξ + xη) 8 and, therefore, would contain terms with more than four ξ's and four η's, in contradiction with half-BPS condition. The relation O [START_REF] Beem | The N = 4 Superconformal Bootstrap[END_REF] 20 ′ = 0 can be also verified by a direct calculation. 4 Thus, the operator O 20 ′ (x, Y ) contains exactly 4 fermion ζ-modes and is given by

O 20 ′ (x, Y ) = 1 g 2 128ρ 4 [ρ 2 + (x -x 0 ) 2 ] 4 × Y AB Y CD (ζ 2 ) AC (ζ 2 ) BD .
(2.20)

This feature allows us to verify the known properties of correlation functions of half-BPS operators.

We recall that in order for a correlation function to be different from zero, the product of operators should involve terms containing sixteen fermion modes. Since the product of two and three half-BPS operators have 8 and 12 fermion modes, respectively, the corresponding correlation functions

O 20 ′ (1)O 20 ′ (2) and O 20 ′ (1)O 20 ′ (2)O 20 ′ (3)
do not receive instanton corrections in the semiclassical approximation. This result is in agreement with the known fact that the above mentioned correlation functions are protected from quantum corrections and are given by their Born level expressions. The simplest correlation function that receives instanton correction involves four half-BPS operators. We shall return to this correlation function in Sect. [START_REF] Sen | Dyon -monopole bound states, selfdual harmonic forms on the multimonopole moduli space, and SL(2,Z) invariance in string theory[END_REF].

Let us finally examine the operator K ′ defined in (2.10). Replacing the scalar fields by their expressions (2.3), we find in a similar manner

K ′ = 1 g 4 tr ([Z, X][Z, X]) = 1 g 4 K ′ (8) + K ′ (12) + K ′ (16) , (2.21)
where the lowest term involves eight fermion modes, ) , X (2) ][Z (2) , X (2) ] , K ′ (12) = 2 tr [Z (6) , X (2) ][Z (2) , X (2) ] + 2 tr [Z (2) , X (6) ][Z (2) , X (2) ] .

K ′ (8) = tr [Z (2
As in the previous case, the dependence on fermion modes enters into these expressions through their linear combination ζ A α defined in (2.5). According to (2.4), the scalar field φ (2),AB is given by the product of two fermion modes φ (2),AB ∼ ζ A ζ B . Taking into account that Z = φ 14 and X = φ 24 , we obtain that K ′ (8) ∼ (ζ 4α ζ 4 α ) 2 = 0, which is consistent with the findings of [START_REF] Kovacs | On instanton contributions to anomalous dimensions in N=4 supersymmetric Yang-Mills theory[END_REF]. In a similar manner, we can show that K [START_REF] Arutyunov | Perturbative and instanton corrections to the OPE of CPOs in N=4 SYM(4)[END_REF] depends on all 16 fermion ξ-and η-modes. However their product ξ 8 η8 is the SU (4) singlet and can not contribute to K ′ (16) due to mismatch of the SU (4) quantum numbers leading to K ′ (16) = 0. We therefore conclude that K ′ = 0 in the instanton background and, as a consequence, all correlation functions involving this operator vanish in the semi-classical approximation.

′ (12) ∼ (ζ 4α ζ 4 α )ζ 4 β = 0. The top component K ′ (

Instanton corrections to Konishi operator

We are now ready to evaluate the leading instanton corrections to correlation functions involving the Konishi operator for the SU (2) gauge group. We start with the two-point function

K(1)K(2) inst = dµ phys e -S isnt K(1)K(2) = e 2πiτ g 4
dµ phys K (8) (1)K (8) (2) ,

where the SU (2) integration measure dµ phys is defined in (1.3). Here in the first relation we replaced operators by their instanton profile (2.15) and in the second relation retained terms involving 16 fermion modes. The terms with higher number of fermion modes do not contribute. Replacing K (8) with (2.17) we get

K(1)K(2) inst = 81 16π 10 g 4 e 2πiτ d 4 x 0 dρ ρ 5 ρ 12 (ρ 2 + x 2 10 ) 6 (ρ 2 + x 2 20 ) 6 × d 8 ξ d 8 η 4 A=1 ζ A 1 (x 1 )ζ A 2 (x 1 )ζ A 1 (x 2 )ζ A 2 (x 2 ) , (2.23)
where ζ A α (x) = ξ A α + x α α η αA . Integration over Grassmann variables in the second line yields (x 2 12 ) 4 leading to

K(1)K(2) inst = 81 16π 10 g 4 e 2πiτ d 4 x 0 dρ ρ 5 (x 2 12 ) 4 ρ 12 (ρ 2 + x 2 10 ) 6 (ρ 2 + x 2 20 ) 6 , (2.24)
where the integration is over the size and position of the instanton. We verify that this expression has the expected dependence (2.13) on the coupling constant.

The integral on the right-hand side of (2.24) develops a logarithmic divergence that comes from integration over instantons of small size, ρ → 0, located close to one of the operators, |x 10 |→ 0 or |x 20 |→ 0. It indicates that the instanton corrections modify the scaling dimension of the operator. It is convenient to regularize the integral by modifying the integration measure over x 0

d 4 x 0 → d 4-2ǫ x 0 . (2.25)
The resulting integral in (2.24) is well-defined for ǫ < 0 and develops a simple pole as ǫ → 0. Combining (2.24) with the Born term (2.12) (evaluated for the SU (2) gauge group) we obtain

K(1)K(2) = 36 (4π 2 ) 2 (x 2 12 ) 2 1 - 9 5ǫ (x 2 12 ) -ǫ g 2 4π 2 2 e 2πiτ . ( 2

.26)

Following a standard procedure, we apply the dilatation operator i (x i ∂ i ) to both sides of this relation, and find that the correction to the scaling dimension of K is given by the residue at the pole. Thus, we conclude that the leading instanton correction to the anomalous dimension of the Konishi operator in the SU (2) gauge group is given by

γ K = - 9 5
g 2 4π 2 2 e 2πiτ + e -2πiτ , (2.27)
where we have added a complex conjugated term to take into account the contribution from the anti-instanton. Notice that γ K has negative sign. The result (2.27) holds for the SU (2) gauge group. Its generalization to SU (N ) gauge group will be discussed in Sect. 2.6.

For the three-point function of a Konishi operator and two half-BPS operators we can proceed analogously. We find

O 20 ′ (1)O 20 ′ (2)K(3) inst = e 2πiτ g 2 dµ phys O 20 ′ (1)O 20 ′ (2)K (8) (3) = - 9 32π 10 g 2 e 2πiτ ×I B × I F , (2.28)
where in the second relation we replaced operators by their expressions on the instanton background, Eqs. (2.17) and (2.20), and introduced a short-hand notation for the integrals over bosonic and fermion collective coordinates

I B = d 4 x 0 dρ ρ 5 ρ 14 (ρ 2 + x 2 10 ) 4 (ρ 2 + x 2 20 ) 4 (ρ 2 + x 2 30 ) 6 , I F = d 8 ξ d 8 η Y 1,AB Y 1,CD (ζ 2 (x 1 )) AC (ζ 2 (x 1 )) BD × Y 2,EF Y 2,KL (ζ 2 (x 2 )) EK (ζ 2 (x 2 )) F L 4 A=1 ζ A 1 (x 3 )ζ A 2 (x 3 ) , (2.29)
with ζ 2 (x) given by (2.6) and (2.5). Both integrals are well-defined and their dependence on x-and Y -variables is uniquely fixed by conformal and R-symmetry, respectively. Going through the calculation we get

I B = π 2 90 1 (x 2 13 x 2 23 ) 3 x 2 12 , I F = 9(x 2 13 x 2 23 ) 2 (y 2 12 ) 2 , (2.30)
where y 2 12 is defined in (2.12). Plugging (2.30) into (2.28) we obtain the final result for the instanton contribution. Taking into account the anti-instanton contribution and combining this with the Born term we obtain

O 20 ′ (1)O 20 ′ (2)K(3) = 6 (y 2 12 ) 2 (4π 2 ) 3 x 2 12 x 2 23 x 2 31 1 - 3g 2 10π 2 e 2πiτ + e -2πiτ . (2.31)
We observe that, in agreement with (2.27), the three-point function does not receive logarithmically divergent corrections due to anomalous dimension of the Konishi operator. The latter corrections scale as O(g 4 e 2πiτ ) and are subleading in (2.31). As a consequence, the leading instanton correction to the OPE coefficient

C K = O 20 ′ O 20 ′ K /( O 20 ′ O 20 ′ KK 1/2 )
K by the same factor that enters the right-hand side of (2.31)5

C K /C (0) K = 1 - 3g 2 10π 2 e 2πiτ + e -2πiτ .
(2.32)

Instanton corrections to twist-two operators

The twist-two operators provide the leading contribution to four-point correlation functions in the light-cone limit x 2 12 → 0. In N = 4 SYM these operators belong to the same supermultiplet which allows us to restrict our consideration to a particular twist-2 operator

O S = 1 g 2 tr(ZD S + Z) + . . . , (2.33)
which is built from the complex scalar field Z = φ 14 and the light-cone component of the covariant derivative D + = (nD) (with n 2 µ = 0). The ellipses denote similar terms with covariant derivatives distributed between the two scalar fields, their form fixed by conformal symmetry. For S = 0 the operator (2.33) coincides with the half-BPS operator O 20 ′ (x, Y ) for a special choice of Y -variables. For S = 2 it is a superconformal descendant of the Konishi operator K. As in the case of Konishi operator, the additional factor of 1/g 2 was introduced in (2.33) in order for the two-point correlation function of O S to scale as O(g 0 ) in the Born approximation.

In order to study instanton corrections to O S it is convenient to switch to nonlocal light-ray operators

O(z) = 1 g 2 tr [Z(0)E(0, z)Z(nz)E(z, 0)] = S≥0 z S S! [O S (0) + . . . ] , (2.34)
where two scalar fields are separated along the light-ray defined by a null vector n and lightlike Wilson lines were introduced to restore gauge invariance, (2.34) in order to project it onto the contribution of the operator O S (0). As before, we have to examine the product of operators in the background of the instanton and identify the contributions involving 16 fermion modes. Since the product of two half-BPS operators contains 8 modes, the remaining 8 modes should be soaked up by O(z). Replacing Z = Z (2) + Z (6) + . . . and A = A (0) + A (4) + . . . in (2.34) we find the corresponding contribution is given by O (8) (z) = tr Z (2) (0)E (0) (0, z)Z (6) (nz)E (0) (z, 0) + Z (6)

E(z 1 , z 2 ) = P exp i z 2 z 1 dt n • A(nt) . ( 2
(0)E (0) (0, z)Z (2) (nz)E (0) (z, 0) + i z 0 dt Z (2) (0)E (0) (0, t) n • A (4) (nt)E (0) (t, z)Z (2) (nz)E (0) (z, 0) + i 0 z dt Z (2) (0)E (0) (0, z)Z (2) (nz)E (0) (z, t) n • A (4) (nt)E (0) (t, 0) , (2.36)
where the subscript indicates the number of fermion zero modes and E (0) (z 1 , z 2 ) depends on A (0) .

Taking into account (2.4) and recalling that Z = φ 14 , we can evaluate O (8) (z) explicitly. Going through a lengthy calculation (the details are presented in [START_REF] Alday | Instanton corrections to twist-two operators[END_REF]) we obtain

O 20 ′ (1)O 20 ′ (2)O(z) inst = e 2πiτ g 2 dµ phys O 20 ′ (1)O 20 ′ (2)O (8) (z) = - 27 32π 10 g 2 e 2πiτ y 2 12 y 2 13 y 2 23 z 2 (nx 2 )x 2 1 -(nx 1 )x 2 2 2 D 4433 (x 1 , x 2 , 0, nz) , (2.37)
where y 2 ij are defined in (2.12) with all Y 3,AB vanishing except Y 3,14 = -Y 3,41 = 1/2. Here the product of y-variables keeps track of the R-charges of the operators whereas the nontrivial dependence on x-variables is described by D-function defined in Appendix B.

To extract the correlation function O 20 ′ (1)O 20 ′ (2)O S (0) , we expand the expression in the second line of (2.37) in powers of z and decompose it over the conformal partial waves. In this way we find that (2.37) receives a nonvanishing contribution from only one partial wave with S = 2. In other words, the three-point correlation function in the semiclassical approximation is different from zero only for twist-two operators with spin S = 2:

O 20 ′ (1)O 20 ′ (2)O S (0) inst = -δ S,2 9g 2 10π 2 e 2πiτ + e -2πiτ y 2 12 y 2 13 y 2 23 (4π 2 ) 3 x 2 12 x 2 1 x 2 2 2(nx 1 ) x 2 1 - 2(nx 2 ) x 2 2 2 , (2.38)
where we have added the contribution from the anti-instanton. For S = 0, this correlation function is protected from quantum corrections. For higher spin S > 2, the instanton corrections to (2.38) scale as O(g 4 e 2πiτ ) at least. For S = 2 we verified that the relation (2.38) divided by the Born level result coincides with the analogous expression for the Konishi operator (2.31). This is not surprising given the fact that the two operators O S=2 and K belong to the same N = 4 supermultiplet, but rather serves as a nontrivial check of our calculation.

It is straightforward to extend the above considerations to the two-point correlation function of twist-two operators. Computing the leading instanton correction to the two-point correlation function of light-ray operators (2.34) and projecting them onto operators O S with a help of (2.34), we obtain (see [START_REF] Belitsky | [END_REF] for details on the projection procedure)

O S (x) ŌS ′ (0) inst = -δ SS ′ δ S,2 81 5ǫ (x 2 ) -ǫ g 2 4π 2 2 e 2πiτ + e -2πiτ [2(xn)] 2S (4π 2 ) 2 (x 2 ) 2+2S .
(2.39)

In other words, the instanton corrections vanish for all spins except for S = 2. In the latter case, they generate the same correction to the scaling dimension of the operator O S=2 as to the Konishi operator (2.27). We recall that the two operators belong to the same supermultiplet and their anomalous dimension ought to coincide.

Generalization to the SU(N) gauge group

Having determined the contribution of a single (anti)instanton to correlation functions in N = 4 SYM for the SU (2) gauge group, we can now generalize the above results to the SU (N ) gauge group and, in addition, take into account the contribution of an arbitrary number of (anti)instantons at large N . The instanton for the SU (N ) gauge group has 8N fermion modes. Among them there are 16 exact supesymmetric and superconformal zero modes, ξ and η, respectively. The remaining 8N -16 'nonexact' fermion modes do not correspond to any symmetry and the corresponding SU (N ) instanton action S inst develops a nontrivial dependence on these modes. This leads to significant simplification in computing the correlation functions. As in the previous case, for the instanton correction to be different from zero, all fermion modes should be saturated. Then, in the semiclasscial approximation, the exact modes are absorbed by the instanton profile of the operators whereas the nonexact modes are saturated by S inst . As a consequence, the contribution of the nonexact modes to the correlation functions factorizes into a univeral N -dependent factor [START_REF] Dorey | Yang-Mills instantons in the large N limit and the AdS / CFT correspondence[END_REF][START_REF] Dorey | The Calculus of many instantons[END_REF] 6

O(1) . . . O(n) SU(N), 1-inst = κ N O(1) . . . O(n) SU(2), 1-inst , (2.40)
where κ N takes into account both the embedding of the SU (2) instanton in SU (N ) and integration over the nonexact modes

κ N = 2 3-2N Γ(2N -1) Γ(N )Γ(N -1) = 2Γ(N -1/2) √ π Γ(N -1) N →∞ → 2 √ N √ π . ( 2

.41)

In application to the Konishi operator, we can use (2.40) together with (2.27) and (2.32) to get its anomalous dimension and OPE coefficient for the SU (N ) gauge group

γ K = - 27κ N 5(N 2 -1)
g 2 4π 2 2 e 2πiτ + e -2πiτ , C K = C (0) K 1 - 9κ N 10(N 2 -1) g 2 4π 2 e 2πiτ + e -2πiτ .
(2.42)

Here we inserted an additional factor of 3/(N 2 -1) to account for N -dependence of twoand three-point correlation functions in the Born approximation (see Eq. (2.12)). The relations (2.42) can be further generalized to include the contribution of multiinstantons. As was shown in [START_REF] Dorey | Multi-instanton calculus and the AdS / CFT correspondence in N=4 superconformal field theory[END_REF], the calculation simplifies dramatically in the large N limit due to the fact that the integration over the moduli space of n-instantons is dominated by the saddle-point. Repeating the analysis of [START_REF] Dorey | Multi-instanton calculus and the AdS / CFT correspondence in N=4 superconformal field theory[END_REF], we find that in this limit the profile of the Konishi operator in the n-instanton background is proportional to its one-instanton expression. This makes the evaluation of instanton corrections to the Konishi operator very similar to that performed in [START_REF] Dorey | Multi-instanton calculus and the AdS / CFT correspondence in N=4 superconformal field theory[END_REF] for the half-BPS operator. In this way, going through the calculation we find the generalisation of (2.42)

γ K n-inst = - 54 N 3/2 g 2 4π 2 2 n 2 (Q n + Qn ) , C K n-inst = C (0) K 1 - 9 5N 3/2 g 2 4π 2 n 3 (Q n + Qn ) , (2.43)
where we added the contribution of n anti-instantons and introduced

Q n = n -7/2 √ π e 2πinτ d|n 1 d 2 , (2.44)
where the sum runs over the positive divisors of n. We would like to emphasize that the relations (2.43) hold up to corrections suppressed by powers of 1/N and g 2 . The latter come from taking into account quantum fluctuations around the instanton configuration.

3 Instanton corrections to higher spin operators from crossing symmetry

In the previous section we have computed the instanton corrections to the scaling dimensions of the Konishi and twist-two operators, as well as to the OPE coefficients defining their contribution to the product of two half-BPS operators. Using these results we can determine the leading instanton contribution to the four-point correlation function

G 4 = O 20 ′ (1) . . . O 20 ′ (4) (3.1)
at short distances, x 1 → x 2 , and in the light-like limit, x 2 12 → 0. In this section we use this information, together with crossing symmetry, in order to compute instanton corrections to twist-four operators with large spin.

Properties of the correlation function

The four-point correlation function of half-BPS operators in N = 4 SYM with SU (N ) gauge group has the following structure [START_REF] Eden | Partial nonrenormalization of the stress tensor four point function in N=4 SYM and AdS / CFT[END_REF]

G 4 = 2(N 2 -1) (4π 2 ) 4 y 4 12 y 4 34 x 4 12 x 4 34 × (G short + G long ) , (3.2)
where G short and G long denote the contributions from (semi-)short multiplets and from long multiplets, respectively. The former contribution does not depend on the coupling constant, whereas the latter can be expressed in terms of a single function A long (u, v) of the conformal cross-ratios

G long = (z -α)(z -ᾱ)(z -α)(z -ᾱ) (αᾱ) 2 × A long (u, v) v 2 . ( 3.3)
The prefactor carries the R-charge dependence of the operators and we have introduced the following notation

u = x 2 12 x 2 34 x 2 13 x 2 24 = zz , v = x 2 23 x 2 14 x 2 13 x 2 24 = (1 -z)(1 -z) (3.4)
and similar for αᾱ = y 2 12 y 2 34 /(y 2 13 y 2 24 ) and (1α)(1ᾱ) = y 2 23 y 2 14 /(y 2 13 y 2 24 ). The function A long (u, v) admits a decomposition in terms of super-conformal blocks [START_REF] Dolan | Superconformal symmetry, correlation functions and the operator product expansion[END_REF]

A long (u, v) = v 2 S, ∆ a ∆,S u (∆-S)/2 g ∆+4,S (u, v) , (3.5)
where the sum runs over superconformal primary operators (and hence in the singlet of SU (4)) with even Lorentz spin S and scaling dimension ∆ ≥ S + 2 and a ∆,S is the square of the canonically normalised OPE coefficient. The contribution from super-conformal descendants is taken into account by the super-conformal blocks, where

g ∆,S (u, v) = - 1 2 S 1 z -z z S+1 k ∆+S (z)k ∆-S-2 (z) -zS+1 k ∆+S (z)k ∆-S-2 (z) , (3.6)
with k β (z) = 2 F 1 (β/2, β/2, β; z) and complex z and z variables defined in (3.4).

It is convenient to decompose A long (u, v) into the free-theory result A Born plus the quantum (coupling dependent) contribution

A A long (u, v) = A Born (u, v) + A(u, v) . (3.7)
The explicit expression for A Born (u, v) is not needed for our purposes, but it can be derived from the analysis of [START_REF] Dolan | Conformal partial wave expansions for N=4 chiral four point functions[END_REF]. At weak coupling, the expansion of A(u, v) runs in powers of 't Hooft coupling constant a = g 2 N/(4π 2 ) and (anti) instanton weight factors, e 2πiτ and e -2πiτ .

To leading order in these parameters we have

A(u, v) = - a 4 uv D1111 (u, v) + 15 κ N 2(N 2 -1) e 2πiτ + e -2πiτ u 2 v 2 D4444 (u, v) + . . . , (3.8)
where the D-functions are introduced in appendix B. Here the dots denote subleading terms suppressed by powers of the expansion parameters. Higher order perturbative corrections to A(u, v) were found in [START_REF] Eden | Hidden symmetry of four-point correlation functions and amplitudes in N=4 SYM[END_REF][START_REF] Eden | Constructing the correlation function of four stress-tensor multiplets and the four-particle amplitude in N=4 SYM[END_REF][START_REF] Drummond | Leading singularities and off-shell conformal integrals[END_REF].

Invariance of G 4 under the exchange of any pair of points leads to the crossing symmetry relations 7

A(u, v) = A(v, u) = v 2 A u v , 1 v . (3.9)
Each term on the right-hand side of (3.8) satisfies this relation. 7 While the free theory contributions G short and A Born (u, v) mix with each other.

-17 -

Instanton corrections to light-cone asymptotics

In the light-like limit x 2 12 → 0, or equivalently u → 0, the leading asymptotic behavior of (3.7) comes from the contribution of twist-two operators with scaling dimension ∆

S = 2 + S + γ S A long (u, v) = uv 2 S/2∈Z + a ∆ S ,S u γ S /2 f ∆ S +4,S (v) + O(u 2 ) , (3.10)
where the collinear conformal block f ∆,S (v) describes the small u limit of (3.6)

f ∆,S (v) = g ∆,S (0, v) = v -1 2 S 2 F 1 ∆ + S 2 , ∆ + S 2 , ∆ + S; 1 -v . (3.11)
The first term on the right-hand side of (3.10) with S = 0 corresponds to the Konishi supermultiplet. It gives the leading asymptotic behaviour of A(u, v) at short distances, x 12 → 0, or equivalently u → 0 and v → 1.

According to (3.7) and (3.8), the instanton correction to A long (u, v) takes the form

A inst (u, v) = A (1) inst + g 2 4π 2 A (2) inst + O(g 4 ) , (3.12)
where A

inst ∼ u 2 v 2 D4444 is given by the second term on the right-hand side of (3.8) and A

(2) inst is the first subleading correction that we shall discuss in a moment. Since A

inst (u, v) scales as O(u 2 ln u) in the light-cone limit, it does not affect the leading asymptotic behaviour (3.10). This is in agreement with the known fact that the scaling dimensions and the OPE coefficients of the Konishi and twist-two operators do not receive instanton corrections at leading order, see e.g. [START_REF] Intriligator | Bonus symmetry and the operator product expansion of N=4 SuperYang-Mills[END_REF] and [START_REF] Arutyunov | Perturbative and instanton corrections to the OPE of CPOs in N=4 SYM(4)[END_REF] 8 .

In the previous section we have shown that, among all twist-two operators, only those belonging to the Konishi supermultiplet receive O(g 2 ) instanton corrections to their OPE coefficients and O(g 4 ) corrections to their scaling dimensions. Together with (3.10) this allows us to fix the small u behaviour of the subleading instanton corrections to (3.12):

A (2) inst (u, v) = a (inst) 2,0 uv 2 f 6,0 (v) + O(u 2 ) , = 30 a (inst) 2,0 uv 2 3 -3v 2 + v 2 + 4v + 1 log v (v -1) 5 + O(u 2 ) , (3.13)
where we neglected O(g 4 ) corrections to the scaling dimension ∆ S=0 = 2 + γ K , since they do not enter at this order. Here a 8 Although the claim in [START_REF] Intriligator | Bonus symmetry and the operator product expansion of N=4 SuperYang-Mills[END_REF] is that this is true in general, their general arguments are only valid in the semiclassical approximation. 9 The OPE coefficient for the Konishi operator in the free theory is a (0)

2,0 = 1/6, but we will not need its explicit value in this section.

The ratio a (inst) 2,0 /a (0) 2,0 coincides with the instanton correction to (C K /C (0) K ) 2 (see Eq. (2.42)). Note that A

(2) inst (u, v) does not contain O(u log u) terms. This should of course be the case, as there are no instanton corrections to the anomalous dimensions of twist-two operators at this order. Such corrections first appear at order O(g 4 ). In the next subsection we use crossing symmetry of A inst (u, v) together with the small u behaviour (3.13) in order to compute instanton corrections to certain higher spin operators.

Crossing symmetry and higher spin operators

Before proceeding, let us make an important comment. As follows from the light-cone asymptotic behaviour of A [START_REF] Beisert | Review of AdS/CFT Integrability: An Overview[END_REF] inst (u, v) ∼ u 2 log u, the anomalous dimension of operators of twist four and higher do receive instanton corrections at the leading O(g 0 ) order [START_REF] Arutyunov | Perturbative and instanton corrections to the OPE of CPOs in N=4 SYM(4)[END_REF]. The conformal partial wave analysis shows [START_REF] Alday | Lessons from crossing symmetry at large N[END_REF] that only operators with spin zero receive such corrections. We show below that the situation is very different at order O(g 2 ).

Let us examine A inst (u, v) in the double light-cone limit u, v → 0. Combining the leading asymptotics (3.13) with the crossing relation A inst (u, v) = A inst (v, u) we infer that in the small u, v limit A inst (u, v) should contain the following term

A inst (u, v) = -30 g 2 4π 2 a (inst) 2,0 v u 2 log u + • • • . (3.15)
Following [START_REF] Komargodski | Convexity and Liberation at Large Spin[END_REF][START_REF] Fitzpatrick | The Analytic Bootstrap and AdS Superhorizon Locality[END_REF] we then try to answer how to get such asymptotics from a conformal partial wave expansion. The crucial observation is that, due to the presence of u 2 log u term, this must come from twist-four operators with anomalous dimension γ 4,S = O(g 2 e 2πiτ + e -2πiτ ). More precisely, the expansion (3.5) should contain a contribution from a tower of twist four operators with Lorentz spin S such that

v 2 S=0,2,••• a 4,S u 2 u γ 4,S /2 f 8+S,S (v)| u 2 log u = -30 g 2 4π 2 a (inst) 2,0 vu 2 log u + O(v 2 ) , (3.16)
or equivalently

1 2 S=0,2,••• a 4,S γ 4,S f 8+S,S (v) 1/v = - 30 v g 2 4π 2 a (inst) 2,0 . (3.17)
It is important to emphasise that, as opposed to twist-two operators, for a given spin S the sum on the left-hand side receives the contribution from many twist-four operators. To simplify formulae we do not add an additional index to distinguish such operators.

A very important point about (3.17) is that, given that each term on the left-hand side diverges only logarithmically as v → 0, we need an infinite number of them in order to reproduce the power law divergence 1/v on the right-hand side of (3.17). Furthermore, the divergence will come from the region with large spin, S ≫ 1. The corresponding OPE coefficients a 4,S in the free theory were found in [START_REF] Dolan | Conformal partial wave expansions for N=4 chiral four point functions[END_REF]. In the large spin limit they reduce to

a 4,S = a (0) 2,0 √ π(N 2 -1)S 5/2 2 S+8
+ O(g 2 ) .

(3.18)

The leading asymptotic behaviour of the left-hand side of (3.17) for v → 0 can be computed following [START_REF] Alday | From correlation functions to Wilson loops[END_REF] (see footnote 19 there). Matching the leading 1/v terms on both sides of (3.17) we obtain

γ 4,S = - α S 2 , α = 480 N 2 -1 g 2 4π 2 a (inst) 2,0 a (0) 2,0 . (3.19)
Replacing a (inst) 2,0 with (3.14), we arrive at our final expression for the large spin behaviour of the anomalous dimension of twist-four operators

γ 4,S = 864 κ N (N 2 -1) 2 g 2 4π 2 (e 2πiτ + e -2πiτ ) 1 S 2 .
(3.20)

We remind that since twist-four operators are degenerate, this anomalous dimension should in principle be understood as an average weighted by the tree-level OPE coefficients.

The following comments are in order. Firstly, the large spin asymptotics γ 4,S ∼ 1/S 2 is consistent with the expected behaviour of anomalous dimensions of double trace operators [START_REF] Alday | Comments on operators with large spin[END_REF]. Moreover, γ 4,S is suppressed by the factor of (N 2 -1) as compared with the anomalous dimension of the Konishi operator (see Eq. (2.42)) which is also a characteristic feature of double trace operators. Secondly, while at order O(g 0 ) only twist four operators with zero spin receive instanton corrections, at order O(g 2 ) operators with arbitrarily high spin do. It would be very interesting to compute these corrections directly. Finally, the relation (3.20) describes the contribution of one-(anti)instanton in the SU (N ) gauge group. Making use of (2.43), it is straightforward to generalize it to multi-instantons in the large N limit.

Conclusions

In the present paper we have computed, in the semi-classical approximation, instanton corrections to various correlation functions, involving the half-BPS operator O 20 ′ and the Konishi operator K. Our main results are the explicit expressions (2.42) and (2.43) for the leading instanton contribution to the anomalous dimension of the Konishi operator, as well as for the OPE coefficient of the Konishi operator with two half-BPS operators.

In addition, we considered twist-two operators of general spin S, and showed that the only operators that receive the leading instanton corrections are those which carry spin S = 2 and belong to the Konishi supermultiplet. Using this information, we derived the asymptotic lightcone behaviour of the correlation function of four half-BPS operators and, then, employed the crossing symmetry to determine the instanton contribution to the anomalous dimension of twist-four operators, in the limit of large spin.

Our computations show a very interesting interplay between semi-classical vs quantum instanton corrections and the symmetries of N = 4 SYM. An instance of this interplay arises when comparing instanton corrections to the scaling dimensions of Konishi operator K and its dimension-four supersymmetric descendant K ′ ∼ δ 2 Q δ 2 QK . Superconformal symmetry implies that these two operators should have the same anomalous dimensions. On the other hand, while the leading non-vanishing instanton correction to ∆ K comes from a semi-classical computation, from the perspective of ∆ K ′ computing the same correction would require going through a highly nontrivial analysis of quantum fluctuations! Something similar happens when considering the four-point correlation function of half-BPS operators. Subleading instanton corrections to this correlator involve including quantum fluctuations, but on the other hand, the asymptotic behaviour of such corrections at short distances and on the light-cone is controlled through the OPE by two-and three-point correlation functions, in which instanton corrections at the same order in the coupling constant can be derived from a semi-classical computation. All this seems to hint at existence of some hidden structure underlying instanton corrections in N = 4 SYM.

There are several directions in which this work can be extended. As was mentioned in the Introduction, the spectrum of the dilatation operator in N = 4 SYM should be invariant under the S-duality. Viewed as functions of the coupling constant, the scaling dimensions of operators carrying the same quantum numbers with respect to global symmetries cannot cross each other and they should be invariant under modular transformations. Since modular invariant functions independent on θ-angle ought to be constant, the scaling dimensions should have a nontrivial dependence on θ. One of the direct consequences of our study is that it opens up the possibility to construct such functions in N = 4 SYM by taking into account the instanton corrections. Our results represent the first explicit calculation of the instanton correction to Konishi operator which is the lowest unprotected operator at weak coupling. The obtained expression (2.42) can be thought of as representing the first term in the expansion of the modular invariant function ∆ K (τ, τ ) at weak coupling. It would be interesting to try to determine the modular properties of ∆ K for an arbitrary coupling. As a first step in this direction, one can use the results of this paper to improve an interpolating procedure proposed in [START_REF] Beem | Resummation and S-duality in N=4 SYM[END_REF][START_REF] Alday | Modular interpolating functions for N=4 SYM[END_REF].

Another interesting question concerns S-duality properties of the OPE coefficients. In distinction with the scaling dimensions, they do not have to satisfy von Neumann-Wigner non-crossing rule [START_REF] Korchemsky | On level crossing in conformal field theories[END_REF] and, as a consequence, they may transform nontrivially under the modular transformations of the coupling constant. In general, the properties of structure constants under S-duality are poorly understood. Our results may provide the first hints in this direction, for the simplest case, corresponding to the OPE coefficient of two half-BPS operators and a unprotected operator. For the Konishi operator the leading instanton correction is given by (2.42), finding higher order corrections is an open problem.

It would also be interesting to understand better the interplay between semi-classical and quantum instanton corrections mentioned above. Our results seem to hint at unexpected simplifications when considering quantum corrections around instanton backgrounds in N = 4 SYM. In addition to this, we demonstrated that the semiclassical result for the Konishi operator combined with the crossing symmetry leads to a definite prediction for the instanton correction to the scaling dimensions of twist-four operators with large spin. Computing this correction directly would require to include a quantum fluctuations. This remains a largely unexplored subject. where x i0 = x ix 0 . Using Schwinger parameterization together with the Symanzik star formula we get

D ∆ 1 ∆ 2 ∆ 3 ∆ 4 = 2K ∞ 0 i dt i t ∆ i -1 i e -i<j t i t j x 2 ij = K 4 i<j dδ ij 2πi Γ(-δ ij )(x 2 ij ) δ ij , (B.2)
where K = π 2 Γ 1 2 ∆ i -2 /(2 i Γ(∆ i )) and the integration in the last relation goes parallel to the imaginary axis. The variables δ ij are not independent and satisfy the relations

δ ij = δ ji , 4 j =i δ ij = -∆ i , (B.3)
leaving only two variables independent. It is convenient to choose the latter as δ 12 = j 1 and δ 23 = j 2 . In this way, we obtain where the D-function only depends on conformal cross-ratios (3.4)

D ∆ 1 ∆ 2 ∆ 3 ∆ 4 = K (x 2 13 ) ∆ 1 (x 2
D∆ 1 ∆ 2 ∆ 3 ∆ 4 = i∞ -i∞ dj 1 dj 2 (2πi) 2 u j 1 v j 2 Γ(j 1 + j 2 + ∆ 2 )Γ(j 1 + j 2 + ∆ -∆ 4 )
× Γ(-j 1 )Γ(-j 2 )Γ (-j 2 + ∆ -∆ 2 -∆ 3 ) Γ (-j 1 -∆ + ∆ 3 + ∆ 4 ) , (B.5) with ∆ = 4 i=1 ∆ i /2. Here the integration contours are chosen in such a way that the poles generated by the product of gamma-functions in the first and second lines are located on the different sides.

  .35) The light-ray operator O(z) serves as a generating function of twist-two operators, which are generated by expanding in powers of the light-cone separation z. The dots on the right-hand side of (2.34) denote the contribution of the conformal descendants of the form z ℓ (n∂) ℓ O S (0). To find instanton corrections to the three-point function O 20 ′ (1)O 20 ′ (2)O S (0) , we can first evaluate O 20 ′ (1)O 20 ′ (2)O(z) and, then, apply

9 a

 9 correction to the OPE coefficient of the Konishi operator. It can be written in terms of the free-theory coefficient

  24 ) ∆ 2 D∆ 1 ∆ 2 ∆ 3 ∆ 4 (u, v) , (B.4)

	x 2 13 x 2 34	∆ 1 -∆ 3 2	x 2 14 x 2 34	∆ 2 -∆ 4 2
	x 2 14		x 2 13

It remains unclear, however, whether the scaling dimensions of unprotected operators receive nonzero corrections running in powers of e

2πi(τ -τ ) due to the contribution of instanton-anti-instanton configurations. For some evidence on the string theory side see[START_REF] Green | Duality and higher derivative terms in M theory[END_REF]. We thank Nick Dorey and Pierre Vanhove for discussions about this point.

We recall that the anti self-dual part of the instanton strength tensor vanishes, F α β = 0.

Indeed, as we show in the next subsection, the correlation functions involving the operator K ′ do vanish in the semiclassical approximation.

Indeed, tr(φ[START_REF] Montonen | Magnetic Monopoles as Gauge Particles?[END_REF],AB φ (6),CD ) is proportional to ζ 8 ǫ ABCD and, therefore, O(8)20 ′ ∼ Y AB Y CD ǫ ABCD = 0.

Here we took into account that the two-point correlation function of K does not receive instanton correction at order O(g 2 e 2πiτ ), see(2.13).

Strictly speaking, this relation holds for the so-called minimal correlation functions[START_REF] Green | Instanton induced Yang-Mills correlation functions at large N and their AdS(5) x S**5 duals[END_REF], for which the product of n operators soaks up the sixteen exact fermion modes. It is straightforward to verify that the correlation functions (2.38) and (2.39) are indeed minimal.

Acknowledgments

We are grateful to Agnese Bissi, Nick Dorey, Jaume Gomis, Stefano Kovacs and Pierre Vanhove for useful discussions. L.F.A. and G.P.K. acknowledge Nordita and Perimeter Institute, respectively, for hospitality where part of this work has been done. The work of L.F.A. was supported by ERC STG grant 306260. L.F.A. is a Wolfson Royal Society Research Merit Award holder. This work of G.P.K. was supported in part by the French National Agency for Research (ANR) under contract StrongInt (BLANC-SIMI-4-2011).

A Conventions

In this appendix we summarize our conventions for dealing with spinor indices. For an arbitrary four-dimensional vector x µ we associate a 2 × 2 matrix

where σ µ = (1, σ) is defined in terms of Pauli matrices. The indices are raised and lowered with a help of the completely antisymmetric tensor

with ǫ αβ ǫ αγ = δ γ β and ǫ α β ǫ α γ = δ γ β . Using these definitions we find

For derivatives we have in the similar manner

The gauge field strength tensor F µν can be decomposed into (anti) self-dual part as

B Definition of the D-functions

The D-functions are defined as