\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$C_2$$\end{document}C2 generalization of the van Diejen model from the minimal \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(D_5,D_5)$$\end{document}(D5,D5) conformal matter

We study superconformal indices of 4d compactifications of the 6d minimal \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(D_{N+3},D_{N+3})$$\end{document}(DN+3,DN+3) conformal matter theories on a punctured Riemann surface. Introduction of supersymmetric surface defect in these theories is done at the level of the index by the action of the finite difference operators on the corresponding indices. There exist at least three different types of such operators according to three types of punctures with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A_N, C_N$$\end{document}AN,CN and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( A_1\right) ^N$$\end{document}A1N global symmetries. We mainly concentrate on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$C_2$$\end{document}C2 case and derive explicit expression for an infinite tower of difference operators generalizing the van Diejen model. We check various properties of these operators originating from the geometry of compactifications. We also provide an expression for the kernel function of both our \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$C_2$$\end{document}C2 operator and previously derived \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A_2$$\end{document}A2 generalization of van Diejen model. Finally, we also consider compactifications with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A_N$$\end{document}AN-type punctures and derive the full tower of commuting difference operators corresponding to this root system generalizing the result of our previous paper.


Introduction
Intricate connection between supersymmetric gauge theories and integrable systems plays an important role in modern theoretical physics and mathematics. From physics point of view whenever some sector of the gauge theory is related to an integrable system, observables in this sector can be computed using rich variety of integrability techniques. Canonical example of this situation is the integrability of N = 4 Super Yang-Mills (SYM) theory [1,2]. In particular it connects calculation of the planar scaling dimensions of N = 4 SYM to the integrable system of spin chains and allows to compute these scaling dimensions at arbitrary coupling. On the other hand exploration of such connections can also shed light on some questions about integrable systems and even lead to construction of new classes of such systems making this kind of studies interesting from mathematics point of view.
In our work we are exploring a particular class of connections between six dimensional (1, 0) superconformal field theories and elliptic quantum mechanics Hamiltonians in the spirit of Bethe/gauge correspondence of Nekrasov and Shatashvili [3][4][5][6]. In particular we consider four-dimensional theories with four supercharges obtained by compactification of a 6d SCFT on a punctured Riemann surface. In order for an integrable model to emerge in this setting we have to introduce surface defects into our 4d theory and study its superconformal index.
This construction was first established in the context of compactifications of 6d (2, 0) of ADE type [7,8]. In particular the corresponding 4d superconformal indices with the defect were found to be closely related to the Ruijsenaars-Schneider (RS) elliptic analytic finite difference operators (A∆Os). Later these results were extended to many other cases: class S k 4d theories [9], compactifications of A 2 and D 4 minimal 6d SCFTs [10] and compactifications of rank one E-string theories [11]. In some of these cases the obtained A∆Os were already known in the literature. For example in E-string compactifications van Diejen (vD) model [12,13] was observed. But in some cases operators were previously unknown as in the case of the minimal 6d SCFTs compactifications. Study of such novel operators constitutes an interesting field of research with some initial steps already taken in this direction [14,15].
From the point of view of physics the connection of superconformal indices in the presence of defects with the integrable quantum mechanics Hamiltonians allows one to bootstrap index of an arbitrary theory obtained by compactifying corresponding 6d theory on a punctured Riemann surface. In particular if we know the eigenfunctions of the corresponding A∆O we can compute the index of any theory obtained in such compactifications including non-Lagrnagian theories for which there are no other methods of computing the index. As mentioned previously, the first setting where these ideas were tested are 4d class-S theories obtained in the compactifications of 6d (2, 0) theory [7]. In this case the corresponding integrable system was given by elliptic RS model. Eigenfunctions of the elliptic RS Hamiltonians are not known in general but some of their limits are well studied in the mathematical literature. These limits were used in order to prove previously established relations [16][17][18] of superconformal indices with the Macdonald and Schur polynomials allowing one to compute indices of class S non-Lagrangian theories.
The construction outlined above relies on the intermediate 5d layer. For things to work there should exist an effective 5d gauge theory obtained by compactifying original 6d SCFT on a circle with a choice of holonomies for its global symmetries. In particular different 5d compactifications lead to different types of punctures on the Riemann surfaces used to obtain 4d theories. One of the interesting problems related to this fact is the compilation of the dictionary between known compactifications of various 6d SCFTs and elliptic integrable systems. In our previous paper [19] we have considered compactifactions of the 6d minimal (D N +3 , D N +3 ) conformal matter theories [20,21]. In particular we have derived A∆Os corresponding to the intermediate 5d SU(N + 1) gauge theory or equivalently the A N type puncture on the compactification surface. These elliptic A∆Os appeared to be previously unknown A N generalizations of vD model. On the other hand there are at least two more 5d effective descriptions of 6d SCFTs corresponding to USp(2N ) and SU(2) N gauge theories giving rise to the punctures with the same C N and (A 1 ) ⊗N global symmetries. These two descriptions should lead to additional higher-rank generalizations of the vD model.
In our present paper we follow this line of research and closely study compactifications of the minimal (D N +3 , D N +3 ) conformal matter theory on a Riemann surface with the C Ntype punctures. In particular we concentrate on the next to simplest case of N = 2 1 and derive corresponding infinite tower of A∆Os. We also devote part of the paper to the study and proof of some remarkable properties of these operators that follow from the geometry of corresponding compactifications. In particular we prove the novel kernel property for two operators of different type but same rank. We also discuss commutation property of obtained operators.
The paper is organised as follows. In section 2 we review our previous results for A N type operators. In addition to the previous results we also derive full tower of such A∆Os which was not obtained previously. In Section 3 we derive in details novel generalization of vD operators corresponding to C 2 root system. In section 4 we discuss properties of derived C 2 A∆Os. In particular we pay special attention to a new kernel function for simultaneously A 2 and C 2 operators and give fully analytic proof of the corresponding kernel equation. We also briefly discuss commutation relations and check them perturbatively in expansion. In section 5 we briefly summarize our results and discuss plans for the future research further developing these results. Finally the paper has a number of Appendices collecting useful formulas as well as technical details of various calculations in our paper.
2 A N operators.
In this section we will briefly review derivation of the A N generalization of the van Diejen operator. So called basic version 2 of this operator has been derived in our previous paper [19]. Here we will extend this result to the full tower of operators.
Just as in [19] in order to derive A∆O we start with the 4d three-punctured sphere theory which was first introduced in [21]. This theory is obtained in the compactification of the 6d minimal (D N +3 , D N +3 ) conformal matter on a sphere with two maximal SU(N +1) punctures and one minimal SU(2) puncture. The quiver of this theory is shown in Figure 1. In addition to the global symmetry, the punctures are characterized by the moment map operators. In particular both SU(N + 1) maximal and SU(2) minimal punctures are characterized by (2N + 4) mesonic and 2 baryonic moment maps: where a i , u, v, w are fugacities of the Cartans of the 6d global SO(4N + 8) symmetry. Subscripts of the moment maps written above denote their charges w.r.t. to these symmetries. Further S-gluing two such trinions along the maximal punctures we can obtain fourpunctured sphere with zero flux two maximal and two minimal punctures. To obtain the corresponding 4d gauge theory we should just take two copies of trinion theories shown in Figure 1 and then identify and gauge corresponding global symmetries of the maximal (a)A N trinion theory. punctures. Here and everywhere else in the paper all operations with gauge theories are expressed in terms of the superconformal indices. In this langauge, the gluing procedure takes the following form: where K A 3 (x, y, z) is the index of the trinion with y being fugacity of the global SU(N + 1) symmetry of the puncture we glue along andK A 3 is the index of the conjugated trinion. Finally κ N is the usual constant given by: Performing this S-gluing operation we obtain the four-punctured sphere theory shown in Figure 1 with the corresponding superconformal index specified in (B.1). This theory was previously obtained by us in [19] to derive basic A N A∆Os. Now in order to derive the operator we should close two minimal punctures of this four-punctured sphere. To do it we should break the global symmetry of the puncture. This can be achieved by giving a non-trivial vev ∂ L 12 ∂ K 34 M = 0 to the derivatives of one of the moment map operators. When we close punctures with at least one of K or L not equal to zero, i.e. vev is space-time dependent, we effectively insert defect into the theory [7,22]. At the level of the superconformal index, closing the puncture amounts to giving a corresponding weight to the fugacity of the puncture's global symmetry. Once we do it we hit a pole of the index. Then computation of the residue of this pole results in the superconformal index of the IR theory that the UV theory flowed to due to the introduction of the vev. In our case we choose to close SU(2) z minimal puncture with the defect and SU(2)z puncture without, i.e. choosing L = K = 0 in the vev. In particular let's say that we are going to compute A∆O acting on the puncture with the moment maps of chargesh i and an overall U(1) chargẽ For example the moment maps of the SU(N + 1) x puncture of the trinion shown in Figure 1 correspond to M u operators specified in (2.1) and have the following charges: where we have flipped the chargeh 2N +6 of the last moment map since this parametrization will be more natural for our operators and in this case all of the moment maps of the maximal puncture transform in the fundamental representation of SU(N + 1) x . We can now notice from (2.1) that charges of the minimal puncture moment maps are related to the charges of the moment maps of the maximal punctures by simple relatioñ So we can express everything in terms of only the moment maps of the maximal puncture we act on. Now assume we give vevs to the moment maps (mesonic or baryonic) with charges h −1 ih 1/4 of both SU(2) z and SU(2)z. Moment maps we use should be the same in order to keep total flux of the 6d global symmetries zero. At the level of the index calculations it corresponds to computing the residue of the index of the four-punctured sphere theory at the pole where L, M,L,M are positive integers corresponding to the powers of derivatives inside the vev. As we mentioned previously it is enough to introduce defect only for one of the two punctures. Hence we chooseM =L = 0 and keep L , M general. Then we compute corresponding residues of the index of the four-punctured sphere theory and obtain theory for the tube with two maximal SU(N + 1) punctures and a codimension-two defect. Its superconformal index is given by 3 : Finally in order to obtain desired A∆O we glue our tube with the defect to an arbitrary Riemmann surface with maximal SU(N + 1)x puncture. As the result of this gluing we expect to obtain action of a finite-difference operator on the index I(x) of this 4d N = 1 theory: Details of all the calculations summarized above are given in Appendix B. They result in the following operator: where C l, m, r, s L,M are x-independent constant factors given by, Also ∆s are shift operators defined as follows: The operator contains all shifts of the form q m i p l i x i where m and l are all possible partitions of length N + 1 of M and L correspondingly. At each level, i.e. fixed M and L there are 2N + 6 operators due to 2N + 6 moment maps with the chargesh kh − 1 4 . There are also 2N + 6 other operators obtained by giving vevs to the flipped moment maps of the chargeh −1 kh 1 4 . They have similar form and properties so we do not present them here. Now we will refer to the case M = 1 and L = 0, or vice versa, as basic operators. These basic operators for A N generalizations of vD model were derived by us previously in [19]. The operators above also reproduce our previous results when we fix M = 1, L = 0 or M = 0, L = 1. Further in our paper we will also need the basic operator obtained by closing flipped moment maps. This operator can be found in [19] and has the following form: where the shift part of this operator is given by (2.14) and the constant part is given by: This constant part is elliptic function in each x i variable with periods 1 and p. It has poles in the fundamental domain at the following positions: where we defined 3 Derivation of C 2 operators.
In this section we will discuss derivation of the C-type rank-2 analytic finite difference operators A∆O. In our previous paper [19] we have already derived the basic operator for A N generalization of van Diejen model and the calculation of the full tower of these operators is summarized in Section 2. Unfortunately full C N generalization is still out of reach for us due to technical complications but we can concentrate on the study of rank-2 generalizations with C 2 root system. We will derive C 2 operator corresponding to the insertion of the codimension-two defect due to the nontrivial vev of the holomorphic derivatives of the moment maps ∂ K ± M = 0. To simplify our calculations of A∆O we will start with the derivation of four-and three-punctured spheres with C 2 type punctures.
The basic C N trinion was already derived by us in Appendix B.5 of [19] and is shown in Figure 2(b). Generalizations of this trinion to the cases of three maximal punctures and higher numbers of minimal punctures can also be found in [23]. In order to derive this trinion we used a tube theory with one USp(2N ) and one SU(N + 1) maximal punctures shown in Figure 2(a). This theory was first discussed in [24]. In order to obtain trinion theory shown in Figure 2(b) we start with the three-punctured sphere with two maximal SU(N + 1) and one minimal SU(2) punctures and glue two A N C N tubes to the maximal punctures. The moment maps of the punctures of the resulting trinion are given by (2N + 6) mesons both for maximal and minimal punctures: where w is one of the parameters of Cartans of the global SO(4N + 8) symmetry of the 6d minimal conformal matter theory. In general we will use two types of parametrization of 6d global symmetry in this paper. First one is natural to use in compactifications with C N type punctures and it is given bỹ Another parametrization is useful when we work with A N expressions and is given by: These two parametrizations are related by the following map that we will often need in our calculations: Now in order to derive finite difference operators we should proceed in the same way as in the case of A N operators summarized in Section 2. We start by taking two C N trinion theories T A x,y,z where A is the non-zero flux of the trinion and x, y, z in the subscript are fugacities of the global symmetries of two maximal and one minimal punctures correspondingly. Then we take an arbitrary N = 1 theory obtained in the compactification of the minimal (D N +3 , D N +3 ) conformal matter theory on the Riemmann surface of genus g with s punctures denoted as C g,s with at least one USp(2N )x maximal puncture. Gluing all three surfaces together along the maximal punctures results in a Riemmann surface of the same genus and global symmetry fluxes but two extra minimal punctures. We will be performing all these operations at the level of superconformal indices. There gluing amounts to identifying global symmetries of the punctures we glue and gauging it. So for the indices we can write down the following identity after gluing: where I [C g,s ] is the index of a theory obtained by compactifications of 6d theory on the Riemmann surface C g,s and K C 3 (x, y, z) is the index of the trinion theory shown in Figure 2(b) whileK C 3 (x, y, z) is its conjugate. Definition and properties of the elliptic Γ-function Γ e (x) are given in Appendix A. Then, closing SU(2) z and SU(2)z minimal punctures with or without introduction of defects we can obtain finite difference operators. As discussed previously in order to close minimal punctures we have to give certain weights to corresponding global symmetry fuagcities z 1 and z 2 . As a result we get the follwoing identity lim where O is some finite difference operator, and in case we close punctures without introducing defects we just obtain the identity operator.
Here for convenience we will not compute the full expression (3.5) directly right away. Just as for A N operators derived in Section 2 it appears to be much easier to perform calculation in a slightly different way. First we derive the index of the four-punctured sphere with zero flux, two maximal USp(2N ) and two minimal SU(2) punctures. For this purpose we perform S-gluing of two trinions T AC x,y,z with SU(N + 1) and USp(2N ) maximal punctures and SU(2) minimal punctures each. This kind of trinions can be derived by appropriate gluing of A N C N tube theory shown in Figure 2 (b) Trinion with two maximal USp(2N ) and one minimal SU(2) puncture. in [21]. Superconformal index of the resulting four-punctured sphere theory is given by: where K AC 3 (x, y, z) is the index of the trinion with y being SU(N + 1) fugacity and κ N is the usual constant given in (2.3). Geometrically this operation is shown in Figure 5(a).
Next we close on of the two minimal punctures by giving nonzero vev to one of the M z moment maps given in (3.1). At the level of the index computations this means that the corresponding global symmetry fugacity should be consistent with the vev, i.e. in our case we should fixz fugacity tõ whereã i can be chosen arbitrary since all the expressions are symmetric w.r.t. permutations ofã i . Integers K and M correspond to an order of the derivative of the moment map in 34 and 12 planes correspondingly, i.e. we give vev to ∂ M 12 ∂ K 34 M i = 0. Physically this corresponds to introducing various codimension two defects into 4d theory. For the first SU(2)z minimal puncture we choose to close it without introducing any defect, which in turn corresponds to K = M = 0 choice in (3.8). At this value ofz, the superconformal index (3.7) of the four-punctured sphere has a pole. Computing the residue at this pole we obtain the index of the three-punctured sphere: Here subscript (3; i, 0) refers to the fact that we obtain three-punctured sphere by closing minimal puncture of the four-punctured sphere choosingã i for the vev in (3.8) and not introducing a defect, which corresponds to the choice K = M = 0 in the same equation. Calculation of (3.7) and (3.9) results in the theory shown in Figure 3. Detailed derivations of this section can be found in Appendix C. Now that we have obtained the desired three-punctured sphere theory we are ready to derive A∆O. For this purpose we should close the remaining minimal SU(2) z puncture by giving nontrivial vev ∂ M 12 ∂ K 34 M = 0 to a derivative of one of its moment maps. In order to get zero total flux through the resulting two-punctured sphere we should choose the same moment map as we did closing SU(2)z minimal puncture previously. The difference is that now we have to choose non-trivial derivative that is at least one of K and M numbers is not zero. Physically this corresponds to introducing codimension two defect into the tube theory with two maximal punctures as shown in Figure 5(b). At the level of the index we should give the following value to the z-fugacity: Then, according to (3.6), capturing the corresponding pole and gluing the resulting tube with two maximal USp(2N ) punctures to an arbitrary theory with at least one puncture of this type as shown in Figure 5(c), we obtain an A∆O. However, technically it sometimes appears not to be as straightforward. In particular if we perform these operations with the expression (C.12) instead of A∆O we obtain some integral-finite difference operator. It is highly possible that in fact this operator can be written in the form of A∆O. However technical issues make it too difficult task and we leave it for future investigation.
Here we will instead concentrate on the derivation of next to the lowest rank C 2 operator. This case is simpler to analyze. For example, in our trinion theory shown in Figure 3, when N = 2 the AS multiplet of the right SU(3) gauge node becomes just antifundamental, and a chain of duality transformations can be used to simplify the theory. These calculations are summarized in Appendix C. As a result we obtain single-node SU(6) gauge theory shown in Figure 4 with the corresponding index specified in (C.13). Finally, we can close SU(2) z puncture. In particular we give the following weight consistent with (3.8) to the fugacity z: Performing this closure we obtain the index of the theory for two punctured sphere with two USp(4) maximal punctures The index itself is specified in (C.18). This time it does not have natural gauge theory interpretation in case of general K due to the presence of the codimension two defect. As a final step of our derivation, we glue the obtained tube with the defect to an arbitrary N = 1 theory with at least one maximal USp(4) puncture and obtain A∆O as follows: Closing minimal punctures of the four-punctured sphere. x Gluing tube with the defect to an arbitrary Riemann surface. Details of this calculation can be found in Appendix C. It leads to the following expression for A∆O: where the constantC k K is given bỹ We have also introduced notation h i , i = 1, . . . , 10 for U (1) charges of the moment maps we act on. Notice that according to (3.1) moment maps of both minimal and maximal punctures have the very same charges, so that in fact for both charges of the maximal punctures we act on and minimal punctures we close. We have also introduced here an overall U (1) charge h. The operator itself is labelled by an index k according to the choice of U (1) charge h k of the moment map we give vev to in order to close the minimal puncture and to introduce defects into the theory. The first sum in the expression is performed over all possible partitions K = (k 1 , · · · , k 6 ) of the integer K. Finally the products in (3.14) should be understood as ordered ones, i.e., one should assume n 2 i=n 1 is just 1 if n 2 < n 1 since there are no terms in the product. The coefficientC k K is chosen so that the x-independent factors of the highest order shift terms ∆ ±K q (y 1,2 ) are just one. Notice that the expression above is not the most general since we could have closed SU(2) z minimal puncture using general K and M integers in (3.10), but we will concentrate here only on the case of non-trivial K while putting M = 0. As we see the operator appears to be quite complicated and it is hard to analyze it. Instead it is useful to write down basic operator corresponding to the choice of K = 1, i.e. the simplest non-trivial case. Taking into account six possible partitions K we can write down explicit form of A∆O: One thing to be noticed is that just as in the case of vD model, the constant term choice is not unique. What really defines this constant term are the following properties. First of all it can be noticed that W (h k ;1,0) (x, h) is an elliptic function in both x 1 and x 2 with periods 1 and p. In the fundamental domain this elliptic function has poles located at: which are in fact exactly the same as the poles of vD model. Corresponding residues are given by: This concludes our derivation of the C 2 A∆O. All possible details of these derivations can be found in the Appendix C.

Properties of the C 2 operators.
There is a number of interesting properties that the operators we derived should posses by construction. In this section we discuss checks and, where it is possible, give proofs of these properties using explicit expressions we have derived.
The most important but also most complicated property to discuss is the so-called kernel property of our operators. The main idea behind it is that the superconformal index of any N = 1 theory obtained in the compactifications of 6d minimal (D 5 , D 5 ) conformal matter theory is a kernel function of our operators. The kernel function is defined by the following mathematical identity: (4.1) Physically we consider the superconformal index of the theory obtained in the 6d compactification on the Riemmann surface C g,s [z, u] with at least two maximal punctures parametrized by the fugacities z and u correspondingly. In general there can be as many punctures as we want. The claim here is that any such superconformal index play the role of the kernel function for the derived A∆Os according to (4.1). Namely we can act with our operators on different punctures and we should always obtain the same result. This property is expected to hold due to an argument coming from the geometry of compactification shown in Figure  6. Equivalently it can be understood from the invariance of the superconformal indices under S duality transformations. In principle we can choose any sphere with multiple punctures to check the kernel property. Natural simplest candidate would be WZW model for the two-punctured sphere with two maximal USp(2N ) punctures shown in Figure ??. However, we will check a more interesting and trickier case where two punctures are of different types. Namely we take the tube theory shown in Figure 2(a) which has one SU(N + 1) and one USp(2N ) maximal puncture. Then for N = 2 the kernel property (4.1) in this case reads: where K AC 2 (x, y) is the index of the tube specified in (C.5) with x and y being fugacities of USp(4) and SU (3) .13) and (3.17) correspondingly. Notice that in order for the kernel property to work we need to close minimal punctures in the same way on two sides of the equation, i.e. we have to give vev to the same moment map of the minimal puncture. In particular in A 2 case we should give vevs to the minimal puncture moment map of the chargeh −1 ih 1 4 whereh i are given in (2.5). At the same time, in C 2 case minimal puncture is closed by giving vev to the moment map of charge h i specified in (3.16). Using the map (3.4) we can see that the only map that matches between the two is If we would like to check the kernel property for other ways of closing minimal punctures we should consider A 2 operators (2.10) with the non-flipped moment maps so that in this case we close the minimal puncture with the vevs of the moment map with chargeh ih − 1 4 . Because of this complication here we consider only closing with the vev given to the moment map of charge (4.3). The proofs for other operators should work identically. Also in our proof we restrict ourselves to only basic operators (3.17) and (2.13), i.e. we fix r = 1 , m = 0 (or equivalently m = 1 , r = 0 on both sides of the kernel equation (4.2). We have to do it since higher operators specified in (3.14) and (2.10) are too complicated for the analysis. However the kernel property (4.2) should also work for these higher operators as well.
Since we precisely know explicit expressions for both operators (2.13) and (3.17) as well as supposed kernel function (C.5) it is straightforward to check this kernel identity. Acting with finite difference operators of two kinds on the kernel we get two algebraic functions. Equality of these functions implies validity of the kernel property (4.2). In Appendix D we give details of the analytic proof of this identity.
Another important property of the derived A∆Os is their commutation with each other which also directly follows from the S duality of our compactification construction as shown in Figure 7. Since it does not depend which duality frame we close the minimal punctures in, it also does not matter in which order two different A∆Os act on the index. Hence we  (2) puncture (shown as red crosses) with the fugacity z. For example in case of the A 2 C 2 tube and kernel property (4.2) one maximal puncture has USp(4) global symmetry while the second one has SU (3) symmetry. Then we close the minimal puncture by giving space dependent vev ∂ m 12 ∂ k 34 M = 0 to one of the moment maps of this puncture. As discussed in the paper this corresponds to the introduction of the codimension-two defect into 4d effective description and results in the action of the A∆O on one of the maximal punctures of the index of the theory corresponding to the Riemann surface with only two maximal punctures. This operation can be performed in different duality frames. In each frame we obtain operators acting on one of the two punctures. Since the result of the calculation should not depend on the duality frame we choose we conclude that the action of two, possibly different operators, on two different punctures leads to two expressions equal to each other. Hence we arrive to the kernel property equation (4.1).
conclude that all of operators we derived should commute with each other: Here h a and h b are moment maps we give vev to in order to close two minimal punctures and K 1,2 , M 1,2 are numbers of holomorphic/antiholomorphic derivatives. Unfortunately it is very hard to prove or even to check these relations for general operators (3.14) of the full tower. However we can perform checks of the commutation relations (4.4) for the basic operators, i.e. for the case when K 1,2 and M 1,2 are taking values 0 and 1. In particular we show that:  (4) puncture (orange disc) with fugacity x and two minimal SU(2) punctures (red crosses) with fugacities z 1 and z 2 correspondingly. Then we introduce codimension two defects into this theory giving space-dependent vevs to some of the moment maps M h1 and M h2 correspondingly. Giving this vev is equivalent to closing puncture and at the level of the index computations leads to an action of two different A∆Os. We can perform this operation in any duality frame and the result should not depend on the choice of a particular frame. Hence the action of operators on the index does not depend on their order. Since the compactification surface, and hence the effective 4d descriptions were chosen arbitrarily we conclude that the operators themselves also commute.
In Appendix E we check these relations. In particular, using ellipticity properties of θ p (x) function we prove that the action of the third type of the commutators on an arbitrary trial function is zero. Hence these commutators are zero themselves. For the first and second type of commutators it is hard to prove analytically that their action on a trial function is zero. Instead we check these identities perturbatively in p and q expansion. These checks performed up to a sufficiently high order suggest that both of the commutators are indeed zero.
5 Discussion and outlook.
In this paper we have considered 4d description of the compactification of the minimal (D 5 , D 5 ) conformal matter theory on a Riemann surface with USp(4) maximal punctures. Using this 4d description and introducing codimension-two defect in corresponding theory we derive C 2 generalization of the van Diejen model. In particular we obtain an infinite set of analytic difference operators acting on the maximal puncture of C 2 type. Different operators of our set correspond to different ways of closing minimal punctures on the compactification surfaces. They are organized in a decuplet of basic operators and an infinite tower of operators on top of each of the basic ones. The operators we obtain are supposed to satisfy interesting and important properties following directly form the geometry of compactifications. First such property is the commutation of all operators. In our paper, using combination of residue computations and perturbative expansions we show that at least all basic operators indeed commute with each other. Second important property is that the superconformal index of any compactification of the minimal (D 5 , D 5 ) conformal matter theory on a surface with several USp(4) punctures is the kernel function of our difference operators. As an example of such kernel function we considered a tube theory corresponding to the compactification on a sphere with one USp(4) and one SU(3) puncture. In order to prove the kernel property we should act on two punctures with two different operators: A 2 operator derived in [19] from one side and C 2 operator derived in the present paper from the other. This fact makes corresponding kernel function very interesting for the study. In the paper we managed to prove this kernel property fully analytically for the basic A 2 and C 2 operators.
Another result of our paper is derivation of the full tower of A N generalization of van Diejen operator orginating from the compactification of the minimal (D N +3 , D N +3 ) conformal matter theory on a Riemann surface with SU(N + 1) maximal punctures. We have already derived basic operators of this kind in our previous paper [19]. However the tower of operators was missing and we filled this gap in the present paper.
The present paper is another brick in our program of establishing dictionary between compactifications of 6d SCFTs and integrable analytic difference operators and studying this dictionary in details. There are plenty of ways this research can be continued in. In particular there should be separate difference operator for each possible 5d compactification of 6d theory, or equivalently for each possible puncture type from 4d prospective. So far several results have been observed in this frame. First of all BC 1 van Diejen model was observed in E-string compactifications in several ways [11,19]. Also previously unknown A 2 and A 3 operators were derived using compactifications of the minimal conformal matter theories of types SU(3) and SO(8) [10]. Finally we have considered compactifications of the minimal (D N +3 , D N +3 ) conformal matter theories on a Riemann surfaces with A N and C 2 (in case of (D 5 , D 5 ) 6d theory) type punctures. However there are plenty of examples of 6d SCFT compactificatoins that are known up to date [25][26][27][28][29][30][31][32][33][34]. Our main goal is to extend our results to more examples of compactifications and ideally establish a dictionary between them and integrable models. Such dictionary can from one point of view shed light on physics of 6d, 5d and 4d supersymmetric gauge theories. From the other point of view this research program can lead to important results in the field of integrable systems.
First of all as a continuation of the present paper it is natural to consider compactifications of the minimal (D N +3 , D N +3 ) conformal matter on Riemann surfaces with other types of punctures. In first place these are of course punctures of C N type for general N . In the present work we failed to derive explicit form of the corresponding A∆O because of the technical difficulties but it is worth trying to overcome them. Second candidate for the study is (A 1 ) N puncture. Corresponding 4d description was derived in [20]. It would be interesting to derive explicit form of the corresponding difference operators and study their properties. For example we immediately know examples of the kernel functions of these operators even before deriving them explicitly.
Second candidate for the studies are non-minimal (D N +3 , D N +3 ) conformal matter theories, which are obtained as worldvolume theories of the stack of k M5 branes probing D N +3 singularity [26]. In this case the known 5d gauge theory, and hence the puncture type on Riemann surface, is a direct generalization of (A 1 ) N case of minimal conformal matter compactification discussed above. It corresponds to the (A k−1 ) 2 × (A 2k−1 ) N × (A k−1 ) 2 gauge theory in 5d and hence the puncture type with the same global symmetry. Trinion theories with two maximal and one minimal punctures, which are building blocks in the construction of A∆Os, are also known [35]. So construction of the correspodning difference operators should be straightforward though can happen to be technically complicated.
Finally one more thing to be made in this direction is the study of the rank Q E-string theory compactifictaions along the same lines of research. Unfortunately so far corresponding trinion theories were not obtained but the tube theories were already derived in [36,37]. Despite we can not derive corresponding A∆O using only tube theory in this case we can naturally conjecture it to be BC n van Diejen model. We can at least test this conjecture by checking corresponding kernel property.
Another possible direction of the research is related to the study of the properties of the operators we derived. Most interesting questions here are related to the eigenfunctions of these operators. Deriving full eigenfunctions is of course difficult and most probably impossible. Realistically we can try to derive certain limits of these eigenfunction. In case we succeed these eigenfunctions can play the same role for indices of N = 1 theories as Macdonald and Schur polynomials played for the indices of N = 2 theories [16][17][18].
It would also be interesting to establish connections of our research program with other integrable models emerging in supersymmetric gauge theories. For example integrable systems derived using compactifications of the worldvolume theory of k M5 branes probing A N −1 singularity were shown to be related to the set of transfer matrices [9,[38][39][40]. It would be interesting to understand if similar connection emerges in case of D N +3 singularity.
Also, a relation of E-string theory to the van Diejen model was observed in [41] using quantization of the corresponding Seiberg-Witten curve. Later, authors of this paper generalized their result to 6d (1, 0) SCFTs with SO(N ) gauge group and (N − 8) fundamental flavors obtaining yet another set of difference operators [42]. It would be interesting to clarify precise relation of our construction with the methods and results of these papers.

Acknowledgments
We would like to especially thank Shlomo Razamat for numerous discussion, collaboration on different stages of this projects and for reading the draft of the paper. AN is also thankful for the hospitality of the Theoretical Physics Division of Uppsala University where parts of this work were done. The work of BN is supported in part by Israel Science Foundation under grant no.

A Special functions
Here we summarize some definitions and properties of special functions used in the paper.
Elliptic Gamma function is defined through the following infinite product: It can be easily seen that the poles of this function are located at the following values of the argument: The following relation will be useful in our calculations: Also we will often deal with the elliptic beta integral formula Here κ is defined to be The Theta function is defined as follows: where (z; p) ∞ is the usual q-Pochhammer symbol defined as follows: Following properties of theta function will be useful to us We will also use the following duality identity from [43]:

B Derivation of higher tower A N operators
In this appendix we will give details of the derivation of A N A∆Os that we summarize in Section 2. We start with the four-punctured sphere theory whose quiver is shown in Figure  1. This theory and its superconformal index were derived in [19]. The index is given by the following expression: This is an index of an SU(N + 3) SQCD with 2N + 6 flavors and a superpotential. Variables x,x, z,z are fugacities of the global symmetry of two maximal SU(N + 1) x,x and two minimal SU(2) z,z punctures correspondingly. We close the two SU(2) punctures by setting z andz to the values in (2.7). Performing calculations for general values of the charges is complicated. So we will perform here the calculation only for one of the charges fixingh i =h 2N +6 = u N w 2 −2N −4 . Corresponding positions of the poles are then given by Also as discussed in Section 2 we will further considerL =M = 0.
The pole inz is an explicit simple pole so we can just take the residue. The pole in z however, is due to contour pinching. To see this clearly, it is useful to first perform a Seiberg duality leading to an index of SU (N + 2) gauge thoery with 2N + 5 flavors, with the followng index, where we have already dropped irrelevant overall factors and substituted values of z from (B.2). The subscript of the index as usually contains all the information of the underlying theory. In particular 2 stands for the number of punctures. Index 2N + 6 stands for the index of the moment map we give vev to in order to close the puncture (moment map of chargẽ h 2N +6h − 1 4 in this case). Finally (L, M ) denotes derivative powers in the vev and are the same as L and M integers in (B.2). As can be seen from the expression above the contour pinching is due to collision of poles at the following values: where m i and l i are partitions of M and L respectively. There are (N + 2)! such poles coming from permutations of t i but they all give the same result and we are ignoring overall factors. It can be checked that the two lines in (B.4) give the same result so we will consider only the first set of poles. After computing the residue we get the following expression for particular partitions l, m, To get the A∆O acting on SU(N + 1) x puncture, we should S glue this tube to the index of an arbitrary theory I(x) with SU (N + 1)x puncture and sum over all partitions, Note that there is a zero in (B.5) coming from Γ e (pq) but this is canceled against the pinching ofx at the following values, up to permutations of x i which give the same result due to the Weyl symmetry of A N root system. If we specifyx i to the values written above we obtain double pole due to the contour pinching of (B.6) but as we mentioned one of them is cancelled by the zero of Γ e (pq).
Computing the residue of the remaining pole we obtain the following contribution of each individual partition: Summing over all partitions and using properties of Γ and θ functions (A.9) we get the difference operator given in (2.10) with theh k =h 10 . Operators obtained using another ways to close minimal punctures can be derived in absolutely identical manner and lead to the same A∆O (2.10).

C Derivation of C 2 operator
In this section we give details of the derivations of results summarized in Section 3 and show in all details how to derive C 2 type generalization of the van Diejen model. Let's start with the derivation of the trinion that has one maximal USp(2N ) puncture, one maximal SU(N + 1) puncture and one minimal SU(2) puncture. It is important since it will be later used by us to derive four punctured sphere with two maximal USp(2N ) and two minimal SU(2) punctures as specified in (3.7).
To obtain this kind of three punctured sphere theory we start with the trinion T A x,y,z with two maximal SU(N + 1) x,y and one minimal SU(2) z punctures. This trinion was derived in [21] and is shown on the Figure 1. Corresponding superconformal index is given by: where u, v, w and a i parametrize Cartans of the 6d SO(4N + 8) global symmetry. SU(N + 2) gauge symmetry is parametrized by t i 's with the relation Global SU(N + 1) symmetries of the maximal punctures are parametrized by x i and y i satisfying Each puncture has (2N + 6) moment map operators with the charges specified in (2.1). Now in order to obtain desired trinion we should turn one of the SU(N + 1) punctures into USp(2N ) type puncture. For this purpose we use an A N C N tube theory introduced in [24] and shown in Figure 2(a). Chirals of these theory correspond to the moment maps of punctures. For USp(2N ) puncture these are just moment maps specified in (3.1). For SU(N + 1) puncture as specified in (2.1) but with the last moment map flipped so that all of them transform in the same fundamental representation of SU(N + 1): The index of this tube theory is given by: i<j Γ e pq uw −1 4(N +2)x ixj , Figure 8. Spiridonov-Warnaar-Vartanov (SWV) duality. On the picture we also specify superpotentials, UV on the left and dynamically generated on the right sides of dualitites. The solid line starting and ending on the same node denotes multiplet in rank-two antisymmetric representation (AS). Now we take this theory and glue it along the SU(N + 1) puncture to the A N trinion shown in Figure 1. Consistent gluing in this case is the mixture of Φ and S gluings. In particular we Φ glue all moment maps except 1 (u N w 2 ) 2N +4 which is S glued. At the level of the index this operation corresponds to Substituting expressions (C.1) for A N trinion and (C.5) for the A N C N tube we immediately see that the SU(N + 1)x gauge theory appears to be S-confining. Corresponding S-confining theory is depicted in Figure 8. This duality was first introduced by Spiridonov and Warnaar [44] as a mathematical identity for superconformal indices of the corresponding theories: s i . Later Spiridonov and Vartanov [45] discussed physical implications of the relation written above. We refer to this duality as Spiridonov-Warnaar-Vartanov (SWV) duality. More detailed discussion about it can also been found in the Appendix B.5 of our previous paper [19]. Using this duality for the SU(N + 1)x node in (C.6) we can finally write down the index of the desired three-punctured sphere as follows: where we have omitted some of the USp(2N ) singlets thus redefining type of the puncture. Corresponding trinion theory is shown in Figure 9.
Now we can use this A N C N type trinion theory in order to derive four-punctured sphere theory with two USp(2N ) maximal punctures. For this purpose we S glue two copies of this trinion along A N type puncture. At the level of the index this corresponds to the operation specified in (3.7). Substituting trinion index (C.8) into this expression we obtain the following four-punctured sphere index: Now if we look on the SU(N + 1) y gauged node we see that it corresponds to the theory with N + 2 hypermmultiplets, meaning the node is S-confining and can be integrated out using the standard Seiberg duality. Performing this simple operation we land on the theory shown on the Figure 10 with the following superconformal index: where we have also used dictionary (3.4) in order to express everything in terms of C N parametrzation (3.2) since now we have only this type of maximal punctures and hence latter parametrization is more natural. Now we finally close one of the two SU(2) minimal punctures without introducing the defect. As discussed in Section 3 this amounts to computing the residue of (C.10) at the pointz = (pq) − 1 2 w 2 (N +2) 2 N +3ã i . Indeed it can be seen that at this value ofz there is pole coming xx t t Figure 10. Four-punctured sphere theory with two maximal USp(2N ) and two minimal SU(2) punctures. The four-punctured sphere is derived by gluing two trinion theories shown on the Figure 9 and performing some duality transformations specified in the text.
from the contour pinching att wheret N +2 variable is chosen without loss of generality. In fact we should compute such pinchings for eacht i variables and sum the results. But due to the Weyl symmetry of the SU(N + 1) root system all such contributions are the same and summing them results in an overall factor of (N + 1)!. We omit such overall factors since they are irrelevant for the structure of the A∆O we get in the end. Computation of the corresponding residue at the pinching point leads to the following expression for the three-punctured sphere theory index: (C.12) Corresponding quiver diagram is shown in Figure 3. In case N = 1 it can be shown that the theory we derived reduces to the one shown in Figure 22 (b) in our previous paper [19]. There we used it to derive C 1 A∆O which appeared to be van Diejen operator. Studying general N case appears to be too complicated since there is no clear way to further simplify theory (C.12). Instead we concentrate here on C 2 case. When N = 2 AS representation of SU(N + 1) node becomes anti-fundamental representation and the index of the theory reads: with the corresponding quiver of the theory shown in Figure 11 (a). At the next step we perform Seiberg duality on the SU(3) gauge node of the theory above. This node has 9 falvors so after dualising it we obtain SU(6) node with the superconformal index of the resulting theory given by ltj Γ e (pq) and corresponding quiver is shown in Figure 11 (b). We see that SU(4) gauge node of this theory appears to be S-confining due to SWV duality. Hence using SWV index identity (C.7) we arrive to the gauge theory with the following superconformal index: lti Γ e (pq) (c) Figure 11. Chain of duality transformations of three-punctured sphere theory in C 2 case. We start with the theory derived by gluing two three-punctured sphere theories T AC along A-type puncture, integrating out A 2 node and closing SU(2)z minimal puncture. This theory is shown on the Figure  ( with the quiver shown in Figure 11 (c). Now we are ready to close second minimal puncture leaving us with the sphere with only two maximal USp(2N ) punctures left. As discussed in the Section 3 for this purpose we fix z fugacity to be Carefully studying expression (C. 15) we see that at this value of z superconformal index indeed has pole coming from the contour pinching at: where {k 1,± , k 2,± , k 5 , k 6 } ≡ K is partition of the integer K. Of course any permutation of (C.17) would also result in the same contour pinching and hence the pole of the superconformal index. Due to the Weyl group symmetry contributions of all these permutations into our operator is the same and taking them into account just results in an overall factor which we omit. Substituting values (C.17) into three punctured sphere index (C.15) we get the following result for the index of the two-punctured sphere: where C ( K;i) is the following constant: here subscript contains all the information on how the tube was obtained. Inside this subscript 2 means that we obtain two-punctured sphere (the tube) in the end. Index i stands for the choice of the moment map we give vev to according to (3.8) and (3.10). Finally K denotes the power of holomorphic derivative of the moment map we give vev to. Summation in (C.18) goes over all possible partitions of K integer. Due to the presence of the defect this index expression unlike our previous four-and three-punctured sphere does not have nice gauge theory interpretation. Now gluing tube (C.18) to an arbitrary N = 1 theory with one maximal USp(2N ) puncture we can obtain desired A∆O. At the level of the superconformal index the gluing is performed as specified in (3.13): Notice that expression (C.18) for the tube index has zeroes coming from the last two Γfunctions. Due to these zeroes only the poles of the integral in (C.20) coming from other Γ functions will contribute to the final expression. In particular these poles are coming from the Γ e x ±1 k x j q −k j,− and Γ e x ±1 k x −1 j q −k j,+ terms and are located at where σ(i) is permutation of i. In the end we should sum over such permutations as well as over all combinations of ±1 powers in the expression above. Since the tube expression (C.18) is symmetric w.r.t. permutations ofx 1 ↔x 2 andx i →x −1 i every contribution of the pinchings specified above is the same and summation just gives an irrelevant overall factor. So for simplicity we can consider only one of the combinations The condition on m i comes from the fact that both types of specified Γ functions should have poles. Half of these terms should lead to integration contour pinching in (C.20) while other half will cancel zeroes coming from the last two Γ functions in the tube index expression (C.18). Now computing the contribution of pinchings (C.22) into gluing (C.20) and using Γ functions identity (A.9) we can directly obtain the full tower of operators (3.14).

D Proof of the kernel property of A 2 C 2 tube theory
Here we give a proof of the kernel property (4.2) where we act with A 2 operator (2.13) and C 2 operator (3.17) on the index (C.5) of the tube with SU(3) and USp(4) maximal punctures.
Let's start the proof of the kernel property by explicitly calculating action of our A∆Os on the tube index and write two sides of the equation in terms of algebraic expressions. First of all we should fix the parametrization of the Cartans of global 6d symmetry. For this it will be more convenient to use C 2 parametrization given in (3.2). We should also use the dictionary (3.4) between A 2 and C 2 parametrizations in order to rewrite everything in terms of the latter one. Fixing N = 2 in the tube index (C.5) and writing everything in terms of C 2 parametrization we obtain: where as previously x 1,2 are fugacities for the global symmetry of the USp(4) puncture, while y 1,2,3 are fugacities of the SU(3) puncture. Latter ones are as usually constrained by the identity 3 i=1 y i = 1. Now we can start with the r.h.s. of (4.2) describing the action of C 2 operator. Main ingredient we need to find the action of the full operator is the action of the shift ∆ q (x i ) on the tube index. Studying this action on (D.1) we obtain: where the function D l (x, y) is given by: Then on the r.h.s. of (4.2) we obtain: Now we consider l.h.s. of the equation above. Similarly we need to understand action of ∆ lm (y) operator on the A 2 C 2 tube (4.2): where D lm (x, y) = Hence kernel property (4.2) can be reduced to the following algebraic identity: It can be checked that both F (A 2 ;h 10 ;1,0) (x, y) and F (C 2 ;h 10 ;1,0) (x, y) functions defined above are elliptic w.r.t both x i and y i variables with periods 1 and p 4 . In order to prove the identity we now need to check poles and residues in the fundamental domain on two sides of equation 5 . We already know the poles and residues of the constant parts W (A 2 ;h −1 10 ;1,0) (y) and W (C 2 ;h 10 ;1,0) (x) summarized in sections 2 and 3 correspondingly. Now let's study functions F (A 2 ;h 10 ;1,0) (x, y) and F (C 2 ;h 10 ;1,0) (x, y) coming from the action of the shift parts.
We start with F (A 2 ;h 10 ;1,0) (x, y) function given in (D.7). This elliptic function appears to have poles at the following positions: where s = ± and P a ≡    10 P a (D.14) Res 10 P a Looking on these residues we immediately see that on the l.h.s. of (D.8) there are vast cancellation of residues between F (A 2 ;h 10 ;1,0) (x, y) and constant part W (A 2 ;h −1 10 ;1,0) (y). Latter ones can be found from (2.18) upon charges identification given in (3.4). In particular the poles at y a = qy b , y a = q −1 y b , y a = sq where as usually s = ±1. The residues at the poles x i = sq ± 1 2 and x i = sq ± 1 2 p 1 2 are given by It can be easily seen that these residues are cancelled by the corresponding residues (3.19) of the constant part W (C 2 ;h 10 ;1,0) . Hence . Corresponding residues of F (C 2 ;h 10 ;1,0) and hence of the full function F (C 2 ;h 10 ;1,0) + W (C 2 ;h 10 ;1,0) are given by 10 . However as it often happens in our calculations accurate derivation of the corresponding residues results in zeros for all of them 7 . Hence these are not really poles and our function M (x, y) is an entire function for all its variables. Since function is entire elliptic function we conclude that it is just a constant.
In order to define which particular constant it is we just need to put in some values for x and y variables. Here it is worth noticing that M (x, y) function is the sum of two products of θ function. Without loss of generality let's choose the following values of y variables: To calculate these residues it is crucial to use the A2 constraint Now since both functions F (C 2 ;h 10 ;1,0) (x, y) + W (C 2 ;h 10 ;1,0) (x) and F (A 2 ;h 10 ;1,0) (x, y) + W (A 2 ;h −1 10 ;1,0) (y) are elliptic with the coinciding poles and residues in the fundamental domain we conclude that they can differ at most by additive constant. And just as we did it with M ab (x, y) function in order to find this constant we need to choose convenient values of x and y variables. Here it is a bit harder to do then in the case of M ab function. However first thing we can notice is that constant part W (C 2 ;h 10 ;1,0) (x) of C 2 operator (3.17) has zero at the following value of variables: x 1 = (pq) where we have chosen x 1 variable without loss of generality (it can as well be x 2 variable) and the other x 2 variable is kept arbitrary. At the same time constant part W (A 2 ;h −1 10 ;1,0) (y) of A 2 operator given in (2.15) has zero at the following value of variables:

E Commutators of C 2 operators
In this appendix we discuss commutation relations (4.5) of the basic C 2 operators (3.17). Let's start with checking the following commutation relation:  .17) we can immediately conclude that corresponding contribution to the commutator action (E.6) is zero. 7) Terms with ∆ −2 q (x i )I(x). These terms can be directly obtained from Type 2 terms containing ∆ 2 q (x i )I(x) by x i → x −1 i transformation of variables. Hence by the same symmetry argument as in the previous case we can immediately conclude that corresponding contribution to the commutator action (E.6) is zero. 8) Terms with ∆ −1 q (x i )I(x). Contributions to the commutator action (E.6) of this type are also zero which follows from the same argument as previous two terms. In this case it can be obtained by using x i → x −1 i transformation in Type 4 terms. Thus all eight types of terms present in the commutator (E.6) give zero contribution and the full commutator action, and hence the commutator itself, appears to be zero. So far we have shown it only in p and q expansion. But in principle this can be also proven analytically by computing positions of the poles and corresponding residues of all equalities we have checked in expansion. This way of proof is similar to the one we used in Appendix D to prove kernel property (4.2).
The last type of commutator identities we would like to prove is O (C 2 ;ha;0,1) by the exchange p ↔ q. Hence all the arguments above work also for this commutator which as result is indeed equal to zero.
To summarize, in this Appendix we gave arguments in favor of all commutation relations (4.5). For the third commutator we gave full analytic proof while for the first and second we checked commutation relations perturbatively in p and q expansion to sufficiently high order.