Toda conformal blocks, quantum groups, and flat connections

This paper investigates the relations between the Toda conformal field theories, quantum group theory and the quantisation of moduli spaces of flat connections. We use the free field representation of the $\mathcal{W}$-algebras to define natural bases for spaces of conformal blocks of the Toda conformal field theory associated to the Lie algebra ${\mathfrak s}{\mathfrak l}_3$ on the three-punctured sphere with representations of generic type associated to the three punctures. The operator-valued monodromies of degenerate fields can be used to describe the quantisation of the moduli spaces of flat $\mathrm{SL}(3)$-connections. It is shown that the matrix elements of the monodromies can be expressed as Laurent polynomials of more elementary operators which have a simple definition in the free field representation. These operators are identified as quantised counterparts of natural higher rank analogs of the Fenchel-Nielsen coordinates from Teichm\"uller theory. Possible applications to the study of the non-Lagrangian SUSY field theories are briefly outlined.


Introduction
Relations between conformal field theories (CFTs), quantum groups and the quantisation of moduli spaces of flat connections have been investigated extensively in the past in connection with the Chern-Simons theory with compact gauge groups. It can be expected that this subject has a somewhat richer counterpart in the cases where the gauge group becomes non-compact. Interest in the non-compact cases is motivated in particular by recent development in supersymmetric field theory, as will be briefly outlined next.
Field theory X(C, sl 2 ) Liouville theory Moduli of flat connections on C Algebra A X generated by Algebra A V of Verlinde Quantised algebra A M of SUSY loop observables line operators functions on M flat (SL(2), C) Instanton partition Liouville conformal Natural bases for functions blocks modules of A M Table 1: Isomorphic algebras appear in different contexts.

Motivation from N = 2 supersymmetric field theory
An important source of motivation for the study of the non-compact cases are the relations between N = 2 supersymmetric field theories in four dimensions and two-dimensional conformal field theories discovered in [AGT]. This development is deeply related to the study of the relations between d = 4, N = 2 supersymmetric field theories and quantum integrable models initiated in [NS]. An interesting class d = 4, N = 2 supersymmetric field theories, nowadays often referred to as class S, consists of theories X(C, g) labelled by a Riemann surface C and a Lie algebra g of ADE-type [Ga, GMN]. The theory X(C, g) on a four-manifold M 4 is believed to be the effective description for the six-dimensional (2, 0) theory associated to g on M 4 × C when the area of C is small. In the simplest case where g = sl 2 an essential role in this story is played by the algebra A X of supersymmetric Wilson and 't Hooft loops representing an important sub-algebra of the algebra of all observables in X(C, g). The algebra A X turns out to be isomorphic to the algebra A V of Verlinde line operators in Liouville CFT [AGGTV, DGOT, GOP, IOT], which is furthermore isomorphic to the algebra A M of quantised functions on the moduli space of flat SL(2)-connections on Riemann surfaces [TV], as summarised in Table 1.
One can argue [TV] that the isomorphisms between the algebras A X , A V and A M imply the relations between instanton partition functions for X(C, sl 2 ) and conformal blocks of Liouville CFT observed in [AGT] and discussed in many subsequent works.
In the cases where g is a Lie algebra of higher rank one still expects to find a similar picture with Liouville CFT replaced by the conformal Toda CFTs associated to finite-dimensional Lie algebras g. However, even in the next simplest case g = sl 3 one encounters considerable additional difficulties obstructing direct generalisations of the results known for g = sl 2 . On the side of N = 2 SUSY field theory one generically finds interacting field theories not having a weak coupling limit with a useful Lagrangian description. It is not known how to calculate the partition functions in these cases. This difficulty has a counterpart on the side of Toda CFT where the spaces of conformal blocks are known to be infinite dimensional generically, but otherwise poorly understood up to now.
This paper is a step towards the generalisation of the correspondences between the theories X(C, g) and Toda CFT to the cases where no Lagrangian description is known. Motivated by Conformal field theory Quantum group theory Moduli of flat connections Invariants in tensor products Invariants in tensor products of States of W-algebra representations quantum group representations (conformal blocks) Table 2: Comparison of modular functors from conformal field theory, quantum group theory (Reshetikhin-Turaev construction) and quantisation of the moduli spaces of flat connections on Riemann surfaces.
the six-dimensional description of X(C, g), we take the quantised algebra of functions A M as an ansatz for the algebra A X . This is perfectly consistent with previous work on the cases which do have a Lagrangian description [Bu,GLF1,GLF2,LF]. The relations between representations of A M and the algebra of Verlinde line operators in Toda CFT investigated in [CGT] suggest that suitable Toda conformal blocks represent natural candidates for the as yet unknown partition functions of the strongly coupled theories X(C, g). This paper provides groundwork for such a program in the simplest nontrivial case g = sl 3 .
Our work is furthermore motivated by the papers [BMPTY, MP, IMP] where it was proposed that the conformal blocks of Toda CFT can be obtained in a certain limit from the partition functions of topological string theory on certain toric Calabi-Yau manifolds used for the geometric engineering of X(C, g). With the help of the refined topological vertex one may represent the relevant topological string partition functions as infinite series admitting partial resummations. If the geometric engineering limit of these partition functions can be taken, it is expected to represent conformal blocks in Toda CFT. However, at the moment it is not known if this is the case and which basis for the conformal blocks in Toda CFT the topological string partition functions will correspond to. Only when this question is answered one can extract predictions for the Toda three point functions from the results of [BMPTY, MP, IMP]. We'll address the connections to topological string theory in subsequent papers.

Context
One of our goals is to establish a precise relation between conformal blocks in Toda CFT and specific states in the space of states obtained by quantising the moduli spaces of flat connections on Riemann surfaces C. The use of pants decompositions reduces the problem to the basic case C = C 0,3 . Some guidance is provided by the paradigm well-studied in the case of rational conformal field theories summarised in Table 2.
However, the Toda CFTs of interest in this context are not expected to be rational CFTs. We'll therefore need a non-rational analog of the set of relations indicated in Table 2 generalising the situation encountered in the case of Liouville theory, see [TV] and references therein. The experiences made in this case suggest that a very useful role will be played in the non-rational cases by the Verlinde line operators, a natural family of operators V γ on the spaces of con-Moduli of flat connections on C Algebra A V of Verlinde Quantised algebra A M of line operators functions on M flat (SL(2), C) Generators Verlinde line operators V γ Quantised trace functions L γ Module Spaces of conformal blocks Spaces of states formal blocks labelled by closed curves γ on the Riemann surface C. The relevant spaces of conformal blocks can be abstractly characterised as modules for the algebra A V of Verlinde line operators [TV]. This means that these spaces of conformal blocks have a natural scalar product making the Verlinde line operators self-adjoint. Natural representations of A V diagonalise maximal commutative sub-algebras of A V . Such sub-algebras of A V are labelled by pants decompositions. Different pants decompositions give different representations of A V intertwined by unitary operators representing the fusion, braiding and modular transformations.
A key step towards understanding the relation between conformal field theory and the quantum theory of flat connections is therefore to establish the isomorphism between the algebra of Verlinde line operators A V and the algebra of quantised functions on the corresponding moduli spaces of flat connections. In the case of Liouville theory this was established in [TV] using the explicit computations of Verlinde line operators performed in [AGGTV, DGOT]. As discussed in [TV] one finds that the Verlinde line operators generate a representation of the quantised algebra A M of functions on M flat (PSL(2, R)). The generator of A M corresponding to the Verlinde line operator V γ is the quantised counterpart L γ of a trace function on M flat (PSL(2, R)).
The situation is schematically summarised in Table 3.
An important feature of the representations associated to pants decompositions is the fact that the Verlinde line operators get expressed in terms of a set of more basic operators representing quantised counterparts of certain coordinates on M flat (PSL(2, R)) which are close relatives of the Fenchel-Nielsen coordinates for the Teichmüller spaces [TV]. The relevance of such coordinates had previously been emphasised in a related context in [NRS]. The distinguished role of these coordinates in our context is explained by the fact that they have two important features: (i) They represent the algebraic structure of the moduli spaces of flat connections in the simplest possible way, and (ii) they are compatible with pants decompositions.

Contents of this paper
The free field representation of the W-algebras allows us to define convenient bases for spaces of conformal blocks. Further developing the known connections between free field representations of conformal field theories and quantum group theory we will relate the spaces of conformal blocks to spaces of intertwining maps between tensor products of quantum group repre-sentations and irreducible representations. This reduces the computation of monodromies to a problem in quantum group theory. The monodromy transformations get represented in terms of operator-valued matrices with matrix elements represented as finite-difference operators acting on the multiplicity labels for bases of conformal blocks. Our explicit computations allow us to identify basic building blocks for the operator-valued monodromies which can be identified with (quantised) coordinates for moduli spaces of flat connections generalising the Fenchel-Nielsen coordinates mentioned above.
Coordinates for M flat (PSL(3, R)) satisfying the conditions (i), (ii) characterising coordinates of Fenchel-Nielsen type have been previously been introduced in [Go, Ki]. A systematic geometric approach to the definition of higher rank analogs of the Fenchel-Nielsen coordinates was proposed in [HN, HK]. It would be very interesting to understand the relation between the coordinates defined in these references and the ones presented in this paper.
Some of the braiding matrices in Toda CFT have been calculated in [GLF1,GLF2,LF]. However, in all these cases the two representations associated to the punctures involved in the braiding operations are degenerate or semi-degenerate, corresponding to theories X(C, g) which have a Lagrangian description.
Section 2 introduces relevant background from conformal field theory and explains how to use conformal blocks with degenerate field insertions to define operator-valued analogs of monodromy matrices. A construction of conformal blocks for the W 3 algebra using the free field representation is introduced in Section 3. It is explained how this construction relates bases for the spaces of conformal blocks on C 0,3 to bases in the space of Clebsch-Gordan maps for the quantum group U q (sl 3 ). The connection to quantum group theory provided by the free field representation is used in Section 4 to calculate basic cases of the monodromies of degenerate fields around generic chiral vertex operators in sl 3 Toda CFT. The results of this computation are then used in Section 5 in order to relate the basic building blocks of the operator valued monodromy matrices to coordinates of Fenchel-Nielsen type on the moduli spaces of flat SL(3)-connections. There exist two different limits in Toda CFT corresponding to a classical limit for A V , b → 0 and b → i, which both yield a deformation parameter q 2 = e −2πib 2 equal to 1. We point out relations to the Yang's function for the Hitchin system, and to the Riemann-Hilbert problem for holomorphic SL(3)-connections on C 0,3 arising in the two limits, respectively. In the concluding Section 6 we discuss in particular the extension of our results to continuous families of conformal blocks as expected to be relevant for a full understanding of Toda CFT.

Conformal blocks and quantum monodromies
The A N −1 Toda CFTs are defined by the Lagrangian where ϕ is a N − 1 component scalar field, the vectors e i are the simple roots of the Lie algebra sl N , and (., .) denotes the standard scalar product in R N −1 .

Chiral algebra W N and its free field construction
The chiral symmetry algebra of A N −1 Toda CFT is a W N -algebra generated by holomorphic currents W d (z), d = 2, . . . , N . A quick way to introduce the algebras W N starts from a collection of N − 1 free chiral bosons, with modes satisfying the commutation relations can then be defined through a deformed version of the Miura transformations, where :: denotes Wick ordering and h i = ω 1 − i−1 j=1 e j are the weights of the fundamental representation of sl N with highest weight ω 1 . The chiral algebra W N defined in this way will be the chiral symmetry of Toda CFT if the parameter Q in (2.5) is related to the parameter b in the Lagrangian (2.1) as The currents T (z) = W 2 (z) generate a Virasoro sub-algebra within W N with central charge c = (N − 1) 1 + N (N + 1)Q 2 . (2.6) Highest weight representations of the algebras W N are generated from highest weight vectors e α labelled by an N − 1 component vector α satisfying where w d (α) are polynomials of α of degree d, with being the conformal dimension. We are here using the notation ρ W for the Weyl vector of sl N , which is equal to half the sum of all positive roots.
The space of states of the Toda CFTs will then be of the form H = S dα V α ⊗V α , whereV α is a representation for the algebra W N generated by the anti-holomorphic currentsW d (z). The set of labels S for the representations appearing in H is expected to be the set of vectors α of the form α = Q + iP , P ∈ R 2 , Q = ρ W (b + b −1 ), ρ W being the Weyl vector of sl N .

Definition of conformal blocks
One of the basic problems in Toda CFT originates from the fact that the spaces of conformal blocks corresponding to the three-punctured spheres are infinite-dimensional, as will briefly be reviewed below. This implies that the Toda CFT three-point function can be represented in the following form (2.9) In (2.9) we are using the notation V α (v ⊗ w; z,z) for the vertex operator associated to a state v ⊗ w ∈ V α ⊗ V α by the state-operator correspondence. It can be represented as a normal ordered product of the primary field V α (z,z) = V α (e α ⊗ e α ; z,z) associated to the highest weight state e α ⊗ e α with differential polynomials in the fields representing the W 3 -algebra. The superscript ρ refers to the triple of labels ρ = (α 3 , α 2 , α 1 ) for the representations of the Walgebra associated to the vertex operators appearing in the correlation function (2.9). f ρ k andf ρ k are bases for the relevant (sub-)spaces of the spaces of conformal blocks on the three-punctured spheres C 0,3 = P 1 \ {z 1 , z 2 , z 3 } having W-algebra representations with labels α i assigned to the punctures z i for i = 1, 2, 3, respectively.
We will now review the basic definitions and results on the conformal blocks of W-algebras that will be used in the following. The presentation will be rather brief, the formulation being to a large extend analogous to the one presented in [T17] for the case of the Virasoro algebra.
The conformal blocks on C 0,n = P 1 \ {z 1 , . . . z n }, are defined as the solutions to the Ward identities expressing the W-algebra symmetry of A N −1 Toda CFT under the W N -algebra. For any collection of meromorphic (1 − d)-differentials ξ d = ξ d (y)(dy) 1−d , d = 2, . . . , N , with poles only at y = z r , r = 1, . . . , n, we may define operators W d [ξ d ] on the tensor product d,k being the coefficients of the Laurent expansions of ξ d around z r , r = 1, . . . , n, and W d,k (z r ) representing W d,k on the r-th tensor factor of R associated to the puncture z r , (2.12) Using these notations we will define conformal blocks 1 as linear maps f : R → C satisfying for all ξ d as introduced above. The set of linear equations (2.13) defines a subspace CB(C, R) in the dual R of the vector space R called the space of conformal blocks.
We will often find it convenient to assume that z n = ∞ which does not restrict the generality in any serious way. By using a part of the defining invariance conditions (2.13) one may reduce the computation of the linear functionals f (v), v ∈ R, to the special values 2 (2.14) There is a one-to-one correspondence between functionals f : R → C satisfying (2.13), and functionals g : V αn → C satisfying a reduced system of invariance conditions following from (2.13) and (2.14). We will therefore also call such functionals g conformal blocks.
In this paper will mostly discuss the case N = 3 corresponding to sl 3 Toda theory in the following, where we will use the simplified notations We expect to be able to reduce the case of an n-punctured sphere to the one for the threepunctured sphere by means of the gluing construction. For this case it is straightforward to show that the defining invariance property (2.13) allows one to compute the values of f (v) on . This means that any conformal block f on C 0,3 is uniquely defined by the sequence of complex numbers F f = (F f,l ) l∈Z ≥0 . As opposed to the case of the Virasoro algebra we therefore find an infinite-dimensional space of conformal blocks in the case of W 3 .
The representation of the W 3 -algebra attached to the puncture z r generated by the operators W d,n (z r ) induces an action on the spaces of conformal blocks CB(C, R) via (2.16) The conformal block W d,n (z r )f obtained by acting on f with W d,n (z r ) is characterised by a sequence (F W d,n (zr)f,l ) l∈Z ≥0 which may be computed by specialising the definition (2.16) to vectors of the form v = e α 1 ⊗ e α 2 ⊗ (W −1 ) l e α 3 and using (2.13) to express the right hand side of (2.16) as a linear combination of the F f,l , l ∈ Z ≥0 .
Before one addresses the problem to compute the Toda three-point functions C k,k (α 1 , α 2 , α 3 ) using an expansion of the form (2.9) one needs to find useful bases {f µ k ; k ∈ I} for the space of conformal blocks appearing in (2.9), with I being a suitable index set. According to the discussion above one might be tempted to consider the basis {f l ; l ∈ Z ≥0 } for CB(C, R) defined such that F f l ,n = δ l,n . However, such a basis does not appear to be the most useful basis for the purpose to construct physical correlation functions in the form (2.9). Bases which appear to be more useful in this regard will be defined using the free field representation for the W-algebras.
It will also be useful to keep in mind the one-to-one correspondence between conformal blocks for three-punctured spheres C 0,3 = P 1 \ {0, z 2 , ∞} and chiral vertex operators V ρ (v 2 , z 2 ) : (2.17) using an invariant bilinear form . , . Vα 3 on V α 3 . One may, more generally, relate conformal blocks f V to families of multi-local chiral vertex operators V (w; z) : V α 1 → V αn , w ∈ V α 2 ⊗ · · · ⊗ V α n−1 , z = (z 1 , . . . , z n−1 ), defined such that Such multi-local vertex operators can be constructed as compositions of vertex operators V ρ (v, z). We may therefore construct large families of conformal blocks by constructing the chiral vertex operators V ρ (v, z), or equivalently the conformal blocks associated to C 0,3 . This will be our main goal here.

Degenerate representations
A special role is played by the vertex operators associated to degenerate representations of the W 3 -algebra. We will only need the examples where the representation label α is equal to either −bω 1 or −bω N −1 , with ω 1 and ω N −1 being the weights associated to the fundamental and antifundamental representation of sl N , respectively. In these cases one will find relations called null vectors among the vectors generated by the algebra W N from the highest weight vector e α . In the case N = 3 one will find the following three null vectors in the representation V −bω 1 : n 2 : 0 = (L 2 −1 + κ 2,1 L −2 + κ 2,2 W −2 )e −bω 1 , (2.20) Explicit formulae for the coefficients can be found e.g. in [FL]. The relations in the representation V −bω N −1 take a similar form, obtained by replacing κ i,i by −κ i,i , i = 1, 2, 3.
The null vectors in degenerate representations have well-known consequences for the vertex operators associated to such representations. If α 2 = −bω 1 , for example, it follows from the relations above that V ρ (v 2 , z 2 ) : V α 1 → V α 3 can only be non-vanishing if α 3 = α 1 − bh ı where h ı are the weights of the fundamental representation of sl 3 for ı = 1, 2, 3. For α 2 = −bω N −1 one will similarly find α 3 = α 1 + bh ı for ı = 1, 2, 3. We may therefore abbreviate the notations for the corresponding vertex operators to D ı (y) andD ı (y), respectively. They map It is furthermore well-known that these relations imply differential equations satisfied by the conformal blocks if some of the tensor factors of R = n r=1 V αr are degenerate in the sense above. To be specific we will in the following consider conformal blocks on a sphere with n + 2 punctures z 1 = 0, z 2 , . . . , z n−1 , y, y 0 , ∞, with representation V −bω N −1 associated to y 0 and representation V −bω 1 associated to y. As in (2.18) we may then consider the conformal blocks F V (y, y 0 ; z)ī ı = e n ,Dī(y 0 )D ı (y) V (e α 2 ⊗ · · · ⊗ e α n−1 ; z) e α 1 Vα n . (2.23) The differential equations following from the null vectors include, in particular, (2.25) Note that the null vectors n 1 allow us to express W −1 (y 0 ) in terms of L −1 (y 0 ) ∼ ∂ ∂y 0 . We may furthermore use the null vector n 2 together with the Ward identities to express W −2 (y 0 ) in terms of a differential operator W −2 (y 0 ) constructed out of derivatives ∂ ∂y 0 and L −1 (z r ). However, even in the case n = 3 one can not use the invariance conditions defining the conformal blocks to express all W −i (z r )F V (y; y 0 ; z)ī ı , i = 1, 2 in terms of derivatives of F V (y; y 0 ; z)ī ı , in general.
Not having a closed system of differential equations represents an obstacle for the use of the differential equations (2.24). Our goal is to define the analytic continuation of F V (y; y 0 ; z)ī ı with respect to y and compute the resulting monodromies. Without a more explicit description of the action of the algebra W 3 on spaces of conformal blocks one can not use the differential equations discussed above to reach this goal. This is the main difference to the case where N = 2. Other methods are needed to reach our goal, as will be introduced next.

Free field construction of conformal blocks
We propose that useful bases for the spaces of conformal blocks can be defined using the free field representation of the algebra W 3 . This section gives an outline of the necessary constructions which is sufficient for the formulation of our main results. A more detailed description will be given in Section 3 below.
Key elements are the normal ordered exponentials V α (z) = e (α,ϕ(z)) . (2.26) One of the most important properties for the following will be the exchange relations where the left hand side is understood as the analytic continuation from a region |z | > |z|, and z moves in a counter-clockwise direction around z. The special fields S i (z) = V be i (z) will be referred to as the screening currents. From the screening currents S i (z) one can construct screening charges Q i (γ) as γ being a suitable contour. Apart from the screening charges associated to the simple roots e i , i = 1, 2, it will frequently be useful to consider the composite screening charge Q 12 defined as with γ y being a suitable contour encircling the point y which will be specified more precisely below.
As will be discussed in more detail in Section 3 we may use these ingredients to construct conformal blocks on C 0,3 by expressions of the following form The precise definition of the expression on the right hand side of (2.30) requires further explanations. The integrand is the multi-valued analytic function obtained by standard normal ordering of the exponential fields appearing in the matrix element on the right hand side of (2.30). Both the contours of integration and the precise choice of branch defining the integrand still need to be specified. All this will be discussed in more detail in Section 3 below.
The vertex operators appearing in (2.30) map from a module V α 1 of the W 3 algebra to the module F α 3 , where α 3 is related to α 2 , α 1 and n = (m, n, p) by the conservation of the free field momentum, Fixing the triple ρ = (α 3 , α 2 , α 1 ) of representation labels only fixes the two combinations s 1 = m + n and s 2 = n + p of the three parameters n = (m, n, p). It is therefore natural to parameterise n in terms of s 1 , s 2 and an additional multiplicity label k such that For fixed external parameters ρ = (α 3 , α 2 , α 1 ) one may therefore use (2.30) to construct a one-parameter family of conformal blocks F ρ k .

Quantum monodromies
A useful probe for the conformal blocks F ρ k will be the 3 × 3 matrix of conformal blocks on C 0,3+2 having matrix elements denoted by F ρ k (y; y 0 )ī ı , ı,ī = 1, 2, 3, defined as (2.33) We will regard y 0 as a base-point, typically assumed to be near ∞ on C 0,3 and sometimes use the abbreviated notation F ρ k (y)ī ı ≡ F ρ k (y; y 0 )ī ı . We are going to demonstrate in Section 4 below that F ρ k (y)ī ı has an analytic continuation with respect to y generating a representation of the fundamental group π 1 (C 0,3 ) of the form (2.34) We will see that only finitely many terms are nonzero in the summation over l. It will be useful to regard F ρ k (y)ī ı as the components of vectors F ρ (y)ī ı in a vector space with respect to a basis labelled by k, and accordingly the coefficients (M γ ) ,k ı,l as matrix elements of operators (M γ ) ı  , ı,  = 1, 2, 3, allowing us to write (2.34) in the form We may anticipate at this point that the matrices (M γ ) ı  will turn out to be "quantum" analogs of the monodromy matrices of a flat SL(3)-connection on C 0,3 .
The derivation of the relations (2.35) described below will allow us to show that the elements (M γ ) ı  of the quantum monodromy matrices can be expressed as functions generating a Weyl algebra W characterised by the commutation relation uv = q −2 vu, using the notation q = e −iπb 2 . It will turn out that the M γ (u, v) ı  are Laurent polynomials of low orders in u, v. This property will be essential for the interpretation of u and v as quantised analogs of Fenchel-Nielsen type coordinates on moduli spaces of flat connections to be discussed in Section 5.

Free field construction of chiral vertex operators
We will now describe how to define the chiral vertex operators (CVOs) precisely, and how to compute the monodromies resulting from the analytic continuation of the location of a degenerate chiral vertex operator around the support of a generic CVO. One needs to specify a suitable set of contours for the integration of the screening currents in order to complete the construction of the chiral vertex operators. Monodromies may then be calculated by deforming the contours appropriately. The main goal of the following two sections will be to derive the claims on the structure of the quantum monodromy matrices formulated at the end of the previous section. The consequences of these results concerning the relations to the moduli spaces of flat connections will be explained in Section 5 below.
The approach described below is inspired by the works [GS, FW] where the connection between the free field representation of the Virasoro algebra and the quantum group U q (sl 2 ) was investigated. Relations between free field representations of other conformal field theories and quantum groups of higher ranks have been studied in [RRR, BMP, BFS]. However, the aims and scope of these works are different from ours, making it difficult to extract the results relevant for us from the references above. We therefore give a self-contained treatment below.

Basic definitions
For both tasks it will be useful to represent conformal blocks as linear combinations of a family of auxiliary matrix elements of the following form over the contoursΓ = Γ nm m × . . . × Γ n 1 1 depicted in Figure 1. The crosses in Figure 1 indicate normalisation points on the contourΓ where the multi-valued integrand in (3.2) is real if z r ∈ R + , r = 2, . . . , m, and 0 < z 1 < · · · < z m . We are using the notations n l = (n l,1 , n l,12 , n l,2 ) for l = 1, . . . , m, and and assume that The contours used to define the composite screening currents S 12 (y) in (2.29) should be chosen such that the only singular point of the integrand in (3.2) encircled by γ y is the point y.
The functions W αm,...,α 1 nm,...,n 1 do not by themselves represent conformal blocks, as is easily seen as follows. The commutators between screening currents S i (z) and the generators of the W 3algebra are total derivatives with respect to the variable z. When the screening currents S i (z) ) S 12 ( y r ) Figure 1: Nested non-intersecting and non-self-intersecting multi-contours for integration over screening currents associated with different vertex operators. The multi-contours around a puncture and which are associated with a particular screening charge are depicted collectively by one loop.
are integrated over open contours as depicted in Figure 1, one will therefore find boundary terms in the commutators between screening charges and the generators of the W 3 -algebra supported at the base-point z 0 . The boundary terms will spoil the validity of the Ward identities for the W 3 -algebra, in general.
However, we will see that the boundary terms all cancel for specific choices of the coefficients C in (3.1). Representing compositions of chiral vertex operators in the form (3.1) is useful since both the problem to find the coefficients C in (3.1) and the calculation of monodromies will be shown to reduce to purely algebraic problems in quantum group theory. In order to explain this relation we shall begin by collecting the definitions and notations concerning the quantum group U q (sl 3 ) relevant for our goals.

Quantum group background
The quantum group U = U q (sl N ) is a Hopf algebra, generated as an associative algebra over C(q) by e i , f i and k ± i for i = 1, . . . , N − 1 which satisfy the relations with |i − j| = 1. Let us now restrict to the case N = 3. Defining it follows from the Serre relations that e 1 e 12 = q −1 e 12 e 1 , e 2 e 12 = qe 12 e 2 , (3.8) (3.9) Using these relations it is possible to show that arbitrary monomials formed out of the basic generators f i , i = 1, 2 can be represented as linear combinations of the ordered monomials f n 1 1 f n 12 12 f n 2 2 with (n 1 , n 12 , n 2 ) ∈ N 3 .

Representations
The irreducible highest weight representation of U q (sl 3 ) with highest weight λ will be denoted by R λ . The highest weight vector v λ ∈ R λ is annihilated by the generators (3.10) The representations R λ are irreducible for generic weights λ. We call a weight λ generic if it is not the highest weight of a finite-dimensional highest weight representation of sl 3 . We will use a basis for the U q (sl 3 ) Verma module R λ which consists of the following vectors e λ n = f n 1 1 f n 12 f n 2 2 v λ , n = (n 1 , n, n 2 ) .
( 3.11) The vector e λ n has weight ν = ν(λ, n) = λ − (n + n 1 )e 1 − (n + n 2 )e 2 . (3.12) It will sometimes be useful to parameterise the labels n = (n 1 , n, n 2 ) in terms of another pair of data n = [δ, n] consisting of the shift of weights δ = −s 1 e 1 − s 2 e 2 together with a multiplicity label n such that s 1 = n 1 + n and s 2 = n 2 + n. For a fixed weight ν = λ + δ parameterised by positive integers s 1 and s 2 , one can only have finitely many values of n between 0 and d ν − 1 = min(s 1 , s 2 ). The multiplicity space of a weight ν has dimension d ν and will be denoted by M ν λ . The action of the generators x of U q (sl 3 ) is represented in the basis (3.11) by the matrices R λ,n (x),n defined by π λ (x)e λ n = n R λ,n (x),n e λ n .
( 3.13) Tables for the matrix elements of R λ,n (x),n for x = e i , f i can be found in appendix A. Here we will only note the following important feature concerning the dependence on the multiplicity labels k. The action of π λ (x) on the basis vector e λ n can for x = e i , f i , i = 1, 2, be represented in the following form (3.14) δ,k+1] , and the coefficient functions E i , E i , F i and F i are Laurent polynomials in the variable u k . This simple feature will play a crucial role for us in the following.
Apart from the representations associated to generic weights λ we will be interested in particular in the fundamental representation generated from the highest weight vector v ω 1 and having a basis {v ω 1 , f 1 v ω 1 , f 12 v ω 1 }.

Tensor products of representations
The coproduct ∆ of U q (sl 3 ) is defined by The coproduct is an algebra homomorphism, so equations (3.15) imply, for example, The Clebsch-Gordan (CG) maps describe the embedding of irreducible representations into tensor products of representations. In the case of the highest weight representations introduced above, one may represent the CG maps C λ 3 Here and below we use the notation A set of equations for the Clebsch-Gordan coefficients (CGC) λ 3 n 3 λ 2 λ 1 n 2 n 1 follows from the intertwining property for any element x ∈ U q (sl 3 ). When x = e i ∈ {e 1 , e 2 } and n 3 = 0, this equation implies that for every fixed pair (n 1 , n 2 ) Equations (3.19) can be regarded as a set of recursion relations. For generic values of λ i , i = 1, 2, 3 one can show that these equations determine λ 3 0 λ 2 λ 1 n 2 n 1 uniquely in terms of λ 3 0 λ 2 λ 1 n 2 0 . Indeed, it is not hard to see (see Appendix B for further details) that one can use equations (3.19) to derive a recursion relation expressing λ 3 is not further constrained by (3.19). A basis for the space of CG maps C λ 3 λ 1 λ 2 is therefore obtained from the CGC λ 3 n 3 λ 2 λ 1 n 2 n 1 k characterised by the property that (3.20) As k 2 can only take a finite set of values determined by δ 2 we see that the coefficients λ 3 n 3 λ 2 λ 1 n 2 n 1 k define a basis for the space of CG maps. For generic weights λ 1 , λ 2 , λ 3 there is an isomorphism between the space of CG maps and the subspace with weight λ 3 − λ 1 within R λ 2 , or equivalently the subspace with weight λ 3 − λ 2 within R λ 1 .
It will be important to keep in mind that the λ 3 0 λ 2 λ 1 n 2 n 1 k defined using (3.20) depend on the multiplicity labels only through the matrix elements R λ,n (e i ),n . The features concerning the dependence of R λ,n (e i ),n on the multiplicity labels observed in (3.14) will be inherited. Note furthermore that the recursive procedure determining the CGC λ 3 0 λ 2 λ 1 n 2 n 1 k proceeds by induction in s = s 1 + s 2 if n 1 = [−s 1 e 1 − s 2 e 2 , k 1 ]. In each step one generically increases the range of values for the multiplicity label k 2 in n 2 = [δ 2 , k 2 ] with non-vanishing λ 3 0 λ 2 λ 1 n 2 n 1 k by one unit. In the case n 1 = [−s 1 e 1 − s 2 e 2 , k 1 ] with s 1 + s 2 ≤ 2, for example, we find non-vanishing λ 3 0 λ 2 λ 1 n 2 n 1 k only for k − 2 ≤ k 2 ≤ k. Further details are given in Appendix B.

R-matrix
The universal R-matrix is an invertible element R ∈ U q (sl N ) ⊗ U q (sl N ) with the property where P ∈ U q (sl N ) ⊗ U q (sl N ) is the permutation that acts as P(x ⊗ y) = y ⊗ x. It satisfies the Yang-Baxter equation R 12 R 13 R 23 = R 23 R 13 R 12 and the action of the coproduct on it is It can be shown that the conditions (3.23) have a unique solution of the form R = q tR with t = ij (κ −1 ij )h i ⊗ h j , with κ ij being the Cartan matrix, if k i = q h i , and R = 1 + (q − q −1 )(e 1 ⊗ f 1 + e 2 ⊗ f 2 − q −1 e 12 ⊗ f 12 ) + (q − q −1 ) 2 e 2 e 1 ⊗ f 2 f 1 + . . . (3.25) The omitted terms in the expression (3.25) are of higher order in e i , i = 1, 2.
The universal R-matrix is related to the braid matrices B λ 1 λ 2 : is the evaluation of the universal R-matrix in the tensor product of representations R λ 1 ⊗ R λ 2 and P is the permutation of tensor factors. Our notation for the matrix elements of B λ 1 λ 2 is defined such that B λ 1 λ 2 e λ 1 n 1 ⊗ e λ 2 n 2 = m 1 ,m 2

Braiding of screened vertex operators
The spaces of functions generated by the functions W αm,...,α 1 nm,...,n 1 (z m , . . . , z 1 ) introduced in (3.2) carry a natural action of the braid group which turns out to be related to the representation of the braid group on tensor products of quantum group representations. We will now briefly explain this first crucial link between quantum group theory and free field representation to be used in this paper.
Let us consider W α 2 α 1 n 2 n 1 (z 2 , z 1 ) defined using the contours from Figure 1. The braiding operations are defined using the analytic continuation of this function with respect to z 2 and z 1 along a path exchanging their positions. This analytic continuation can be defined by integrating along suitably deformed contours as indicated in Figure 2. By means of contour deformations one can represent the result of this operation as a linear combination of integrals over the contours introduced in Figure 1.
This can be done by the following recursive procedure. In a first sequence of steps, performed in an induction over n 2 = (n 2,1 , n 2,21 , n 2,2 ) one may deform the contours on the right of Figure  2 into linear combinations of the contours depicted in Figure 3 .
The induction step is indicated in Figure 4.
In a second sequence of steps one may then deform the contours surrounding both z 1 and z 2 in Figure 3 into a linear combination of contours surrounding only z 1 or z 2 at a time. The result of this procedure can be represented as linear combination of the original contours introduced in z 0 z 2 n 21 z 1 n 1 Figure 3: A set of auxiliary contours for the braiding calculation.

Figure 1, leading to a relation of the form
We will explain in more detail in Appendix C below how such computations can be done and calculate the matrix elements of b α 1 α 2 for the important special case where α 2 = −bω 1 with α 1 generic. The result can be compared to the action of the braid group on tensor products of quantum group representations defined in (3.27). It turns out that the matrix b α 1 α 2 appearing in (3.28) coincides with the matrix B λ 1 λ 2 in (3.27) provided that λ l = −α l /b for l = 1, 2. This result can easily be generalised to the cases where α 2 = −bλ 2 with λ 2 being the weight of any finite-dimensional representation of U q (sl 3 ) by noting that any finite-dimensional representation of U q (sl 3 ) appears in the iterated tensor products of fundamental representations, and that the co-product of the quantum group has a simple representation in terms of the free field representation, as indicated in Figure 5 below.
2n 2,1 +n 2,12 −n 2,2 Figure 5: Decomposition of a simple contour associated with the root e 1 around two punctures into contours enclosing only one of the punctures. This decomposition represents the co-product As representations of the braid group one may therefore identify the vector spaces generated by the functions W αm,...,α 1 nm,...,n 1 (z m , . . . , z 1 ) with the tensor product of quantum group representations m l=1 R λ l assuming that the weights are related as λ l = −α l /b, l = 1, . . . , m.

Relation to quantum group theory
In order to see the relation between the problem to find linear combinations of the form (3.1) representing conformal blocks and the Clebsch-Gordan problem one needs to notice that the boundary terms appearing in the commutators of the screening charges with the generators of the algebra W 3 turn out to be related to the action of the generators e i on tensor products of quantum group representations defined by means of the co-product. In order to formulate a more precise statement, let us introduce some useful notations. Let V αm,...,α 1 nm,...,n 1 (z m , . . . , be the operator appearing in the matrix elements (3.2) defining the functions W αm,...,α 1 nm,...,n 1 . The correspondence with vectors in tensor products of quantum group representations observed above can be schematically represented as allowing us to introduce the notation ∆(x)V αm,...,α 1 nm,...,n 1 (z m , . . . , z 1 ) (3.32) for the vertex operator associated to ∆(x)e λ 1 n 1 ⊗ . . . ⊗ e λm nm for x ∈ U q (sl 3 ) via the correspondence (3.31). In the computation of the commutator of the vertex operators V αm,...,α 1 nm,...,n 1 with the generators of the algebra W 3 one may distinguish terms from the commutators of the screening currents from the terms coming from commutators with other exponential fields. The contribution of the former is then found to be a linear combination of terms which are proportional to (∆(e i )V αm,...,α 1 nm,...,n 1 )(z m , . . . , z 1 ). Considering, for example, the commutator with L −1 we find L −1 , V αm,...,α 1 nm,...,n 1 (z m , . . . , z 1 ) = m l=1 ∂ z l V αm,...,α 1 nm,...,n 1 (z m , . . . , z 1 ) + (boundary terms), (3.33) with boundary terms proportional to (∆(x)V αm,...,α 1 nm,...,n 1 )(z m , . . . , z 1 ) for (3.34) The derivation of these statements is outlined in Appendix D. We conclude that the boundary terms in the commutators with generators of the algebra W 3 will vanish if the coefficients C in (3.1) are taken to be the expansion coefficients of the images of highest weight vectors under the CG maps.

Independence of the choice of base-point
The linear combinations (3.1) representing conformal blocks do not depend on the position of the base point z 0 indicated in Figure 1. Indeed, the derivative of W αm,...,α 1 nm,...,n 1 with respect to z 0 yields boundary terms which cancel in the linear combinations (3.1) by the same arguments as used above.
The independence of z 0 can be used for our advantage in two ways. One may note, on the one hand, that the expressions for conformal blocks resulting from (3.1) can be replaced by expressions where different base points z 0 are associated to each individual CG coefficient appearing in expressions like (3.29). In (3.29), for example, one could have a base-point z 0 for the integrations with screening numbers n 1 and n 2 , and another base-point z 0 for those with screening numbers n 12 and n 3 , as depicted in Figure 6. In this form it is manifest that the conformal blocks defined via (3.29) factorise into three-point conformal blocks, as they should. This representation furthermore exhibits clearly the relation between the intermediate representation V α 12 of the algebra W 3 appearing in the conformal blocks on the left side of (3.29) and the intermediate representation R λ 12 appearing in the Clebsch-Gordan coefficients on the right side of (3.29), with representation labels related by α 12 = −bλ 12 .
It can also be useful to consider the limit where z 0 → z 2 . In this case one may use contours supported near the circles around the origin with radius |z 2 |, allowing us represent the CVOs as composite operators of the form Q n 1 1 (γ)Q n 12 12 (γ)Q n 2 2 (γ)V α 2 (z 2 ), with screening charges Q i (z) = γ S i (y) being defined by integration over γ = {y ∈ C; |y| = |z 2 |} for i = 1, 12, 2.

Computation of monodromies
We will now explain how the connection between free field representation and quantum group theory helps us to compute the representation of monodromies on spaces of conformal blocks. Recall that we are interested in calculating the monodromies of the conformal blocks F ρ k (y; y 0 )ī ı = e α 3 ,Dī(y 0 )D ı (y) V ρ k (z) e α 1 , (4.1) regarded as function of y. The task can be simplified slightly by sending the base-point y 0 to infinity, reducing the problem to the computation of the monodronies of the four-point conformal blocks where v ∞ ∈ V α 4 with α 4 being related to the variables α 3 andī in (4.1) as α 4 = α 3 − bhī.
The fundamental group of C 0,3 = P 1 \{0, z, ∞} is generated by the loops γ 0 , γ z , γ ∞ around 0, z and ∞. Since the loop γ 0 • γ z • γ ∞ is contractible, it suffices to compute the monodromies along two out of the three loops. The monodromy of the conformal blocks G ρ k (y) ı introduced in (4.2) around ∞ is a simple diagonal matrix. It will therefore suffice to compute the monodromy of G ρ k (y) ı along a loop surrounding z 2 = z. The free field representation introduced above relates this monodromy to where D is the linear operator defined in (3.29), and B λ 3 λ 2 : R λ 3 ⊗ R λ 2 → R λ 2 ⊗ R λ 3 is the braiding operation.
As explained previously, we will mainly be interested in the case where the representation R λ 3 assigned to z 3 = y is the fundamental representation, λ 3 = ω 1 . The simplifications resulting from λ 3 = ω 1 will be discussed next.

Braiding with fundamental representation
The form of the braid relations also simplifies when λ 2 = ω 1 or λ 3 = ω 1 . Let us first consider the case λ 2 = ω 1 . The braid matrix preserves weight subspaces in the tensor product. Recalling that the only allowed weights in the fundamental representation are h ı , ı = 1, 2, 3 one notes that the subspace in R ω 1 ⊗ R λ having fixed total weight ν is spanned by vectors of the form e ω 1 δı ⊗ e λ n with weight of e λ n being ν − h ı for ı = 1, 2, 3. For fixed ı one gets a subspace of R ω 1 ⊗ R λ isomorphic to the multiplicity space M ν−hı λ . The projection of the braiding operator B ω 1 λ : R ω 1 ⊗ R λ → R λ ⊗ R ω 1 onto a subspace of total weight ν can be represented by an upper triangular matrix b + ≡ b + (λ, ν) of operators b + ı mapping from M ν−hı λ to M ν−h λ for ı,  = 1, 2, 3.

The operators b +
ı can be conveniently represented as follows. The basis (3.11) allows us to introduce a basis for the multiplicity spaces Out of these operators one may then construct the operatorsb + ı aŝ b + 11 = q (h 1 ,ν) ,b + 22 = q (h 2 ,ν+e 1 ) ,b + 33 = q (h 3 ,ν+e 12 ) (4.10) It can easily be checked that the image of the restriction ofb + ı to M ν−hı λ is contained in the subspace of M isomorphic to M ν−h λ . Using these definitions one may represent our results for the operators b + ı as the projections ofb + ı to operators from M ν−hı In the second case λ 3 = ω 1 one may similarly represent the projection of B λω 1 onto the subspace of weight ν by a matrix b − ≡ b − (ν) having matrix elements representable through A feature of particular importance for us is the fact that the braid matrix can be represented as a matrix having elements which are finite difference operators of second order in the multiplicity label. This is a direct consequence of (3.14) given the relation between braiding and quantum group R-matrix, as is easily seen using (3.25).

Form of the monodromy matrix
With both of the braid matrices b + and b − explicitly derived, we now find that the matrix M 2 representing the monodromy around z 2 takes the form M 2 = D −1 b − b + D with D being a lower triangular matrix with matrix elements acting as pure multiplication operators on the multiplicity spaces. It may be useful to note that It has now become straightforward to derive our main claim concerning the form of the monodromy matrix: It has matrix elements which are finite difference operators in the multiplicity labels of low order.

Relations to the quantisation of moduli spaces of flat connections
Our next goal will be to explain how the "quantum monodromy" relations (2.35) can be used to define a quantisation of the moduli spaces of flat sl 3 -connections on C = C 0,3 , and more generally on punctured Riemann spheres C 0,n . After recalling the description of these moduli spaces as algebraic varieties, we will recall the definition of the Verlinde line operators V γ , a natural family of operators acting on the spaces of conformal blocks labelled by closed curves γ on C. It was shown in [CGT] that the Verlinde line operators generate representations of the quantised algebras of functions on the moduli spaces of flat connections with V γ representing the quantised counterpart of a trace function L γ associated to γ. Our main goal in this section is to explain why the operators u, v representing elementary building blocks of the operator-valued monodromy matrices M γ (u, v) ı  are closely related to quantum counterparts of a natural higher rank generalisation of the Fenchel-Nielsen coordinates.

Moduli spaces of flat SL(3)-connections
Let us denote the space of holomorphic connections on P 1 \ {z 1 , . . . , z n } of the form The description can be reduced to the case n = 3 by pants decompositions. As a useful set of generators for M B (SL(3), C 0,3 ) one can choose the following combinations of trace functions [La1,La2], see also [CGT]: Let M i , i = 1, 2, 3 be the monodromies around the three punctures Other trace functions L γ can be expressed in terms of these generators using the skein relations.
It was shown in [La1,La2,CGT] that the generators introduced above satisfy two algebraic relations of the form P i {L γ , γ ∈ π 1 (C 0,3 )} = 0, i = 1, 2. (5. 3) The explicit formulae for the polynomials P 1 and P 2 can be found in the references cited above.
The algebra A 0 = Fun(M B (SL(3), C 0,3 )) is a Poisson-algebra with Poisson bracket defined using the Atiyah-Bott symplectic structure. The Poisson structure is algebraic, it closes among the polynomial functions of (N,N , W,W ) and A = ( La2,CGT]. The functions (A i ,Ā i ) i=1,2,3 generate the center of this Poisson algebra.

Quantisation of the moduli spaces of flat connections from CFT
Spaces of conformal blocks carry two natural module structures. One type of module structure is canonically associated to the definition of conformal blocks as solution to the conformal and W-algebra Ward identities, as defined in Section 2. Of particular interest for us is another type of module structure defined in terms of the so-called Verlinde line operators and their generalisations. This subsection will present the definition of the Verlinde line operators, and explain why the algebra of Verlinde line operators is a quantum deformation of the algebra of regular functions on Fun(M B (SL(3), C 0,3 )) [CGT].
In order to define the Verlinde line operators let us note that the identity operator appears in the operator product expansion , with Y n (y) being combinations of the generators of the W 3 algebra, and P ı ∈ C, ı = 1, 2, 3. The chiral vertex operators V l (y) correspond to the other finite-dimensional representations appearing in the tensor product of fundamental and anti-fundamental representation of sl 3 . This shows that a canonical projection P from the space of conformal blocks on C 0,3+2 spanned by the F ρ k (y)ī ı to the space of conformal blocks on C 0,3 spanned by F ρ k can be defined by setting P(F ρ k (y)ī ı ) = δī ı P ı F ρ k . By taking linear combinations of (5.4) one may eliminate the terms containing the chiral vertex operators V l (y), leading to a relation of the form This allows us to define an embedding E of the space of conformal blocks on C 0,3 into the space of conformal blocks on C 0,3+2 by setting E(F ρ k ) = 3 ı=1 E ı F ρ k (y) ıı . As before we find it convenient to regard the functions F ρ k as components of a vector F ρ with respect to a basis labelled by k. We had observed in Section 2 that the analytic continuation of F ρ (y)ī ,ı along γ ∈ π 1 (C 0,3 ) defines an operator M γ on the space of conformal blocks on C 0,3+2 , The composition P • M γ • E will then be an operator from the space of conformal blocks on C 0,3 to itself which may be represented in the form The operators V γ will be called Verlinde line operators. They can be regarded as "quantum" versions of the traces of the monodromy matrices M γ of a flat SL(3)-connection.
Assuming the validity of the mondromy relations (2.35) it was shown in [CGT] that the algebra generated by the operators V γ is isomorphic to the algebra A of quantised functions on M flat (SL(3), C 0,3 ), the moduli space of flat SL(3)-connections on C 0,3 . This means in particular that the operators V γ satisfy deformed versions of the relations (5.3) with polynomials P i differing from the polynomials P i in (5.3) by fixing an operator ordering and -dependent deformations of the coefficients. The algebra generated by the Verlinde line operators can be identified with the skein algebra defined from quantum group theory using the Reshetikhin-Turaev construction [CGT]. 3 It follows that the operators V γ satisfy qcommutation relations reproducing the Poisson brackets of the corresponding trace functions to leading order in .
We had observed in Section 2.5 that the elements of the quantum monodromy matrices can be expressed as functions of two operators u and v. It follows from (5.7) that the same is true for the Verlinde line operators V γ . We are now going to argue that the classical limits of u and v are closely related to natural higher rank analogs of the Fenchel-Nielsen coordinates.

The notion of Fenchel-Nielsen type coordinates
As a preparation let us first discuss defining properties of a class of coordinate system which will be called coordinates of Fenchel-Nielsen type. In order to motivate our proposal we will briefly review the case of flat sl 2 connections.
Recall that useful coordinates on M B (SL(2), C 0,4 ) are the trace coordinates, associating to a closed curve γ the trace of the holonomy along γ. For N = 2, C = C 0,4 = P 1 \ {z 1 , z 2 , z 3 , z 4 } one may, for example, introduce the holonomies M i around the punctures z i and define the following set of trace functions: In order to describe M B (SL(2), C 0,4 ) as an algebraic Poisson-variety let us introduce the polynomial P 0,4 (L s , L t , L u ) : = L 1 L 2 L 3 L 4 + L s L t L u + L 2 s + L 2 t + L 2 u + L 2 1 + L 2 2 + L 2 3 + L 2 4 (5.11a) The polynomial P 0,4 allows us to write both the equation describing M B (SL(2), C 0,4 ) as an algebraic variety and the Poisson brackets among the regular functions L s , L t and L u in the following form: A set of Darboux coordinates λ, κ for M B (SL(2), C 0,4 ) can then be defined by parameterising L s , L t and L u in terms of coordinates u, v as with coefficient functions C t (u), C u (u), = −1, 0, 1 explicitly given as . These coordinates are relatives of the coordinates previously introduced in [Ji, NRS]. The coordinates defined by representing u, v as u = e iλ , v = e iκ are closely related to the Fenchel-Nielsen coordinates from Teichmüller theory. Restricting λ and κ to the flat PSL(2, R)connections coming from the uniformisation of Riemann surfaces one finds purely imaginary values of these coordinates. 2iλ will coincide with the geodesic length along γ s , and iκ will differ from the Fenchel-Nielsen twist parameter τ only by a simple function of λ, iκ = τ + η(λ). Compared to the Fenchel-Nielsen coordinates the parametrisation introduced above has the advantage that the expressions for the trace functions in terms of λ, κ have more favourable analytic properties. In the case of the Fenchel-Nielsen coordinates one would find square-roots in these expressions.
Requirement (R) means that coordinates of Fenchel-Nielsen type are reflecting the Poissonalgebraic structure of M B (C) in an optimal way. Note that this property is shared by the Fock-Goncharov coordinates [FG1] for M B (C). The main feature referred to in the terminology "Fenchel-Nielsen type" is (P), shared by the Fenchel-Nielsen coordinates for the Teichmüller spaces. Compatibility with a pants decomposition means that appropriate subsets of the coordinates would represent coordinate systems for the surfaces obtained by cutting the surface C along any simple closed curve defining the pants decomposition. The Fock-Goncharov coordinates for M B (C) do not have this property.

Quantum Fenchel-Nielsen type coordinates
In order to relate the operators u and v to quantum counterparts of Fenchel-Nielsen type coordinates for M B (SL(3), C 0,3 ) one mainly needs to observe the following simple consequence of our calculation of quantum monodromy matrices.

(L)
The operators V γ generating a representation of A , = b 2 , on the spaces of conformal blocks can be represented as Laurent polynomials V γ (u, v) of u and v.
In order to support this interpretation let us consider the Poisson-algebra R 0 generated by rational functions of two variables u, v with Poisson-bracket {u, v} = −2π uv. The classical limit A 0 of the algebras A considered above is a Poisson-algebra which can be identified with the sub-algebra of R 0 generated by the Laurent polynomials V 0 γ (u, v) The logarithmic coordinates µ = −i log u, ν = −i log v have Poisson bracket {µ, ν} = 2π, identifying them as Darbouxcoordinates for A 0 . It follows from our observations above that mapping V 0 γ to L γ defines an isomorphism of Poisson algebras. This means in particular that the functions V 0 γ (u, v) satisfy the relations (5.8), and have Poisson brackets that reproduce the Poisson brackets of the trace functions L γ when L γ is replaced by V 0 γ . These properties are completely analogous to the description of M B (SL(2), C 0,4 ), motivating us to call the coordinates µ, ν Fenchel-Nielsen type coordinates for M B (SL(3), C 0,3 ).

Yang's functions and solutions to the Riemann-Hilbert problem from classical limits
We will now point out further relations between Toda CFT and the theory of flat connections on C 0,3 arise in two different limits. The standard classical limit of Toda field theory corresponds to b → 0. Assuming that the vectors α i are of the form α i = −b −1 η i for i = 1, 2, 3 one may expect that the conformal blocks behave in the limit b → 0, α i = −η i /b, k = µ/2πb 2 , with η i , i = 1, 2, 3, and µ finite, as follows: where χ ı (y), andχī(y), ı,ī = 1, 2, 3, are three linearly independent solutions of the differential equations Differential equations of this form are closely related to a special class of flat sl 3 -connections on C 0,3 called opers. The action of the operator-valued matrices M γ (u, v) on F ρ (y 0 , y)ī ı turns into the multiplication with the ordinary matrices M γ (u, v) with u = e iµ and v = e iν in the limit (5.12), where The matrices M γ (u, v) represent the monodromy of the differential equation (5.13). Equation (5.15) describes the image of Op sl 3 (C 0,3 ) within M B (SL(3), C 0,3 ) under the holonomy map. It therefore represents a higher rank analog of the generating function of the variety of sl 2 -opers proposed to represent the Yang's function for the quantised Hitchin system in [T10,NRS,T17].
Another interesting limit is b → i, where c = 2. In this case one finds q 2 = 1, so that the operators u and v commute with each other. The transformation G ρ (y 0 , y; κ, λ)ī ı = n∈Z e πi nλ e − π 2 i n 2 F ρ κ+n (y 0 , y)ī ı (5.16) diagonalises u and v simultaneously with eigenvalues u = e −2πi κ and v = ie −πi(κ+λ) , respectively. It follows that the operator-valued monodromy matrices M γ = M γ (u, v) are transformed into ordinary matrices M γ (u, v). The conformal blocks G ρ represent solutions to the Riemann-Hilbert problem to construct multi-valued analytic functions on C 0,3 with monodromies fixed by the data λ, κ. Comparison with the main result of [ILT] suggests that G ρ (κ, λ) = n∈Z e πi nλ e − π 2 i n 2 F ρ κ+n (5.17) may serve as a partial replacement for the isomondromic tau function in this context.

Conclusions and outlook
Our main result is a remarkably simple relation between the screening charges of the free field representation for the W 3 algebra and the Fenchel-Nielsen type coordinates M B (SL(3), C 0,3 ). This direct link between the free field realisation of CFTs and the (quantum) geometry of flat connections on Riemann surfaces seems to offer a key towards a more direct understanding of the relations between these two subjects. We believe that these relations deserve to be investigated much more extensively.

More punctures, higher rank
The generalisation to conformal blocks on surfaces with higher number of punctures should be straightforward. A simple counting of variables indicates that Fenchel-Nielsen type coordinates for surfaces C 0,n with arbitrary n can be obtained from the Fenchel-Nielsen coordinates of all subsurfaces of type C 0,3 and C 0,4 appearing in a pants decomposition of C 0,n . This paper has introduced all ingredients necessary to compute the monodromies of degenerate fields on any punctured Riemann sphere. Using these ingredients within the set-up used in [ILT] it is not hard to see that the operator-valued monodromy matrices can be represented as Laurent polynomials of basic shift and multiplication operators which can be related to Fenchel-Nielsen type coordinates in the same way as before.
In the case of Toda CFTs associated to Lie algebras of higher ranks an interesting feature not discussed here will play a more important role. Instead of only three types of screening charges, here denoted Q 1 , Q 2 and Q 12 one will in the case of sl N Toda CFT have 1 2 (N − 1)N screening charges Q i,j , j > i. These screening charges are in one-to-one correspondence with the positive roots of sl N , and can be constructed as multiple q-commutators of the simple screening charges Q i,i+1 associated to the simple roots. This implies that different ordering prescriptions for the positive roots of sl N will correspond to different bases in the space of conformal blocks on C 0,3 . The physical correlation functions in Toda CFT can not depend on such choices. This invariance may represent a supplement to the crossing symmetry conditions exploited in the conformal bootstrap approach which may be potentially useful.

Continuous bases for spaces of conformal blocks in Toda field theories
The conformal blocks defined above generate a subspace of the space of conformal blocks which is sufficient to describe the holomorphic factorisation for the subset of the physical correlation functions in sl 3 Toda CFT admitting multiple integral representations of Dotsenko-Fateev type. For the case of Liouville theory it is known that one needs to consider continuous families of conformal blocks in order to construct generic correlation functions in a holomorphically factorised form. The same qualitative feature is expected to hold in Toda CFTs associated to higher rank Lie algebras. Based on the experiences from Liouville theory one may anticipate some essential features of the generalisation from the discrete families of conformal blocks studied in this paper to the cases relevant for generic Toda correlation functions.

Continuation in screening numbers
Recall that the conformal blocks provided by the free field representation can be labelled by the so-called screening numbers, the numbers of screening currents integrated over in the multipleintegral representations. We conjecture that conformal blocks appearing in generic Toda correlation functions are analytic functions of a set of parameters which restrict to the conformal blocks coming from the free field construction when the parameters are specialised to a discrete set of values labelled by the screening numbers.
In order to support this conjecture let us offer the following observations. A variant of the free field construction can be used to construct continuous families of conformal blocks. For real and sufficiently small values of the parameter b one may define, following [T01], operators Q i (γ), i = 1, 2 by integrating along contours γ supported on the unit circle. Following the discussion of the Virasoro case in [T03] one may argue that the operators Q i (γ) are densely defined and can be made positive self-adjoint with respect to the natural scalar product on the direct sum of unitary Fock modules S dβ V β , with S being the set of vectors β of the form β = Q + iP , P ∈ R 2 , Q = ρ W (b + b −1 ), ρ W being the Weyl vector of sl 3 . For purely imaginary values of n 1 , n 12 and n 2 one may therefore define unitary operators denoted by Q n 2 1 (γ), Q n 12 12 (γ), and Q n 1 2 (γ) Conformal field theory Quantum group theory Moduli space Liouville theory Modular double PSL(2, R)-connections, of U q (SL(2, R)) Fuchsian component Table 4: Relations between conformal field theory, quantum group theory and moduli spaces of flat connections relevant for the case of Liouville theory. The Fuchsian component in the moduli space of flat PSL(2, R)-connections is the connected component formed by all connections having holonomies ρ : π 1 (C) → PSL(2, R) such that H/ρ is a Riemann surface of the same topological type as C.
by simply taking Q n i i (γ) = exp(n i log(Q i (γ)), for i = 1, 12, 2. We conjecture that Q n i i (γ) can be defined for even more general values of n i by analytic continuation.
This means in particular that the definition of the conformal blocks F ρ k given in Section 2.5 can be generalised to provide families of conformal blocks for C 0,3 labelled by a set of continuous parameters ρ = (α 3 , α 2 , α 1 ) and k. We conjecture that the dependence of F ρ k on ρ and k is meromorphic. This is consistent with, and supported by the fact that the operator-valued monodromy matrices calculated in this paper have an obvious analytic continuation which is entire analytic in the parameter k and meromorphic in ρ. If the conformal blocks on C 0,3 have an analytic continuation in these parameters as conjectured, the quantum monodromies defined from the analytically continued conformal blocks must coincide with the obvious analytic continuation of the monodromy matrices calculated in this paper.

Relation to quantum group theory and higher Teichmüller theory
There are further encouraging hints that the description of the continuous families of Virasoro conformal blocks relevant for Liouville theory admits a higher rank generalisation representing the analytic continuation of the conformal blocks considered in this paper. Let us recall that the conformal blocks needed for generic correlation functions come in families associated to a continuous series of Virasoro representations [ZZ]. The monodromies of conformal blocks can be represented in term of the 6j-symbols associated to a continuous family of representations of the non-compact quantum group U q (SL(2, R)) [T01]. The relevant family of quantum group representations [PT1,PT2] is distinguished by remarkable positivity properties closely related [BT1] to the phenomenon of modular duality [F95,PT1,F99]. All this is closely related to the quantisation of the Teichmüller spaces [CF,Ka1]: It turns out that the 6j-symbols of U q (SL(2, R)) mentioned above describe the change of pants decompositions on the four-holed sphere in quantum Teichmüller theory [T03,NT]. The situation is summarised in Table 4.
The mathematical structures appearing in Liouville theory have generalisations associated to Lie algebras of higher ranks. The higher Teichmüller theories [FG1] can be quantised [FG2].
A key feature facilitating the quantisation of the higher Teichmüller spaces is the positivity of the coordinates introduced in [FG1] when restricted to the so-called Teichmüller component, a connected component in the moduli spaces of flat PSL(N, R)-connections generalising the Fuchsian component isomorphic to the Teichmüller spaces for N = 2 [Hi]. This generalises similar properties of the ordinary Teichmüller spaces when described in terms of the coordinates going back to [Pe].
The main new feature in the higher rank cases is the appearance of non-trivial spaces of threepoint conformal blocks. We expect these spaces to be isomorphic to the multiplicity spaces of Clebsch-Gordan maps for the positive representations of U q (SL(N, R)). The study of the positive representations of U q (SL(N, R)), N > 2, was initiated in [FI]. The Clebsch-Gordan maps for this class of representations of U q (SL(N, R)) have recently been constructed in [SS2]. For the case N = 3 one gets bases for the space of Clebsch-Gordan maps labelled by one continuous parameter. The constructions in [SS1,SS2] reveal a direct connection between the positive representations of U q (SL(N, R)) and the quantised higher Teichmüller theories generalising what was found for N = 2 in [Ka2,NT].
It should be possible to establish a link between the continuous families of conformal blocks introduced above and a suitable basis for the space of Clebsch-Gordan maps between positive representations of U q (SL(3, R)) by generalising the constructions in [T01] for the case of Liouville theory. The positivity of the screening charges together with the direct relation between Clebsch-Gordan maps and free field representation established in the paper suggests that there is a scalar product on the space of conformal blocks on C 0,3 with ρ = (α 3 , α 2 , α 1 ), α i = Q+iP i for which the conformal blocks F ρ k with k ∈ iR generate a basis.
The e i label can be omitted since it is redundant, being implied by the difference of weights. Similarly, the action of generators f i is represented by

B. Clebsch-Gordan coefficients
Here we sketch how to compute the Clebsch-Gordan coefficients which enter the construction of chiral vertex operators in section 3.2.2, starting from equation (3.19) reproduced below (B.1) To lighten notation, let us drop the first column in the notation (·| · ·), whose first column will be everywhere ( λ 3 0 ) and where λ 3 is fully determined by the other weights.
Inductive step: Let us consider a pair of equations labeled by two integers (s 1 , s 2 ) which enter the weight ν 1 = λ 1 − s 1 e 1 − s 2 e 2 < λ 1 . If we denote by w the sum w = s 1 + s 2 of these labels, which enter the CGC under the weight λ 1 ( λ 2 λ 1 n sm ) k , λ 2 λ 1 n s m+1 k , (B.5) then these coefficients are determined by the coupled pair of linear equations with label (s 1 , s 2 ) in terms of coefficients with w = w − 1 and for which ν 1 = λ 1 − s 1 e 1 − s 2 e 2 + e i > ν 1 . It is thus possible to set up an induction procedure through the pair of equations were not invertible. This is however not the case for generic weights λ 1 , λ 2 , given the explicit expressions (A.16).

C. Braid matrix derivation
In this section we collect a few of the more technical points that occur in the calculations of the braiding of screened vertex operators, as was depicted in Figure 2 and reproduced below in Figure 7. This procedure was described in section 3.3 as an induction process to deform contours around z 2 in Figure 7 to pass to linear combinations of the types of loops illustrated in Figure 8 and then deform the resulting auxiliary loops into a set of basis contours like those in Figure 1. Note: There are two additional technical observations to note here: i) contours associated with the composite screening charge Q 12 = Q 1 Q 2 − qQ 2 Q 1 are most safely treated when considered as linear combinations of contours for simple charges and ii) one must take care of the ordering of the newly generated contours, in particular when they are associated with the screening charges Q 12 and Q 2 . Any reordering of contours implies a permutation of their associated normalization points, which generates braiding factors as illustrated for example by γ dyS 1 (y)V n α (z) = V n+s 1 α (z) vs.

(C.2)
The case where a screening charge Q 2 is added has one additional feature. Changing the order of contours such that the s 2 contour is inside of the loops with labels s 1 means commuting an S 2 screening current past S 1 , which by the definition of composite screening charges generates a current S 12 and thus γ dyS 2 (y)V n α (z) = q n 12 −n 1 V n+s 2 α (z) − q −n 1 [n 1 ]V n+s 12 −s 1 α (z) .

D. Realisation of the generators e i on screened vertex operators
In this section we will indicate how to derive the identification between the action of L −1 on screened vertex operators (3.30), which removes one contour of integration, and the action of e i generators on tensor products of representations (3.31). Using notations similar to the ones introduced in section 3.4.1 for the simpler case of U q (sl 2 ), we find that L −1 , V α n (z) − ∂ z V α n (z) = S(z 0 )(q − q −1 )k −1 eV α n−1 (z), (D.1) by the following considerations. For screening number n = 2, there are only two contours of integration. L −1 acts by the Leibniz rule on each screening current in V α n (z), producing total derivative terms. The contributions from the two boundaries of the integration contour are related by monodromy factors. In this way it is not hard to arrive at the following equation L −1 , V α 2 (z) = dy(1 − e −4πib(α+b) )S(z 0 )S(y) + (1 − e −4πibα )S(y)S(z 0 )V α (z) + ∂ z V α 2 (z) .
In order to see that the boundary terms occurring in the commutators with arbitrary generators of the W -algebra admit a similar representation in terms of the generators e i the main point to observe is that all commutators between screening currents and W -generators can be represented as total derivatives.