The most general λ-deformation of CFTs and integrability

We show that the CFT with symmetry group Gk1×Gk2×⋯×Gkn\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {G}_{k_1}\times {G}_{k_2}\times \cdots \times {G}_{k_n} $$\end{document} consisting of WZW models based on the same group G, but at arbitrary integer levels, admits an integrable deformation depending on 2(n − 1) continuous parameters. We derive the all-loop effective action of the deformed theory and prove integrability. We also calculate the exact in the deformation parameters RG flow equations which can be put in a particularly simple compact form. This allows a full determination and classification of the fixed points of the RG flow, in particular of those that are IR stable. The models under consideration provide concrete realizations of integrable flows between CFTs. We also consider non-Abelian T-duality type limits.


Introduction
Integrability is one of the key features for obtaining exact results in quantum field theory (QFT). The most well-known example where integrability was greatly exploited is that of the maximally supersymmetric gauge theory in four dimensions, that is N = 4 SYM. Thanks to the AdS/CFT correspondence [1] the dynamics of N = 4 SYM can be translated to those of an integrable two-dimensional non-linear σ-model evading thus certain no-go theorems which forbid integrable theories in more than two dimensions. A variety of integrability-based techniques, from the asymptotic Bethe ansatz [2] to the thermodynamic Bethe ansatz [3] and the Y-system [4], were employed in order to determine the planar anomalous dimensions of gauge invariant operators for all values of the't Hooft coupling λ = g 2 Y M N . For further developments on integrability and the AdS/CFT correspondence see [5] and references therein.
It is an intriguing endeavour to formulate new examples of gauge/gravity dualities where although supersymmetry may be broken the integrability of the theories involved will be maintained. The comments above make it clear that one starting point could be new integrable two-dimensional field theories which can serve as the seed for such constructions. One such direction which has attracted attention lately is the λ-deformed models of [6], as well as generalisations in many directions [7][8][9][10].
The aforementioned constructions provide an effective and rather effortless method for obtaining exact results in a general class of two-dimensional QFTs. The starting point of JHEP03(2019)094 the construction is certain conformal field theories (CFTs) of the WZW type perturbed by current bilinear operators with the currents belonging either to the same and/or different groups. The essence of the method relies on the construction of the corresponding all-loop effective actions for the deformed theories by a gauging procedure [6][7][8][9][10]. These effective actions possess non-perturbative symmetries in the space of couplings [9][10][11][12]. 1 Combining low-order perturbation theory with the manifestation of the aforementioned non-perturbative symmetries at the level of correlators one is able to derive exact expressions for them. This program was initiated and implemented in a series of papers in which exact expressions were obtained for many observables. In particular, when gravitational techniques were also implemented the exact in the deformation parameters β-functions were found [9,11,12,14], with expressions for the most general case of anisotropic λdeformations and different WZW levels derived in [15]. 2 In addition, anomalous dimensions of current [20] and primary operators [21,22], as well as three-point correlators involving currents and/or primary operators [21,22] were calculated for the models of [6][7][8][9]. Furthermore, the computation of the C-function of Zamolodchikov [23], exactly in the deformation parameter for the case of isotropic perturbations but to leading order in k, was performed in [24] and further generalized for anisotropic λ-deformations in [25].
In a parallel development of particular importance the extension to the case where the unperturbed CFT of a single WZW model is replaced by a coset CFT was considered in [6,[26][27][28]. The deformed theory was found to be integrable when the coset is chosen to be a symmetric space [26]. Subsequently, the corresponding analysis for the case of supergroups was considered in [27]. Although integrability has not been an essential ingredient in the computation of the β-functions and of the operators anomalous dimensions, in the case of isotropic deformations the above models have been demonstrated to be integrable [6,8,9,[26][27][28]. Note that for the special case of the isotropic deformation based on the SU(2) group the model has been initially proven to be integrable in [29]. Furthermore, integrability was shown to persist in some other cases with more deformation parameters [30,31]. In addition, certain deformed models of low dimensionality have been embedded to type-IIA or type-IIB supergravity [32][33][34][35].
An interesting relation between λ-deformations and η-deformations for group and coset spaces was discussed in [31,[36][37][38][39][40]. In particular, the λ-deformed models are related via Poisson-Lie T-duality, which has been introduced for group spaces in [41] and extended for coset spaces in [42], and appropriate analytic continuations to the η-deformed models. The latter were introduced in [43][44][45] and [46][47][48] for group and coset spaces, respectively. Moreover, the dynamics of scalar fields in some λ-deformed geometries corresponding to coset CFTs has been discussed in [49] while the relation to Chern-Simons theories was discussed in [50,51]. Finally, D-branes regarded as integrable boundary configurations were introduced in the context of λ-deformations in [52]. 1 For isotropic λ-deformations the simplest such non-perturbative symmetry was found initially in [13] using path integral arguments. 2 The exact β-functions for CFTs deformed by operators bilinear in currents and for isotropic cases have been obtained in the past either by field theoretical methods (resummation of the perturbation series or the background field method) [16][17][18][19]. The results are in complete agreement.

JHEP03(2019)094
The models presented in [9,10] have several virtues. They provide concrete realizations of flows between exact CFTs in which the Lagrangian of the theory is known all the way through from the UV to the IR fixed points. Integrability is preserved in the entire flow and therefore these flows are called integrable. Furthermore, the construction of [10] provides the first example in which self-and mutual-interactions between different WZW models are present already at leading order in the deformation parameters. In this work we continue this line of research and construct integrable multi-parameter deformations of CFTs whose Lagrangian and RG flow equations are known to all-orders in perturbation theory. The integrable sector will belong to the most general λ-deformed model which we construct. In this model self-as well as mutual-interactions are present.
The plan of the present paper is as follows: in section 2, we construct the all loop effective action of a general class of models whose UV Lagrangian is the sum of an arbitrary number n of WZW models based on the same group G but at different levels. In general our models depend on n 2 general coupling matrices. In section 3, we consider a consistent truncation of our models that depend on 2(n − 1) couplings, identify the non-perturbative symmetries in the space of couplingsλ i1 andλni and show that the theory is classically integrable by finding the appropriate Lax connection. In section 4, we calculate the exact in the deformation parameters β-functions of our models which can be cast in a particularly simple and compact form. Subsequently, we determine and classify the fixed points of the RG flow in particular those that are IR stable. Finally, in section 5 we consider non-Abelian T-duality type limits in the case when the theory is integrable. In the last section we present our conclusions.

Generalities: Lagrangian and equations of motion
In this section we construct the effective actions of our model and derive the corresponding equations of motion.
Our starting point is to consider group elements g i , i = 1, 2, . . . , n in a semi-simple group G and the corresponding WZW model actions S k i (g i ) at levels k 1 , k 2 , . . . , k n . We add to this the action of n PCMs which are both self and mutually interacting and are constructed using n group elementsg i , i = 1, 2, . . . , n belonging in the same group G. Namely, we have that where the E ij , i, j = 1, 2, . . . , n are generic coupling matrices and the indices a, b = 1, 2, . . . , dim(G). In the spirit of [6,8,10] we gauge the global symmetry acting on the group elements as follows: . . , n. The resulting gauge invariant action reads is the standard gauged WZW action where we have suppressed the group indices. Furthermore, the covariant derivatives were defined asD ±gi = (∂ ± − A (i) ± )g i . One may now fix the gauge in (2.2) choosingg i = to arrive at the following action where we have redefined the coupling matrices appearing in the PCM models as 3 In order to obtain the σ-model and to show integrability in the next section we should integrate out the gauge fields. These are not dynamical and appear only quadratically in the action. We cast the corresponding equations of motion into the following particularly convenient form where as usual we have defined ± , g i ] and where we note that the transpose in λ −T ij , as well as the inverses refer only to the suppressed group indices and not on the space of couplings with indices i, j.
Varying with respect to the group elements one obtains (2.7) To proceed we find it convenient to rescale the gauge fields and definẽ (2.8) 3 To compare with the corresponding actions in [10] (see eqs. (2.1) and (2.4)) corresponding to two WZW models self-and mutually interacting we have

JHEP03(2019)094
Using them and substituting (2.6) in (2.7) we obtain after some algebra that (2.9) To present the σ-model action we introduce representation Hermitian matrices t a and define the currents and primary field to be These will be computed for the particular group elements g i in which case an extra index i will be inserted in the above quantities. We also define the diagonal in the coupling space matrix Then we may solve for the redefined gauge fields using (2.6) to obtaiñ Note that the entries of the matrices λ −T − D and λ −1 − D T are themselves matrices in the group G. Thus, their inversion is to be understood as an inversion in the space of different models keeping in mind that their entries are non-commutative objects.
Finally, substituting the values for the gauge fields in the action (2.4) we obtain the following σ-model This is a particularly compact form resembling the single λ-deformed model of [6]. The above action encompasses all previous ones for λ-deformed models. Furthermore it is the most general action that can be construct using the same group G for all couplings. A further generalization in which a different group is associated to each WZW model can be constructed straightforwardly.
For small entries in the matrices λ ij the action becomes (2.14) Note that λ −1 −1 ij = λ ij since the inverse in λ −1 is taken in the group space, and the dots denote subleading terms in the small λ-expansion.

JHEP03(2019)094
Finally, note that in integrating out the gauge field a dilaton field is generated. The value of the dilaton is proportional to the logarithm of the determinant of the square matrix λ −1 − D T inverted in (2.12). Note that the dimensionally of this matrix is n dim G.
In the next section the inversion of the above matrices in the coupling space will be done explicitly for the case of integrable models.

Integrable deformations
In this section we will attempt to answer the following question. For which choices of the matrix E ab ij appearing in (2.1) is the theory described by (2.13) integrable? To the best of our knowledge we are lacking a general answer to this question even for the case of a single group G. Nevertheless, it has been proven that there are several cases where these theories are integrable for specific choices of the couplings matrices. Among the single λ-deformations the first one is the case of the isotropic λ-deformation, that is when λ ab = λδ ab [6], a second one is the case where the matrix E is of the form E = 1 t (1 − ηR) −1 , where the matrix R satisfies the modified Yang-Baxter equation [44,45]. Another case is that where instead of the group G one has a coset with the coset being a symmetric space [6,26]. 4 Furthermore, integrability has been shown for the models of [8,9] and [10] representing particular cases of self-and mutual-interactions of current algebras based on WZW models.

A truncation of our models
In what follows, we will show that the theory (2.13) is integrable in the case where the matrices λ ij is of the following form  Hence, the corresponding coupling matrices reads As a result it will be shown that one obtains an integrable deformation of the G k 1 × G k 2 × · · · × G kn depending on 2(n − 1) parameters. This provides a partial answer to the question posed in the beginning of this section. We stress that by turning on just one more coupling, in addition to those appearing in (3.2), it will generate other couplings through quantum corrections. Then, the theory will most likely cease to be integrable. However, we have no proof that (3.2) exhausts all integrable cases among our general class of models.

JHEP03(2019)094
Turning to the general action (2.13), after using (3.1) and explicitly inverting the relevant matrix in the coupling space we find that From this point on we will focus for simplicity on the case where the couplings are taken isotropic in the group space, that is One now encounters the following problem. By combining S k i (g i ), i = 2, . . . , n − 1 with the last term in the second line of (3.3) we obtain a sum of WZW model actions S −k i (g −1 i ), i = 2, . . . , n−1 with negative signature (this is also true for the PCM part of the action (2.1) before the gauging procedure is performed). As was pointed out in the case of two interacting WZW models in [10] to remedy this situation one can perform the following redefinitions of the couplings and analytic continuation in the following specified order. First define the hatted couplingŝ and then flip the signs of levels and invert the corresponding group elements as Then the action (3.3) becomes

JHEP03(2019)094
This is the final expression for the all-loop effective action of our σ-model. In what follows we will prove that it is indeed classically integrable. Before that, we comment on the constraints imposed by demanding a non-singular σ-model of Euclidean signature. In order to avoid singularities the couplings λ 11 and λ nn should be such that |λ 11 |, |λ nn | < 1, the reason being that the matrix D is orthogonal and therefore its eigenvalues lie in the unit circle. In addition, similarly to [10], it can be shown that the signature of (3.7) is Euclidean provided that the couplings lie within the ellipsoids defined by It is straightforward to show that the action (3.7) is independently invariant under the following two non-perturbative symmetry operations, as well as under their combination (A) : This is a generalization of the corresponding symmetry found in [10]. Expanding (3.7) and keeping the linear terms in the couplings one obtains (3.10) Notice the last line of (3.7) has disappeared from the small coupling expansion since it is quadratic in the λ's. This term, as well as the full action of (3.7) is generated when quantum corrections at arbitrary order in perturbation theory are taken into account. Note that there is the following interesting truncation of (3.7) which dramatically further simplifies it. If we let λ 11 = λ nn = 0 then we obtain This effective action is exact in λ and is at most quadratic in the couplings and represents only mutual interactions between the various WZW models as depicted. It becomes linear if we make the further truncationλ n1 = 0 orλ i1 = 0.

Proof of integrability
To prove that (3.7) is integrable we write the equations of motion (2.9) for the choice of couplings appearing in (3.1). The first equation of (2.9) can be decomposed to (3.12) while the second one to − , i = 2, . . . , n. The two aforementioned sets of equations are decoupled. We stress that even though in these systems the couplings λ −1 i1 completely disentangled from the λ −1 ni 's this is not the case in the effective action (3.7) where the last term involves both sets of couplings. The reason is that the various gauge fields depend on all group elements g i , i = 1, 2, . . . , n.
Our strategy will be to determine a Lax pair for each set of equations and then show that the charges obtained from the first Lax pair are in involution with those obtained from the second one, thus proving that the theory is integrable. To this end we define the rescaled anti-commutators (3.14) In terms of these redefined quantities the solution to the first set of equations reads − ] * L   .

(3.15)
These imply the existence of a Lax pair which is constructed by first assuming that this is of the form In conclusion we have shown that the equations of motion (3.15) imply the existence of a Lax pair from which an infinite tower of conserved charges can be calculated.
In a similar manner, the second system of equations can be solved for the derivatives of the gauge fields as (3.20) The corresponding Lax pair will be of the form The flatness of this Lax pair is guaranteed when The conserved charges obtained from (3.16) and (3.21) are in involution. This is so because any the gauge fields appearing in (3.16) have zero Poisson bracket with any of the gauge fields appearing in (3.21). To see this one may define, similarly to what was first done in the Hamiltonian treatment of gauged WZW models in [53], the following dressed currents  (i) ± which satisfy the usual Kac-Moody algebra at level k i  (i)

JHEP03(2019)094
Using these definitions and the constraints (2.6) one obtains This concludes the proof that the CFT deformed by the 2(n − 1) deformation parameters λ −1 i1 , i = 1, 2, . . . , n − 1 and λ −1 ni , i = 2, 3, . . . , n is integrable. A final comment is in order. The careful reader may have noticed that the analysis for integrability has been done using the equations of motion (3.12) and (3.13) in which the analytic continuation discussed in (3.6) has not been applied yet. However, one can easily see that the analytic continuation of (3.6) can be straightforwardly be done in the final expressions for the Lax pairs (3.16) and (3.21) once the latter are expressed in terms of the usual gauge fields A (i) ± and the couplingsλ i1 andλ ni appearing in the σ-model (3.7).

The RG flow equations
In this section we will first calculate the running of the couplings in the case of the integrable deformations presented in the previous section. We will need the system of RG flow equations for the running of the couplings of the anisotropic λ-model [12] where t = ln µ 2 , µ being the energy scale and where In what follows each of the capital indices of (4.1) will be split in two, that is A = (ia), where i, j, k = 1, 2, . . . , n − 1 enumerate the different groups while a, b, c = 1, . . . , dim(G) denote as usual group indices. Furthermore, in order to take into account for the different levels k i of the WZW models, one should set k = 1 in (4.1) and rescale the structure constants of each group as f ABC = f (ia)(ib)(ic) = f abc / √ k i . All structure constants that do not have the indices enumerating the different groups equal are set, of course to zero.

JHEP03(2019)094
Finally, the coupling constants matrix λ AB can be read from (3.10). It reads Keeping in mind that the couplings of the first column in (4.3) decouple from the ones in the last row of the same equation we focus on the first column and rename λ 11 by λ 1 and similarlyλ i1 by λ i , where i = 2, 3, . . . , n − 1. 5 Then the non-zero components of N ia,ib jc (Λ) and N 1a,1b ic (Λ T ) are is the same quantity, called ∆ 1 in (3.8), after the renaming of the couplings we mentioned. Furthermore, Notice that N ia,jb kc (Λ) = −N jb,ia kc (Λ) and N ia,jb kc (Λ T ) = −N jb,ia kc (Λ T ). Using the expressions above we finally obtain for the running of the couplings the following formula where We note that despite of the privileged rôle of λ 1 in the action (3.7) (see also (3.10)) the β-functions for all couplings are on equal footing and can be obtained from a single expression (4.7). This fact will allow us to fully determine the fixed points of the RGflow in the next sections. Furthermore, the β-functions in (4.7) are invariant under the non-perturbative symmetry of (3.9). Finally, similar expressions hold for the running of the couplingsλ ni . The same is true for the analysis following in this section.

Fixed points
We now determine the fixed points of the RG flow equations (4.7). Each fixed point belongs to one of the following classes. Each class is characterised by three integers (n 1 , n 2 , n 3 ), obeying the condition n 1 + n 2 + n 3 = n − 1. By n 1 we denote the number of the couplings that are set to zero, that is λ m i = 0, i = 1, 2, . . . , n 1 . These couplings can be distributed randomly among the complete set of the n−1 couplings. By n 2 we denote the number of the couplings that are set to one, that is λ p i = 1, i = 1, 2, . . . , n 2 . Finally, by n 3 we denote the number of couplings that are neither zero nor one, i.e. λ q i = 0, 1, with i = 1, 2, . . . , n 3 . By subtracting pairwise the equations in (4.7) it is straightforward to see that all the λ q i which are neither zero nor one should be equal to each other, that is λ q i = λ * , i = 1, 2, . . . , n 3 . Then the quantity ∆ becomes As a result the vanishing of the bracket in (4.7) implies (4.10) We will assume in the rest of our analysis that n 3 1 since it turns out that this should be the case for the existence of physical IR stable fixed points. Let us note that each of the distinct classes of fixed points characterised by (n 1 , n 2 , n 3 ), that obey the condition n 1 + n 2 + n 3 = n − 1 has (n − 1)! n 1 !n 2 !n 3 ! different members.

The stability matrix
We evaluate the stability matrix for each of the fixed points given by (4.10). This will allow us to determine the relevant and irrelevant directions of each of the fixed points and identify the fixed points which are IR stable. Let's split the index i = (m, a, α) in such a way that λ m = 0 , m = 1, 2, . . . , n 1 , λ a = 1 , a = 1, 2, . . . , n 2 , Then at the FPs we have the following relations Furthermore, the derivatives of ∆ and Z evaluated at the FPs read

JHEP03(2019)094
The non-zero entries of the stability matrix defined as H ij = ∂ j β i | FP are given by (4.14) Note that H αa = 0, so that the stability matrix is not a symmetric one. Clearly the part of this matrix corresponding to the n 1 values λ m = 0 decouples. For the eigenvalue problem for the rest of the stability matrix we clearly obtain a matrix of the form Specifically, we have the matrix elements A ab = H ab − Eδ ab already in diagonal form with E appearing in the left hand side of (4.15) being an eigenvalue of the stability matrix. The matrix C αβ = H αβ − Eδ αβ is of the form c 1 δ αβ − c 2 k β and as a result its determinant is equal to c n 3 −1 1 (c 1 − ζc 2 ). Furthermore, the second matrix in the right hand side of (4.15) is triangular. The vanishing of the determinant of the matrix in (4.15) gives the eigenvalues of the stability matrix, as well as the corresponding degeneracies Deg = n 1 : (4.16) These eigenvalues and their degeneracies will be instrumental below in determining the physical IR stable fixed points of the RG flow.

IR stable fixed points
It turns out that in order to have a Euclidean signature for the metric, the constant ∆ should be positive leading to the following condition As mentioned before we will only consider cases with n 3 1. In addition, physical fixed points cannot have λ 1 = 1 since the σ-model action (3.7) is in that case singular.
We are primarily interested in identifying the physical IR stable points, that is the ones which have all eigenvalues of the stability matrix positive. The condition (4.17) guaranties that H 0 > 0 and H 1 > 0. However, the third eigenvalue in (4.16) is negative, i.e. H * < 0

JHEP03(2019)094
and therefore we should necessarily have n 3 = 1 so that this eigenvalue is non-existing. Finally, the positivity of the last eigenvalue H 1 * is also guaranteed by (4.17).
The condition n 3 = 1 means that there is a single coupling that is not zero or one. Keeping in mind that λ 1 = 1, there are two cases. In the first one let λ 1 = 0 and a single λ a = λ * . Then also ζ = k a . The condition (4.17) implies that ρ < k 1 < ρ + k a . The second case is when λ 1 = λ * . In this case λ 1 = 1 − ρ/k 1 which holds automatically since according to (4.17). Also, in this case ζ = k 1 . We conclude that there is a multitude of fixed points the number of which depends on the relative ordering of the levels k i . It would be certainly interesting to study the structure and properties of the corresponding CFTs.

The non-Abelian T-duality limit
When some of the λ's approach unity then we get a singularity in the manifold. Then a zoom in procedure maybe applied as in [6]. We will not discuss this for the most general action (2.13) but for (3.7) corresponding to integrable models.
In that case near λ 11 = 1 or near λ nn = 1 we get a singularity in the manifold. However, one may zoom in by taking simultaneously the large k 1 and k n limits. Then we expand for k 1 , k n 1 as follows where ζ 1 and ζ n are new coupling parameters. This leads to the following expressions for the currents and the operators D i In this limit the action (3.7) becomes Note that Euclidean signature imposes a constraint on the parameters

JHEP03(2019)094
This σ-model represents the interaction of (n−2) WZW models for a group G with two non-Abelian T-duals of the PCM for the same group which in turn interact among themselves. We note that the σ-model (5.3) being the limit of (3.7) is itself integrable. The original action whose non-Abelian T-dual is (5.3) is that of interacting PCMs. An interesting limit of (5.3) arises if we take the remaining levels k i → ∞ as . . , n − 1 . (5.5) This represents the mutual interactions of (n−2) dim G free fields with the two non-Abelian T-duals of PCM for a group G and is also integrable. We may clearly simply even further by consistently taking ζ i1 = 0 or ζ ni = 0.

Discussion and future directions
We have construct the all loop effective action of a general class of models whose UV Lagrangian is the sum of an arbitrary number n of WZW models based on the same group G, but at different levels. Although the complete effective action can be quite involved, at the linear level the theory is driven away from the conformal point by operators bilinear in the WZW currents. These current bi-linears involve currents belonging to both the same and different CFTs. Hence we have self-as well as mutual interactions of current algebra theories in their most general form. In general our models depend on n 2 general coupling matrices. We considered a consistent truncation of our models that depends on 2(n − 1) couplings and showed that the theory is classically integrable by finding the appropriate Lax connection. Turning on even one more coupling will generate other couplings through quantum corrections and the theory will most likely cease to be integrable. Subsequently, we proved that the theory possesses certain non-perturbative symmetries in the space of couplingsλ i1 andλni and calculated the exact in the deformation parameters β-functions of our models which can be cast in a particularly simple and compact form. This fact allowed us to fully determine and classify the fixed points of the RG flow in particular those that are IR stable. Last but not least we consider non-Abelian T-duality type limits in the case when the theory is integrable.

JHEP03(2019)094
A number of open questions remain to be addressed. Given that our models provide concrete realizations of integrable flows between exact CFTs it would be interesting to elucidate the nature and symmetries of the corresponding IR stable CFTs. Furthermore, one could exploit the aforementioned non-perturbative symmetries that our models enjoy to compute the all-loop anomalous dimensions of current operators, as well as that of primary operators in a similar manner to that in [7,14,[20][21][22]24]. One could also calculate the exact in the deformation parameters C-function of the models presented here as was done in [24] for simpler cases. In that respect the general results in [24] should be a useful starting point.
Another direction would be to consider the case where more of the couplings are nonzero, as well as the case where each of the WZW models is based on a different group. It is notable that in the latter case all formulae of section 2 will still be valid after slight modifications. In addition, one could search for integrable deformations in these more general cases. Compared to the integrable models presented in this work, we expect an even richer structure of the RG equations to be unveiled. Recently, a class of integrable models consisting of N coupled principal chiral models each with a WZW term and based on the same group G was presented in [54]. The construction was based on an association of integrable field theories with affine Gaudin models having an arbitrary number of sites. It would be interesting to see if these models bear any relation to the ones constructed in this work (see also [10]) or to the models with N sites firstly presented in [7].
Our approach is valid to all orders in the deformation parameters but only to leading order in the 1/k expansion, i.e. for large levels of the WZW action. One may wonder of how corrections of order 1/k 2 or higher may be computed. Although constructing the 1 k corrections to the effective action of the present work seems a hard task, one may resort to the higher order derivative β-functions of the σ-model. This will be helpful in computing corrections to the β-functions of the deformation parameters using gravitational methods and also in discovering the corrected version of non-perturbative symmetries of the type (3.9). In addition, we may compute the subleading terms in the 1/k expansion of the β-functions, as well as of the anomalous dimensions of operators using CFT techniques. Finally, one could try to embed our models to solutions of type-IIB or type-IIA supergravity (for n = 2, 3 and small rank semisimple groups) and/or construct the corresponding ηdeformed integrable models.