Classification of Non-Affine Non-Hecke Dynamical R-Matrices

A complete classification of non-affine dynamical quantum R-matrices obeying the Gln(C)-Gervais-Neveu-Felder equation is obtained without assuming either Hecke or weak Hecke conditions. More general dynamical dependences are observed. It is shown that any solution is built upon elementary blocks, wh...

Повний опис

Збережено в:
Бібліографічні деталі
Дата:2012
Автори: Avan, J., Billaud, B., Rollet, G.
Формат: Стаття
Мова:English
Опубліковано: Інститут математики НАН України 2012
Назва видання:Symmetry, Integrability and Geometry: Methods and Applications
Онлайн доступ:http://dspace.nbuv.gov.ua/handle/123456789/148388
Теги: Додати тег
Немає тегів, Будьте першим, хто поставить тег для цього запису!
Назва журналу:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Цитувати:Classification of Non-Affine Non-Hecke Dynamical R-Matrices / J. Avan, B. Billaud, G. Rollet // Symmetry, Integrability and Geometry: Methods and Applications. — 2012. — Т. 8. — Бібліогр.: 25 назв. — англ.

Репозитарії

Digital Library of Periodicals of National Academy of Sciences of Ukraine
id irk-123456789-148388
record_format dspace
spelling irk-123456789-1483882019-02-19T01:25:44Z Classification of Non-Affine Non-Hecke Dynamical R-Matrices Avan, J. Billaud, B. Rollet, G. A complete classification of non-affine dynamical quantum R-matrices obeying the Gln(C)-Gervais-Neveu-Felder equation is obtained without assuming either Hecke or weak Hecke conditions. More general dynamical dependences are observed. It is shown that any solution is built upon elementary blocks, which individually satisfy the weak Hecke condition. 2012 Article Classification of Non-Affine Non-Hecke Dynamical R-Matrices / J. Avan, B. Billaud, G. Rollet // Symmetry, Integrability and Geometry: Methods and Applications. — 2012. — Т. 8. — Бібліогр.: 25 назв. — англ. 1815-0659 2010 Mathematics Subject Classification: 16T25; 17B37; 81R12; 81R50 DOI: http://dx.doi.org/10.3842/SIGMA.2012.064 http://dspace.nbuv.gov.ua/handle/123456789/148388 en Symmetry, Integrability and Geometry: Methods and Applications Інститут математики НАН України
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection DSpace DC
language English
description A complete classification of non-affine dynamical quantum R-matrices obeying the Gln(C)-Gervais-Neveu-Felder equation is obtained without assuming either Hecke or weak Hecke conditions. More general dynamical dependences are observed. It is shown that any solution is built upon elementary blocks, which individually satisfy the weak Hecke condition.
format Article
author Avan, J.
Billaud, B.
Rollet, G.
spellingShingle Avan, J.
Billaud, B.
Rollet, G.
Classification of Non-Affine Non-Hecke Dynamical R-Matrices
Symmetry, Integrability and Geometry: Methods and Applications
author_facet Avan, J.
Billaud, B.
Rollet, G.
author_sort Avan, J.
title Classification of Non-Affine Non-Hecke Dynamical R-Matrices
title_short Classification of Non-Affine Non-Hecke Dynamical R-Matrices
title_full Classification of Non-Affine Non-Hecke Dynamical R-Matrices
title_fullStr Classification of Non-Affine Non-Hecke Dynamical R-Matrices
title_full_unstemmed Classification of Non-Affine Non-Hecke Dynamical R-Matrices
title_sort classification of non-affine non-hecke dynamical r-matrices
publisher Інститут математики НАН України
publishDate 2012
url http://dspace.nbuv.gov.ua/handle/123456789/148388
citation_txt Classification of Non-Affine Non-Hecke Dynamical R-Matrices / J. Avan, B. Billaud, G. Rollet // Symmetry, Integrability and Geometry: Methods and Applications. — 2012. — Т. 8. — Бібліогр.: 25 назв. — англ.
series Symmetry, Integrability and Geometry: Methods and Applications
work_keys_str_mv AT avanj classificationofnonaffinenonheckedynamicalrmatrices
AT billaudb classificationofnonaffinenonheckedynamicalrmatrices
AT rolletg classificationofnonaffinenonheckedynamicalrmatrices
first_indexed 2025-07-12T19:17:40Z
last_indexed 2025-07-12T19:17:40Z
_version_ 1837469914447740928
fulltext Symmetry, Integrability and Geometry: Methods and Applications SIGMA 8 (2012), 064, 45 pages Classification of Non-Affine Non-Hecke Dynamical R-Matrices Jean AVAN †, Baptiste BILLAUD ‡ and Geneviève ROLLET † † Laboratoire de Physique Théorique et Modélisation, Université de Cergy-Pontoise (CNRS UMR 8089), Saint-Martin 2, 2, av. Adolphe Chauvin, F-95302 Cergy-Pontoise Cedex, France E-mail: avan@u-cergy.fr, rollet@u-cergy.fr ‡ Laboratoire de Mathématiques “Analyse, Géometrie Modélisation”, Université de Cergy-Pontoise (CNRS UMR 8088), Saint-Martin 2, 2, av. Adolphe Chauvin, F-95302 Cergy-Pontoise Cedex, France E-mail: bbillaud@u-cergy.fr Received April 24, 2012, in final form September 19, 2012; Published online September 28, 2012 http://dx.doi.org/10.3842/SIGMA.2012.064 Abstract. A complete classification of non-affine dynamical quantum R-matrices obeying the Gln(C)-Gervais–Neveu–Felder equation is obtained without assuming either Hecke or weak Hecke conditions. More general dynamical dependences are observed. It is shown that any solution is built upon elementary blocks, which individually satisfy the weak Hecke condition. Each solution is in particular characterized by an arbitrary partition {I(i), i ∈ {1, . . . , n}} of the set of indices {1, . . . , n} into classes, I(i) being the class of the index i, and an arbitrary family of signs (εI)I∈{I(i), i∈{1,...,n}} on this partition. The weak Hecke- type R-matrices exhibit the analytical behaviour Rij,ji = f(εI(i)ΛI(i) − εI(j)ΛI(j)), where f is a particular trigonometric or rational function, ΛI(i) = ∑ j∈I(i) λj , and (λi)i∈{1,...,n} denotes the family of dynamical coordinates. Key words: quantum integrable systems; dynamical Yang–Baxter equation; (weak) Hecke algebras 2010 Mathematics Subject Classification: 16T25; 17B37; 81R12; 81R50 1 Introduction The dynamical quantum Yang–Baxter equation (DQYBE) was originally formulated by Gervais and Neveu in the context of quantum Liouville theory [18]. It was built by Felder as a quanti- zation of the so-called modified dynamical classical Yang–Baxter equation [5, 6, 15, 16], seen as a compatibility condition of Knizhnik–Zamolodchikov–Bernard equations [7, 8, 17, 19, 22]. This classical equation also arose when considering the Lax formulation of the Calogero–Moser [9, 23] and Ruijsenaar–Schneider model [24], and particularly its r-matrix [1, 4]. The DQYBE was then identified as a consistency (associativity) condition for dynamical quantum algebras. We introduce A as the considered (dynamical) quantum algebra and V as either a finite-dimensional vector space V or an infinite-dimensional loop space V = V ⊗C[[z, z−1]]. We define the objects T ∈ End(V ⊗A) as an algebra-value d matrix encoding the generators of A and R ∈ End(V ⊗V) as the matrix of structure coefficients for the quadratic exchange relations of A R12(λ+ γhq)T1(λ− γh2)T2(λ+ γh1) = T2(λ− γh1)T1(λ+ γh2)R12(λ− γhq). (1.1) As usual in these descriptions the indices “1” and “2” in the operators R and T label the respective so-called “auxiliary” spaces V in V ⊗ V). In addition, when the auxiliary spaces are mailto:avan@u-cergy.fr mailto:rollet@u-cergy.fr mailto:bbillaud@u-cergy.fr http://dx.doi.org/10.3842/SIGMA.2012.064 2 J. Avan, B. Billaud and G. Rollet loop spaces V = V ⊗ C[[z, z−1]], these labels encapsulate an additional dependence as formal series in positive and negative powers of the complex variables z1 and z2, becoming the so-called spectral parameters when (1.1) is represented, e.g. in evaluation form. Denoting by N∗n the set {1, . . . , n}, for any n ∈ N∗ = N \ {0}, both R and T depend in addition on a finite family (λi)i∈N∗n of c-number complex “dynamical” parameters understood as coordinates on the dual algebra h∗ of a n-dimensional complex Lie algebra h. The term “dynamical” comes from the identification of these parameters in the classical limit as being the position variables in the context of classical Calogero–Moser or Ruijsenaar–Schneider models. We shall consider here only the case of a n-dimensional Abelian algebra h. Non-Abelian cases were introduced in [25] and extensively considered e.g. in [10, 11, 12, 14]. Following [14], in addition with the choosing of a basis (hi)i∈N∗n of h and its dual basis (hi)i∈N∗n , being the natural basis of h∗, we assume that the finite vector space V is a n-dimensional diago- nalizable module of h, hereafter refereed as a Etingof-module of h. That is: V is a n-dimensional vector space with the weight decomposition V = ⊕ µ∈h∗ V [µ], where the weight spaces V [µ] are irreducible modules of h, hence are one-dimensional. The operator R is therefore represented by an n2 × n2 matrix. This allows to understand the notation Ta(λ + γhb), for any distinct labels a and b: λ is a vector in h∗ and hb denotes the canonical element of h⊗ h∗ with a natural action of h on any given vector of V . As a matter of fact, for example a = 1 and b = 2, this yields the usual vector shift by γh2 defined, for any v1, v2 ∈ V as T1(λ+ γh2)v1 ⊗ v2 = T1(λ+ γµ2)v1 ⊗ v2, where µ2 is the weight of the vector v2. The shift, denoted γhq, is similarly defined as resulting from the action on hb of φ⊗1, where φ: h −→ A is an algebra morphism, 1 being the identity operator in the space V . If (1.1) is acted upon by 1⊗ 1⊗ ρH , where ρH is a representation of the quantum algebra A on a Hilbert space H assumed also to be a diagonalizable module of h, then ρH(hq) acts naturally on H (in particular on a basis of common eigenvectors of h assuming the axiom of choice) yielding also a shift vector in h∗. Requiring now that the R-matrix obey the so-called zero-weight condition under adjoint action of any element h ∈ h [h1 + h2, R12] = 0 allows to establish that the associativity condition on the quantum algebra (1.1) implies as a consistency condition the so-called dynamical quantum Yang–Baxter algebra for R R12(λ+ γh3)R13(λ− γh2)R23(λ+ γh1) = R23(λ− γh1)R13(λ+ γh2)R12(λ− γh3). (1.2) Using the zero-weight condition allows to rewrite (1.2) in an alternative way which we shall consider from now on; R12(λ+ 2γh3)R13(λ)R23(λ+ 2γh1) = R23(λ)R13(λ+ 2γh2)R12(λ), (DQYBE) where the re-definition Rab −→ R′ab = Ad exp γ(h ·da+h ·db)Rab is performed, h ·d denoting the differential operator n∑ i=1 hi∂λi . Due to the zero-weight condition on the R-matrix, the action of this operator yields another c-number matrix in End(V ⊗ V) instead of the expected difference operator-valued matrix. Note that it may happen that the matrix R be of dynamical zero-weight, i.e. [h · da + h · db, Rab] = 0, in which case R′ = R. Classification of Non-Affine Non-Hecke Dynamical R-Matrices 3 Early examples of solutions in this non-affine case have been brought to light under the hypothesis that R obeys in addition a so-called Hecke condition [20]. The classification of Hecke type solutions in the non-affine case has been succeeded for a long time starting with the pioneering works of Etingof et al. [13, 14]. It restricts the eigenvalues of the permuted R-matrix Ř = PR, P being the permutation operator of vector spaces V ⊗ 1 and 1⊗ V , to take only the value % on each one-dimensional vector space Vii = Cvi ⊗ vi, for any index i ∈ N∗n, and the two distinct values % and −κ on each two-dimensional vector space Vij = Cvi ⊗ vj ⊕ Cvj ⊗ vi, for any pair of distinct indices (i, j) ∈ (N∗n)2, (vi)i∈N∗n being a basis of the space V . The less constraining, so-called “weak Hecke” condition, not explored in [14], consists in assuming only that the eigenvalue condition without assumption on the structure of eigenspaces. In other words, one only assumes the existence of two c-numbers % and κ, with % 6= −κ, such that (Ř− %)(Ř+ κ) = 0. We shall not assume a priori any Hecke or weak Hecke condition in our discussion. However, an important remark is in order here. The weak Hecke condition is understood as a quantization of the skew-symmetry condition on the classical dynamical r-matrices r12 = −r21 [14]. It must be pointed out here that the classical limit of DQYBE is only identified with the consistent associativity condition for the “sole” skew-symmetric part a12 − a21 of a classical r-matrix parametrizing the linear Poisson bracket structure of a Lax matrix for a given classical integrable system {l1, l2} = [a12, l1]− [a21, l2]. Only when the initial r-matrix is skew-symmetric do we then have a direct connection between classical and quantum dynamical Yang–Baxter equation. Dropping the weak Hecke condition in the quantum case therefore severs this link from classical to quantum Yang–Baxter equation and may thus modify the understanding of (1.2) as a deformation by a parameter ~ of a classical structure. Nevertheless it does not destroy any of the characteristic quantum structures: copro- duct, coactions, fusion of T -matrices and quantum trace formulas yielding quantum commuting Hamiltonians, and as such one is perfectly justified in considering a generalized classification of a priori non-weak Hecke solutions in the context of building new quantum integrable systems of spin-chain or N -body type. The issue of classifying non-affine R-matrices, solutions of DQYBE, when the (weak) Hecke condition is dropped, already appears in the literature [21], but in the very particular case of Gl2(C) and for trigonometric behavior only. A further set of solutions, in addition to the expected set of Hecke-type solutions, is obtained. In the context of the six-vertex model, these solutions are interpreted as free-fermion-type solutions, and show a weak Hecke-type, but non-Hecke-type, behavior R12,21 = f(λ1 + λ2), where f is a trigonometric function. We therefore propose here a complete classification of invertible R-matrices solving DQYBE for V = Cn. We remind that we choose h to be the Cartan algebra of Gln(C) with basis vectors hi = e (n) ii ∈Mn(C) in the standard n× n matrix notation. This fixes in turn the normalization of the coordinate λ up to an overall multiplicator set so as to eliminate the prefactor 2γ. This classification is proposed within the following framework. i. We consider non-spectral parameter dependent R-matrices. They are generally called “constant” in the literature on quantum R-matrices but this denomination will never be used here in this sense since it may lead to ambiguities with respect to the presence in our matrices of “dynamical” parameters. This implies that a priori no elliptic dependence of the solutions in the dynamical variables is expected: at least in the Hecke case all dynamical elliptic quantum R-matrices are until now affine solutions. 4 J. Avan, B. Billaud and G. Rollet ii. We assume the matrix R to be invertible. Non-invertible R-matrices are expected to correspond to an inadequate choice of auxiliary space V (e.g. reducible). It precludes even the proof of commutation of the traces of monodromy matrices, at least by using the dynamical quantum group structure, hence such R-matrices present in our view a lesser interest. iii. We assume that the elements of the matrix R have sufficient regularity properties as functions of their dynamical variables, so that we are able to solve any equation of the form A(λ)B(λ) = 0 as A(λ) = 0 or B(λ) = 0 on the whole domain of variation Cn of λ except of course possible isolated singularities. In other words, we eliminate the possibility of “domain-wise zero” functions with no overlapping non-zero values. This may of course exclude potentially significant solutions but considerably simplifies the (already quite lengthy) discussion of solutions to DQYBE. iv. Finally we shall hereafter consider as “(pseudo)-constant” all functions of the variable λ with an integer periodicity, consistent with the chosen normalization of the basis (hi)i∈N∗n . Indeed such functions may not be distinguished from constants in the equations which we shall treat. After having given some preliminary results in Sections 2 and 3 presents key procedures allowing to define an underlying partition of the indices N∗n into r subsets together with an associated “reduced” ∆-incidence matrix MR ∈ Mr({0, 1}) derived from the ∆-incidence ma- trix M. The giving of this partition and the associated matrix MR essentially determines the general structure of the R-matrix in terms of constituting blocks. In Section 4, we shall establish the complete forms of all such blocks by solving system (S). The Hecke-type solutions will appear as a very particular solution1. Section 5 then presents the form of a general solution of DQYBE, and addresses the issue of the moduli structure of the set of solutions. The building blocks of any solution are in particular identified as weak Hecke type solutions or scaling thereof. The continuity of solutions in the moduli space are also studied in details. Finally we briefly conclude on the open problems and outlooks. 2 Preparatory material The following parametrization is adopted for the R-matrix R = n∑ i,j=1 ∆ije (n) ij ⊗ e (n) ji + n∑ i 6=j=1 dije (n) ii ⊗ e (n) jj . A key fact of our resolution is that since the R-matrix is assumed to be invertible, its determinant is non zero. Let n ≥ 2. Since the matrix R satisfies the zero weight-condition, for any i, j ∈ N∗n, the vector spaces Ce(n) ii ⊗ e (n) ii and Ce(n) ij ⊗ e (n) ji ⊕ Ce(n) ij ⊗ e (n) ji are stable. Then its determinant is given by the factorized form det(R) = n∏ i=1 ∆ii n∏ j=i+1 {dijdji −∆ij∆ji}. (det) This implies that all ∆ii are non-zero, and that ∆ij∆ji 6= 0, if dijdji = 0, and vice versa. Using this parametrization, we now obtain the equations obeyed by the coefficients of the R-matrix from projecting DQYBE on the basis (e (n) ij ⊗ e (n) kl ⊗ e (n) mp)i,j,k,l,m,p∈N∗n of n2 × n2 × n2 1For more details, see Subsection 5.5. Classification of Non-Affine Non-Hecke Dynamical R-Matrices 5 matrices. Only fifteen terms are left due to the zero-weight condition. Occurrence of a shift by 2γ (normalized to 1) of the i-th component of the dynamical vector λ will be denoted “(i)”. Distinct labels i, j and k mean distinct indices. The equations then read ∆ii∆ii(i){∆ii(i)−∆ii} = 0 (G0), dijdij(i){∆ii(j)−∆ii} = 0 (F1), djidji(i){∆ii(j)−∆ii} = 0 (F2), dij{∆ii(j)∆ij(i)−∆ii(j)∆ij −∆ji∆ij(i)} = 0 (F3), dji{∆ii(j)∆ij(i)−∆ii(j)∆ij −∆ji∆ij(i)} = 0 (F4), dij(i){∆ii∆ji(i)−∆ii∆ji + ∆ji∆ij(i)} = 0 (F5), dji(i){∆ii∆ji(i)−∆ii∆ji + ∆ji∆ij(i)} = 0 (F6), ∆2 ii(j)∆ij − (dijdji)∆ij(i)−∆ii(j)∆ 2 ij = 0 (F7), ∆2 ii∆ji(i)− (dijdji)(i)∆ij −∆ii∆ 2 ji(i) = 0 (F8), (S) ∆iidij(i)dji(i)−∆ii(j)dijdji + ∆ij(i)∆ji{∆ij(i)−∆ji} = 0 (F9), dij(k)djk(i)dik − dijdjkdik(j) = 0 (E1), djkdik(j){∆ij(k)−∆ij} = 0 (E2), dij(k)dik{∆jk(i)−∆jk} = 0 (E3), dij(k){∆ij(k)∆jk + ∆ji(k)∆ik −∆ik∆jk} = 0 (E4), djk{∆ij(k)∆jk + ∆ik(j)∆kj −∆ij(k)∆ik(j)} = 0 (E5), dij(k)dji(k)∆ik − djkdkj∆ik(j) + ∆ij(k)∆jk{∆ij(k)−∆jk} = 0 (E6). Treating together coefficients of e (n) ij ⊗ e (n) ji and e (n) ii ⊗ e (n) ii as ∆-coefficients is consistent since both tensor products may be understood as representing some universal objects e ⊗ e∗, components of a universal R-matrix R in some abstract algebraic setting. The d-coefficients of e (n) ii ⊗ e (n) jj are in this sense more representation dependent objects and we shall see indeed that they exhibit some gauge freedom in their explicit expression. More generally, in order to eliminate what may appear as spurious solutions we immedi- ately recall three easy ways of obtaining “new” solutions to DQYBE from previously obtained solutions. Let (αi)i∈N∗n be a family of functions of the variable λ. Define the dynamical diagonal operator F12 = eα1(h2)eα2 , where α is the λ-dependent vector α = n∑ i=1 αihi ∈ h. Proposition 2.1 (dynamical diagonal twist covariance). If the matrix R is a solution of DQYBE, then the twist action R′ = F12RF −1 21 is also a solution of DQYBE. Denoting βi = eα i , this is the origin of a particular, hereafter denoted, “twist-gauge” arbitra- riness on the d-coefficients, defined as2 dij → d′ij = βi(j) βi βj βj(i) dij , ∀ i, j ∈ N∗n. Proof. For any distinct labels a, b and c, the operator e±αc commutes with any operator with labels a and/or b and shifted in the space of index c, such as Rab(hc), e±αa(hc) or e±αa(hb+hc). Moreover, the zero-weight condition implies that e±αa(hb+hc) also commute with Rbc. By directly plugging R′ into the l.h.s. of DQYBE and using DQYBE for R, we can write R′12(h3)R′13R ′ 23(h1) = eα1(h2+h3)eα2(h3)R12(h3)e−α2(h1+h3)eα3R13e−α3(h1) 2For more details, see Propositions 4.2, 4.3 and 4.4. 6 J. Avan, B. Billaud and G. Rollet × e−α1eα2(h1+h3)eα3(h1)R23(h1)e−α3(h1+h2)e−α2(h1) = eα1(h2+h3)eα2(h3)eα3R12(h3)R13R23(h1)e−α3(h1+h2)e−α2(h1)e−α1 = eα1(h2+h3)eα2(h3)eα3R23R13(h2)R12e−α3(h1+h2)e−α2(h1)e−α1 = eα2(h3)eα3R23e−α3(h2)R′13(h2)eα1(h2)R12e−α2(h1)e−α1 = R′23R ′ 13(h2)R′12, where the equality e−αa(hb)eαa(hb) = 1 ⊗ 1 ⊗ 1 is used when needed. It is then immediate to check that R′ = n∑ i,j=1 ∆ije (n) ij ⊗ e (n) ji + n∑ i 6=j=1 d′ije (n) ii ⊗ e (n) jj , where the d-coefficients of R′ are given as in the proposition. � Corollary 2.1. Let (αij)i,j∈N∗n ∈ Cn2 be a family of constants, denoting βij = eα ij−αji, there exists a non-dynamical gauge arbitrariness on the d-coefficients as2 dij → d′ij = βijdij , ∀ i, j ∈ N∗n. Proof. Introducing the family (αi)i∈N∗n of functions of the variable λ, defined as αi = n∑ k=1 αikλk, for any i ∈ N∗n, it is straightforward to verify that βij = βi(j) βi βj βj(i) , for any i, j ∈ N∗n. � Remark 2.1. The dynamical twist operator F can be identified as the evaluation representation of a dynamical coboundary operator. Let Raa and Rbb be two R-matrices, solutions of DQYBE respectively represented on auxiliary spaces Va and Vb, being Etingof-modules of the underlying dynamical Abelian algebras ha and hb. Then Va ⊕ Vb is an Etingof-module for ha + hb. Let gab and gba be two non-zero constants, 1ab and 1 ba respectively the identity operator in the subspaces Va ⊗ Vb and Vb ⊗ Va. Define the new object Rab,ab = Raa + gab1 ab + gba1 ba +Rbb ∈ End((Va ⊕ Vb)⊗ (Va ⊕ Vb)), where the sum “+” should be understood as a sum of the canonical injections of each component operator into End((Va ⊕ Vb)⊗ (Va ⊕ Vb)). Proposition 2.2 (decoupled R-matrices). The matrix Rab,ab is an invertible solution of DQYBE represented on auxiliary space Va ⊕ Vb with underlying dynamical Abelian algebra ha + hb. Proof. Obvious by left-right projecting DQYBE onto the eight subspaces of (Va⊕Vb)⊗3 yielding a priori sixty-four equations. The new R-matrix is diagonal in these subspaces hence only eight equations survive. Among them, the only non-trivial equations are the DQYBE for Raa and Rbb, lying respec- tively in End(V ⊗3 a ) and End(V ⊗3 b ), up to the canonical injection into End((Va ⊕ Vb)⊗3), since Raa,bb depends only on coordinates in h∗a,b, and by definition of the canonical injection ha,b acts as the operator 0 on Vb,a. The six other equations are trivial because in addition they contain two factors 1 out of three. � Iterating m times this procedure will naturally produce R-matrices combining m “sub”-R- matrices, hereafter denoted “irreductible components”, with m(m − 1) identity matrices. To this end, it is not necessary to assume that the quantities gab and gba factorizing the identity Classification of Non-Affine Non-Hecke Dynamical R-Matrices 7 operators 1ab and 1ba linking the matrices Raa and Rbb should be constants3. Since, as above, the canonical injection hc acts as 0 on Vb,a, for any third distinct label c, it is sufficient to assume gab and gba to be non-zero 1-periodic functions in coordinates in h∗a and h∗b , the dependence on any coordinate in h∗c remaining free. Finally, a third construction of new solutions to system (S) from already known ones now stems from the form itself of (S). Let R be a matrix, solution of Gln(C)-DQYBE, with Cartan algebra h(n) having basis vectors h (n) i = e (n) ii , for any i ∈ N∗n, and I = {ia, a ∈ N∗m} ⊆ N∗n an ordered subset of m indices. We introduce the matrices eIij = e (m) σI(i)σI(j) ∈ Mm(C), for any i, j ∈ I, and define the bijection σI: I −→ N∗m as σI(ia) = a, for any a ∈ N∗m. Proposition 2.3 (contracted R-matrices). The contracted matrix RI = ∑ i,j,k,l∈I Rij,kle I ij ⊗ eIkl of the matrix R to the subset I is a solution of Glm(C)-DQYBE, with dynamical algebra h(m) having basis vectors h (m) a = e (m) aa , for any a ∈ N∗m. Proof. Obvious by direct examination of the indices structure of the set of equations (S). No sum over free indices occur, due to the zero-weight condition. Both lhs and rhs of all equations in (S) can therefore be consistently restricted to any subset of indices. � Remark 2.2. Formally the matrix eIij consists in the matrix e (n) ij , from which the lines and columns, whose label does not belong to the subset I, are removed. We shall completely solve system (S) within the four conditions specified above, all the while setting aside in the course of the discussions all forms of solutions corresponding to the three constructions explicited in Propositions 2.1, 2.2 and 2.3, and a last one explicited later in Propositions 4.2, 4.3 and 4.4. A key ingredient for this procedure will be the ∆-incidence matrix M∈Mn({0, 1}) of coefficients defined as mij = 0 if and only if ∆ij = 0. 3 The ∆-incidence matrix and equivalence classes We shall first of all consider several consistency conditions on the cancelation of d-coefficients and ∆-coefficients, which will then lead to the definition of the partition of indices indicated above. 3.1 d-indices Two properties are established. Proposition 3.1 (symmetry). Let i, j ∈ N∗n such that dij = 0. Then, dji = 0. Proof. If dij = 0, ∆ij∆ji 6= 0. From (G0) one gets ∆ii(i) = ∆ii and ∆jj(j) = ∆jj . From (F7), one gets ∆ij = ∆2 ii(j) ∆ii . This implies now that ∆ij(i) = ∆ij . (F4) then becomes dji∆ji∆ij(i) = 0, hence dji = 0. � Proposition 3.2 (transitivity). Let i, j, k ∈ N∗n such that dij = 0 and djk = 0. Then, dik = 0. Proof. From dij = 0 and (F9), one now gets ∆ij = ∆ji. From ∆ij(i) = ∆ij now follows that ∆ji(i) = ∆ji, hence ∆ij(j) = ∆ij . 3For a more indepth characterization of the quantities gab, see Proposition 4.5. 8 J. Avan, B. Billaud and G. Rollet From (F8), ∆ij = ∆ji = ∆ii = ∆jj = ∆, where the function ∆ is independent of variables λi and λj . Similarly, one also has ∆jk = ∆kj = ∆kk = ∆jj = ∆, independently of variables λk and λj . Writing now (E6) with indices jki and (E5) with indices jik yields dikdki = ∆ki{∆−∆ki} = ∆ki{∆−∆ik} and dik{∆−∆ik −∆ki} = 0. From which we deduce, if dik 6= 0, that ∆ = ∆ik + ∆ki. Then dikdki = ∆ki{∆−∆ki} = ∆ki∆ik, and det(R) = 0. Hence, one must have dik = 0. � Corollary 3.1. Adding the axiom iDi, for any i ∈ N∗n, the relation defined by iDj ⇔ dij = 0 is an equivalence relation on the set of indices N∗n. Remark 3.1. The D-class generated by any index i ∈ N∗n will be denoted I(i) = {j ∈ N∗n ∣∣ jDi}, and we will introduce the additional subset I0 = {i ∈ N∗n ∣∣ I(i) = {i}} of so-called “free” indices. For any subset I of the set of indices N∗n and any m ∈ N∗, let us also define the set I(m,�D ) = {(ia)a∈N∗m ∈ (N∗n)m | a 6= b ⇒ ia��D ib}. In the following, we will actually consider only the case m ∈ {2, 3}. An element of I(2,�D ) (resp. I(3,�D )) will be refereed as a ��D -pair (resp. ��D -triplet) of indices. 3.2 ∆-indices We establish a key property regarding the propagation of the vanishing of ∆-coefficients. Proposition 3.3. Let i, j ∈ N∗n such that ∆ij = 0. Then, ∆ik∆kj = 0, for any k ∈ N∗n. Proposition 3.4 (contraposition). Let i, j ∈ N∗n. Equivalently, if there exist k ∈ N∗n such that ∆ik∆kj 6= 0, then ∆ij 6= 0. Proof. If ∆ij = 0 then dijdji 6= 0. It follows from Proposition 3.2 that dik 6= 0 or dkj 6= 0, for all k 6= i, j. Assume that dik 6= 0 hence dki 6= 0. (E4) with indices ikj reads dik(j)[∆ik(i)∆kj + ∆ij{∆ki(j)−∆kj}] = 0, hence ∆ik = 0 or ∆kj = 0. If instead dkj 6= 0 hence djk 6= 0. (E5) with indices ijk directly yields ∆ik(j)∆kj = 0 with the same conclusion. � Proposition 3.5. The relation defined by i∆j ⇔ ∆ij∆ji 6= 0 is an equivalence relation on the set of indices N∗n. Moreover, any D-class is included in a single ∆-class. Classification of Non-Affine Non-Hecke Dynamical R-Matrices 9 Proof. Reflexivity and symmetry are obvious. Transitivity follows immediately from Proposi- tion 3.4. If i∆k ⇔ ∆ik∆ki 6= 0 and k∆j ⇔ ∆kj∆jk 6= 0, hence ∆ik∆kj 6= 0 and ∆jk∆ki 6= 0. Then, ∆ij 6= 0 and ∆ji 6= 0, i.e. i∆j. The second part of the proposition follows immediately from (det). � Corollary 3.2. Denote {Jp, p ∈ N∗r} the set of r ∆-classes, which partitions the set of in- dices N∗n. For any p ∈ N∗r, there exist lp ∈ N, a so-called “free” subset I(p)0 = Jp ∩ I0 of free indices (possibly empty), and lp D-classes generated by non-free indices (possibly none), denoted I(p)l with l ∈ N∗lp, such that Jp = lp⋃ l=0 I(p)l is a partition. Finally, iDj, if and only if ∃ l ∈ N∗lp | i, j ∈ I(p)l . 3.3 (Reduced) ∆-incidence matrix The ∆-incidence matrix M = n∑ i,j=1 mije (n) ij ∈Mn({0, 1}) is defined as follows mij = 1 ⇔ ∆ij 6= 0 and mij = 0 ⇔ ∆ij = 0. Let us now use the ∆-class partition and Propositions 3.3, 3.5 and 3.5 to better characterize the form of the ∆-incidence matrixM of a solution of DQYBE. The key object here will be the so-called reduced ∆-incidence matrix MR. Proposition 3.6. Let I, J two distinct ∆-classes such that ∃ (i, j) ∈ I× J ∣∣ ∆ij 6= 0. Then, for any pair of indices (i, j) ∈ I× J, ∆ij 6= 0. Proof. Let i′ ∈ I and j′ ∈ J. Applying Proposition 3.4 to ∆i′i∆ij 6= 0, we deduce ∆i′j 6= 0. Then ∆i′j′ 6= 0, since ∆i′j∆jj′ 6= 0. � Remark 3.2. In the proof of this proposition, note here that nothing forbids i′ = i and/or j′ = j. To facilitate their writing and reading, this convention will be also used in Proposition 4.1, Lemmas 4.5 and 4.6, as well as Theorems 4.2, 4.3 and 4.4. Corollary 3.3. Let I, J two distinct ∆-classes. Then either all connecting ∆-coefficients in ∆ij, with (i, j) ∈ I× J}, are zero or all are non-zero. This justifies that the property of vanishing of ∆-coefficients shall be from now on denoted with overall ∆-class indices as ∆IJ = 0 or ∆IJ 6= 0. This leads now to introduce a reduced ∆-incidence matrix MR = r∑ p,p′=1 mR pp′e (r) pp′ ∈Mr({0, 1}), defined as mR pp′ = 1 ⇔ ∆JpJp′ 6= 0 and mR pp′ = 0 ⇔ ∆JpJp′ = 0. Proposition 3.7. The relation defined by I � J ⇔ ∆IJ 6= 0 is a partial order on the set of ∆-classes. Proof. If I � J and J � I, then ∆IJ 6= 0 and ∆JI 6= 0. Hence, for all (i, j) ∈ I× J, ∆ij∆ji 6= 0, i.e. I = J. If I � J and J � K, then ∆IJ 6= 0 and ∆JK 6= 0. Hence, from Proposition 3.4, ∆ik 6= 0, for all (i, k) ∈ I×K, i.e. I � K. � 10 J. Avan, B. Billaud and G. Rollet Two ∆-classes I and J shall be refereed hereafter as “comparable”, if and only if I � J or J � I, which will be denoted I ≺� J. This order on ∆-classes is of course not total, because there may exist ∆-classes which are not comparable, i.e. such that ∆IJ = ∆JI = 0, being denoted I��≺� J. The order � is to be used to give a canonical form to the matrix MR in two steps, and more particularly the strict order � deduced from � by restriction to distinct ∆-classes. Unless otherwise stated, in the following, the subsets I, J and K are three distinct ∆-classes. Proposition 3.8 (triangularity). The reduced ∆-incidence matrix MR is triangularisable in Mr({0, 1}). Proof. The strict order � defines a natural oriented graph on the set of ∆-classes. Triangularity property of the order implies that no cycle exists in this graph. To any ∆-class I one can then associate all linear subgraphs ending on I as Jp1 � Jp2 � · · · � Jpk � I. There exist only a finite number of such graphs (possibly none) due to the non-cyclicity property. One can thus associate to the ∆-class I the largest value of k introduced above, denoted by k(I). We now label ∆-classes according to increasing values of k(I), with the additional convention that ∆-classes of same value of k(I) are labeled successively and arbitrarily. The labels are denoted as l(I) ∈ N∗r in increasing value, and we have the crucial following lemma. Lemma 3.1. If l(I) < l(J), then ∆JI = 0. Proof. By contraposition, if ∆JI 6= 0 and I 6= J, then J � I. Hence k(I) ≥ k(J) + 1 > k(J), which is impossible if l(I) < l(J), by definition of the labeling by increasing values of k(I). � Let us now introduce the permutation σ: p 7−→ l(Jp) ∈ Sr, its associated permutation matrix Pσ = r∑ p=1 e (r) σ(p)p ∈ Glr({0, 1}), of inverse P−1 σ = Pσ−1 , and the permuted reduced ∆-incidence matrix Mσ R = PσMRP −1 σ . It is straightforward to check that mR,σ pp′ = mR σ−1(p)σ−1(p′). From Lemma 3.1, we deduce that, if p = σ(q) < σ(q′) = p′, then ∆JqJq′ = 0, and mR,σ pp′ = mR qq′ = 0, i.e. the matrix Mσ R is upper-triangular. � Corollary 3.4. Denoting Jσp = Jσ(p), if p < p′, then either Jσp��≺� Jσp′ or Jσp � Jσp′. The characterization of a canonical form for the matrix MR can now be further precise. Proposition 3.9. If I��≺� J and I � K, then J��≺�K. Proof. By assumption, remark that ∆IJ = ∆JI = ∆KI = 0 and ∆IK 6= 0. Let (i, j, k) ∈ I×J×K. Since ∆ij = 0 and ∆ki = 0, from (det), dijdji 6= 0 and dikdki 6= 0, but ∆ik 6= 0. When written with indices ijk, (E4) reduces to dij(k)∆ik∆jk = 0, hence ∆jk = 0. When written with indices ikj, (E4) reduces to dik(j)∆kj∆ik(j) = 0, hence ∆kj = 0. � Proposition 3.10. If I��≺� J and K � I, then J��≺�K. Proof. Identical to the proof of Proposition 3.9 using (E5), written with indices jik and ikj. � Proposition 3.11. If I��≺� J (resp. I ≺� J) and I ≺� K, then J��≺�K (resp. J ≺� K). Proof. From Propositions 3.9 and 3.10, if I��≺� J and I ≺� K, then J��≺�K. If I ≺� J, and assuming that J��≺�K, there is a contradiction with I ≺� K, then J ≺� K. � Corollary 3.5. Let p < p′, p′′. Hence, i. if Jσp��≺� Jσp′ and Jσp � Jσp′′, then Jσp��≺� Jσp′′; ii. if Jσp � Jσp′ and Jσp � Jσp′′, with p′ < p′′, then Jσp′ � Jσp′′. Classification of Non-Affine Non-Hecke Dynamical R-Matrices 11 Proposition 3.12 (block upper-triangularity). The reduced ∆-incidence matrix MR is similar to a block upper-triangular matrix in Mr({0, 1}). That is: there exists a permutation π ∈ Sr and a partition of the set N∗r in s subsets Pq = {pq + 1, . . . , pq+1} (with the convention that p1 = 0 and ps+1 = r), of respective cardinality rq, such that mR,π pp′ = 1 ⇔ ∃ q ∈ N∗s ∣∣ (p, p′) ∈ P(2,<) q . i.e. the matrix Mπ R = r∑ p,p′=1 mR,π pp′ e (r) pp′ is graphically represented by blocks as Mπ R =  Tr1 Or1r2 Or1rs Or2r1 Tr2 Trs−1 Ors−1rs Orsr1 Orsrs−1 Trs  ∈Mr({0, 1}), where the type T , O block matrices are def ined by Tr′ =  1 1 0 0 0 1 ∈Mr′({0, 1}) and Or′r′′ =  0 0 0 0 ∈Mr′,r′′({0}), Or′ = Or′r′∈Mr′({0}). Remark 3.3. For any set I of integers and any m ∈ N∗, by analogy with the definition of the set I(m,�D ), we adopt the notations I(m,<) = {(ia)a∈N∗m ∈ Im | a < b⇒ ia < ib} and I(m,�D ,<) = {(ia)a∈N∗m ∈ I(m,�D ) | a < b⇒ ia < ib}. For example, a pair of labels q, q′ ∈ N∗s such that q < q′ (resp. a ��D -pair of indices (i, j) ∈ I2 such that i < j) belongs to the set N∗(2,<) s (resp. I(2,�D ,<)). Proof. The proof relies on a recursion procedure on the value of the size r of the matrix MR. The proposition being trivial for r ∈ {1, 2}, let us assume that r ≥ 3. 1. Re-ordering from line 1. Starting from the matrixMσ R, whose existence is guaranteed by Proposition 3.8, its upper-triangularity is used following Corollary 3.5. Remember that label ordering and class-ordering run contrary to each other. Note p (1) 1 = 1. Since Jσ1 is always comparable to itself, the set of ∆-classes comparable to Jσ1 is not empty, and we will denote r1 ∈ N∗r its cardinality. If r1 = 1, i.e. if Jσ1 is not comparable to any other ∆-class, line 1 of matrixMσ R consists of an one-label block mR,σ 11 = 1, and the process stops. Assuming that r1 ∈ {2, . . . , r}, consider the subset {p ∈ {2, . . . , r} ∣∣ Jσ1 � Jσp} 6= ∅. This set is naturally totally ordered. Let us then denote its elements as p (1) q by increasing value, where q ∈ {2, . . . , r1}. Then, by convention, (p (1) q , p (1) q′ ) ∈ N∗(2,<) p (1) r1 , if and only if (q, q′) ∈ N∗(2,<) r1 . Moreover, by construction, we have that Jσ1 � Jσ p (1) q and Jσ1 � Jσ p (1) q′ , for any q, q′ ∈ {2, . . . , r1}. Then, from Corollary 3.5, Jσ p (1) q � Jσ p (1) q′ , if and only if (q, q′) ∈ N∗(2,<) r1 , i.e. Jσ1 � Jσ p (1) 2 � · · · � Jσ p (1) r1 , 12 J. Avan, B. Billaud and G. Rollet where no sign � can be reversed. In particular, the ∆-classes in { Jσ p (1) q , q ∈ {2, . . . , r1} } are comparable one-to-one. Since Jσ1 is only comparable to the ∆-classes in { Jσ p (1) q , q ∈ {2, . . . , r1} } , no other ∆-class is comparable to any ∆-class Jσ p (1) q , with q ∈ {2, . . . , r1}. This implies that mR,σ p (1) q p (1) q′ = 1 ⇔ (q, q′) ∈ N∗(2,<) r1 , and that mR,σ p (1) q p = mR,σ pp (1) q = 0, ∀ q ∈ N∗r1 and ∀ p ∈ N∗r \ {p(1) q , q ∈ N∗r1}. Let π1 ∈ Sr be the unique permutation such that π1(p(1) q ) = q, ∀ q ∈ N∗r1 , and that π1 is increasing on N∗r \ {p (1) q , q ∈ N∗r1}. We apply the same reasoning as in the end of the proof of Proposition 3.8. The coefficients of the permuted matrixMπ1◦σ R now satisfy the following equalities mR,π1◦σ pp′ = 1 ⇔ (p, p′) ∈ N∗(2,<) r1 , and mR,π1◦σ pp′ = mR,π1◦σ p′p = 0, ∀ p ∈ N∗r1 and ∀ p′ ∈ {r1 + 1, . . . , r}. Furthermore, the increasing property of the permutation π1 transfers the upper-triangularity of the matrixMσ R to the matrixMπ1◦σ R , which can finally be graphically represented by blocks as Mπ1◦σ R =  Tr1 Or1r′ Or′r1 M′R . 2. Recursion on r. Let assume that the statement is true for any reduced ∆-incidence matrix of size r′ ∈ N∗r−1 associated with a solution of DQYBE. Using the previously defined re-ordering procedure on the first line of a matrix Mσ R ∈ Mr({0, 1}), there exists a upper- triangular reduced ∆-incidence matrix M′R ∈ Mr′({0, 1}) of size r′ = r − r1 < r, which is moreover associated with a solution of DQYBE from Proposition 2.3. The recursion hypothesis can now be applied to the first line of the matrix M′R describing the order of the r′ remaining ∆-classes. 3. Recursive construction of π and {Pq, q ∈ N∗s}. Since the number of ∆-classes is finite, the process described above comes to an end after a finite number s ∈ N∗r of iterations. The qth iteration insures the existence of an integer rq ∈ N∗s and a permutation πq ∈ Sr, built by recursion. Defining pm = m−1∑ m′=1 rm′ ∈ N∗r , for any m ∈ N∗q , the integer rq is the cardinality of the set of ∆-classes comparable to the ∆-class Jπp−1◦···◦π1◦σ pq−1+1 , being the first remaining ∆-class after q − 1 iterations. Introducing the totally ordered set {p(q) q′ , q ′ ∈ N∗rq} of indices of such ∆-classes and putting pq = pq−1 + rq, the permutation πq re-orders the indices as follows πq(p) = p, ∀ p ∈ N∗pq and πq(p (q) q′ ) = pq−1 + q′, ∀ q′ ∈ N∗rq , πq being increasing on {pq + 1, . . . , r} \ {p(q) q′ , q ′ ∈ N∗rq}. Finally, the permutation π = πs ◦ · · · ◦ π1 ◦ σ ∈ Sr leads to the expected permuted matrix Mπ R, and the partition N∗r = s⋃ q=1 Pq stands by construction. � Classification of Non-Affine Non-Hecke Dynamical R-Matrices 13 3.4 Classif ication For convenience, we will now identify in the following the reduced ∆-incidence matrix MR and its associated block upper-triangular matrix Mπ R, as well as the ∆-classes in {Jp, p ∈ N∗r} and the re-ordered ∆-classes in {Jπp , p ∈ N∗r}. Let us conclude this section by fully describing the ∆- incidence structure of a general R-matrix to complete the classification of solutions of DQYBE, together with the required steps to end the resolution of system (S). Theorem 3.1 (∆-incidence matrices). Let n ≥ 2. Then, any R-matrix, solution of DQYBE, is characterized in particular by • an ordered partition of the indices N∗n into r ∆-classes Jp of respective cardinality np, • an ordered partition of the indices N∗r into s subsets Pq = {pq + 1, . . . , pq+1}, of respective cardinality rq (with the convention that p1 = 0 and ps+1 = r), • an ordered partition of each ∆-class Jp into a “free” subset I(p)0 = Jp ∩ I0 (possibly empty) of cardinality n (p) 0 , and lp D-classes I(p)l generated by non-free indices (possibly none) of respective cardinality n (p) l ; such that the following union is an ordered partition of the set of indices N∗n N∗n = s⋃ q=1 Kq = r⋃ p=1 Jp = r⋃ p=1 lp⋃ l=0 I(p)l , denoting Kq = ⋃ p∈Pq Jp the set of ∆-classes, of cardinality Nq = ∑ p∈Pq np ∈ N∗n, associated to each subset Pq. Re-expanding its reduced ∆-incidence matrix MR, the R-matrix has a ∆-incidence matrix M = ∑ p,p′∈N∗r p≤p′ mR pp′e (r) pp′ ⊗ Enpnp′ , which can be graphically represented as M =  T (1) O(1,2) O(1,s) O(2,1) T (2) T (s−1) O(s−1,s) O(s,1) O(s,s−1) T (s)  ∈Mn({0, 1}), where the matrices T (q) = ∑ p,p′∈Pq p≤p′ e (rq) pp′ ⊗ Enpnp′ are graphically represented as T (q) =  Enpq+1 Enpq+1npq+1 Onpq+2npq+1 Onpq+1npq+1 Onpq+1−1npq+1 Enpq+1 ∈MNq({0, 1}), with O(q,q′) = ONqNq′ , and where the type E matrices are defined like the type O matrices except that 0 is replaced by 1, i.e. Er′r′′ ∈Mr′r′′({1}), Er′ = Er′r′ ∈Mr′({1}). 14 J. Avan, B. Billaud and G. Rollet Remark 3.4. By ordered partition, we mean that the indices appear in the partition in the canonical order of integers. For example, a ��D-pair of indices (i, j) ∈ Kq×Kq′ with (q, q′) ∈ N∗(2,<) s satisfies by construction i < j, i.e. is an ordered pair (i, j) ∈ N∗(2,�D ,<) n . Proof. Theorem 3.1 is almost entirely a direct consequence of Proposition 3.12, the only un- proved point being the re-ordering of each subset Jp = lp⋃ l=0 I(p)l as an ordered partition, for any p ∈ Pq, with q ∈ N∗s. To this end, for any l ∈ Nlp , we first denote the elements of the subset I(p)l , when not empty, by increasing values as p (l) i , with i ∈ N∗ n (p) l . If lp = 0, i.e. if n (p) 0 ≥ 1, then Jp = I(p)0 is a single free subset, and is already ordered. If n (p) 0 ≥ 1 and lp ≥ 1, we define the permutation σp ∈ Sn, whose support is a subset of Jp, as σp(p (0) i ) = i+ pq, ∀ i ∈ N∗ n (p) 0 and σp(p (l) i ) = i+ l−1∑ l′=0 n (p) l′ + pq, ∀ i ∈ N∗ n (p) l and ∀ l ∈ N∗lp . If n (p) 0 = 0, i.e. if lp ≥ 1, then Jp = lp⋃ l=1 I(p)l does not contain free indices, and we define the permutation σp ∈ Sn just as above, but omitting the first part of this definition. Therefore, the ∆-class Jp can be written as the following ordered partition Jσp = σp(I (p) 0 ) ∪ σq(Jp \ I(p)0 ) = lp⋃ l=0 σp ( I(p)l ) , where the exponent “σ” indicates that the permutation σp is applied. Moreover, since the supports of the permutations {σp}p∈N∗r are disjoint and since the set {Jp, p ∈ N∗r} has a natural order from Proposition 3.12, the permutation σ = σ1 ◦ · · · ◦ σr ∈ Sn re-orders as expected each element of the set of indices N∗n, i.e. (N∗n)σ = r⋃ p=1 Jσp is an ordered partition. Finally, for convenience, as earlier, we drop the exponent “σ”, and identify the subsets {Jp, p ∈ N∗r} and the set of indices N∗n with the re-ordered ones {Jσp , p ∈ N∗r} and (N∗n)σ. � Corollary 3.6. In addition with the family of diagonal elements (∆ii)i∈N∗n, the associated non- zero R-matrix elements to be determined are the coefficients i) ∆ij for all pairs of indices (i, j), i and j belonging to the same D-class I(p)l . ii) ∆ij, ∆ji, dij and dji, for all ��D -pairs of indices (i, j), i and j belonging to the same ∆-class Jp. This covers the cases of indices i and j • both in the free subset I(p)0 6= ∅, accordingly the corresponding contracted R-matrix will be refereed as “full” since all zero-weight elements are a priori non-zero; • in the free subset I(p)0 6= ∅ and in a D-class I(p)l ; • in two distinct D-classes I(p)l and I(p)l′ , with l < l′. iii) ∆ij, dij and dji, for all ��D -pairs of indices (i, j), i and j belonging to two distinct ∆- classes Jp and Jp′ of the same subset Kq, i.e. with (p, p′) ∈ P(2,<) q . This covers the cases of indices i and j Classification of Non-Affine Non-Hecke Dynamical R-Matrices 15 • in the two free subsets I(p)0 6= ∅ and I(p ′) 0 6= ∅; • in the free subset I(p)0 6= ∅ and in a D-class I(p ′) l , as well as the non-equivalent symmetric case of any pair of indices in a D-class I(p)l and the free subset I(p ′) 0 6= ∅; • in two D-classes I(p)l and I(p ′) l′ . iv) dij and dji, for all ��D -pairs of indices (i, j), i and j belonging to two distinct subsets Kq and Kq′, with (q, q′) ∈ N∗(2,<) s . Proof. This is a simple study of cases, when the indices i and j belong respectively to any possible subsets I(p)l and I(p ′) l′ , with p, p′ ∈ N∗r | p ≤ p′. Cases iii and iv are respectively reduced to (p, p′) ∈ P(2,<) q and (q, q′) ∈ N∗(2,<) s thanks to Proposition 3.12 on the upper-triangularity of the reduced ∆-incidence matrix MR. � 4 Resolution The resolution of the system (S) needs the introduction of the functions “sum” Sij and “deter- minant” Σij defined, for any pair of indices (i, j) ∈ (N∗n)2, as Sij = Sji = ∆ij + ∆ji and Σij = Σji = ∣∣∣∣ dij ∆ij ∆ji dji ∣∣∣∣ 6= 0. 4.1 Preliminaries We will begin these first considerations on the resolution of system (S) by solving cases i and iv of Corollary 3.6. To this end, let q ∈ N∗s and consider the subset Kq. Proposition 4.1 (inside a D-class). Let i ∈ Kq \ (Kq ∩ I0). Then, there exists a non-zero constant ∆I(i) such that the solution of system (S) restricted to the subset I(i) is given by ∆jj′ = ∆I(i), ∀ j, j′ ∈ I(i). Proof. This corresponds to the case i of Corollary 3.6. From the proof of Proposition 3.1, there exists a function ∆I(i) independent of the variable λk, for any k ∈ I(i), such that ∆jj = ∆jj′ = ∆j′j = ∆I(i), for any j, j′ ∈ I(i), with j 6= j′. It remains to prove that ∆I(i) is a 1-periodic function in the variable λk, for any k ∈ N∗n \ I(i). If n = 2, the proof of the proposition ends here. Assuming that n ≥ 3, it is possible to suppose without loss of generality that I(i) ( N∗n, the case of equality having been already treated. Let then k ∈ N∗n \ I(i). By construction, (j, k) ∈ N∗(2,�D ) n and (j′, k) ∈ N∗(2,�D ) n , then (E2) with indices jj′k implies that ∆I(i)(k) = ∆jj′(k) = ∆jj′ = ∆I(i). Reciprocally, it is straightforward to check that this is indeed a solution of system (S) re- stricted to the D-class I(i). The set of solutions of system (S) restricted to a D-class I(i) is exactly parametrized by the constant ∆I(i). � If the subset Kq is reduced to a single D-class, the resolution ends here. We must now consider that the subset Kq is not reduced to a single D-class. In particular, there exists a ��D -pair of indices (i, j) ∈ K(2,�D ,<) q . This also suggests to extend the notation ∆I(i) to ∆ii, even if i ∈ I0. 16 J. Avan, B. Billaud and G. Rollet Corollary 4.1. For any pair of indices (i, j) ∈ K(2,�D ) q satisfying one of the cases ii)–iv) of Corollary 3.6, (F1)–(F9) of system (S) is equivalent to (G0), (F1)⇔ (F2) ∆I(i)(i) = ∆I(i)(j) = ∆I(i) 6= 0 (G′0), (F3)⇔ (F4) ∆I(i){∆ij(i)−∆ij} −∆ji∆ij(i) = 0 (F ′1), (F5)⇔ (F6) ∆I(i){∆ji(i)−∆ji}+ ∆ji∆ij(i) = 0 (F ′2), (F7) ∆I(i)∆ij{∆ii −∆ij} − dijdji∆ij(i) = 0 (F ′3), (S′) (F8) ∆I(i)∆ji(i){∆ii −∆ji(i)} − dij(i)dji(i)∆ij = 0 (F ′4), (F9) ∆I(i){dij(i)dji(i)− dijdji}+ ∆ij(i)∆ji{∆ij(i)−∆ji} = 0 (F ′5). For later purpose, we will now introduce several lemmas, which restrain a priori the depen- dences of the ∆-coefficients (Lemmas 4.1 and 4.2) and d-coefficients (Lemmas 4.3 and 4.4) on the variable λk, for any k ∈ N∗n, as well as their acceptable form. Lemma 4.1. Let i ∈ Kq ∩ I0. Then, ∆I(i) 6= 0 is a constant. Proof. Since i ∈ I0, for any j ∈ N∗n \{i}, system (S′) applies to the pair of indices (i, j) ∈ (N∗n)2 (cf. Corollary 4.1). (G′0) then implies that ∆I(i) is constant. � Lemma 4.2. Let (i, j) ∈ K(2,�D ,<) q be a ��D -pair of indices. Then, ∆ij and ∆ji are 1-periodic functions in the variable λk, for any k ∈ N∗n \ (I(i) ∪ I(j)). Proof. If n = 2, i.e. N∗n \ (I(i) ∪ I(j)) = ∅, the lemma is empty. Assuming that n ≥ 3, let k ∈ N∗n \ (I(i) ∪ I(j)). Since (i, j, k) ∈ N∗(3,�D ) n , (E2) with indices ijk and with indices jik implies that ∆ij(k) = ∆ij and ∆ji(k) = ∆ji. � 4.2 Decoupling procedure This section is dedicated to the specific characterization of the decoupled R-matrices, as defined in Proposition 2.3, the main result being that any R-matrix, solution of DQYBE, characterized by a block-upper triangular matrix reduced ∆-incidence matrixMR with two or more triangular blocks, is in fact decoupled, up to a particular transformation explicited in the following. For the moment, let us focus on two fundamental lemmas, which describe the form of the non-zero d-coefficients. Lemma 4.3. Let s ≥ 2. Then, for any (q, q′) ∈ N∗(2,<) s , there exists a non-zero constant Σqq′ and a family of non-zero functions (gij)(i,j)∈Kq×Kq′ (with the property that gijgji = 1), such that, for any (i, j) ∈ Kq ×Kq′ Σij = Σqq′ , dij = √ Σqq′gij and dji = √ Σqq′gji. Proof. Let (i, j) ∈ Kq×Kq′ , from Corollary 4.1, system (S′) stands for the pair of indices (i, j). Hence, from (F ′5), we get that Σij(i) = dij(i)dji(i) = dijdji = Σij(= Σji(j) by symmetry), i.e. the function Σij is 1-periodic in variables λi and λj . If n = 2, or if n ≥ 2 and Nq = Nq′ = 1, the pair (i, j) is the only such pair of indices to consider. Assuming that Nq ≥ 2 and Nq′ ≥ 1, let k ∈ Kq \ {i}. It follows that ∆ij = ∆ji = ∆kj = ∆jk = 0 and i��D j��D k. From Proposition 4.1 (when kDi) or Lemma 4.2 (when k��D i), the function ∆ik is a 1-periodic function in the variable λj . This implies, when used in (E6), Σij(k) = Σkj = Σkj(k) = Σij , 4 which is the expected result if Nq′ = 1, as far as the determinant 4Theorem 4.1 enunciates a similar result to this part of the reasoning for any pair of indices (i, j) ∈ K(2, �D ,<) q . Classification of Non-Affine Non-Hecke Dynamical R-Matrices 17 is concerned. If Nq = 1 and Nq′ ≥ 2, the symmetrical result is obtained by exchanging the indices i and j, as well as the labels q and q′. Assuming that Nq ≥ 2 and Nq′ ≥ 2, both previous results apply, so that the function Σij does not depend on the index i ∈ Kq nor on the index j ∈ Kq′ , but only on the subsets Kq and Kq′ . There exists then a non-zero function, denoted Σqq′ by language abuse, which is 1-periodic in the variable λk, for any k ∈ Kq ∪Kq′ , such that Σij = Σqq′ , for any (i, j) ∈ Kq ×Kq′ . If s = 2, i.e. N∗n = Kq ∪Kq′ , the subsets Kq and Kq′ are the only subsets to consider. Assuming that s ≥ 3, let q′′ ∈ N∗s \ {q, q′} and k ∈ Kq′′ . Since (i, j, k) ∈ N∗(3,�D ) n , (E1) is non-trivial when written with indices kji and kji. This yields that dki(j)dij(k)dkjdji(k) = dkidijdkj(i)dji(k) = dijdkjdjidki(j) ⇒ dij(k)dji(k) = dijdji, implying that Σqq′ is also 1-periodic in the variable λk, for any k ∈ N∗n \ (Kq ∪ Kq′), hence is constant. There exist then two non-zero functions gij and gji, such that the functions dij =√ Σqq′gij and dji = √ Σqq′gji are the general solution of this equation, with the condition gijgji = 1. � Lemma 4.4. Let (i, j) ∈ K(2,�D ,<) q be a ��D -pair of indices. Then, there exist two non-zero functions d0 ij and d0 ji of the variable λk, for any k ∈ I(i)∪ I(j), and 1-periodic in other variable, and two non-zero functions gij and gji (with the property that gijgji = 1), such that dij = gijd 0 ij and dji = gjid 0 ji. Proof. From Corollary 4.1, system (S′) stands for any pair of indices (i, j) ∈ K(2,�D ,<) q . (F ′1) + (F ′2) implies that Sij(i) = Sij(=Sij(j) by symmetry), i.e. the function Sij is 1-periodic in va- riables λi and λj . When inserted in (F ′2), this yields ∆ji∆ij(i) = −∆I(i){∆ji(i)−∆ji} = ∆I(i){∆ij(i)−∆ij} ⇔ ∆I(i)∆ij = ∆ij(i){∆I(i) −∆ji}. Hence, from (F ′3), and since ∆ij 6= 0, we deduce that ∆ij(i)[{∆I(i) −∆ij}{∆I(i) −∆ji} − dijdji] = 0 ⇔ Σij = ∆I(i){∆I(i) − Sij}, which implies that the function Σij is also 1-periodic in variables λi and λj , because ∆I(i) is constant from Proposition 4.1 (when i ∈ Kq \ (Kq ∩ I0)) or Lemma 4.1 (when i ∈ Kq ∩ I0). Moreover, by exchanging the indices i and j, we get that Σij = ∆I(i){∆I(i) − Sij} = ∆I(j){∆I(j) − Sij}, ∀ (i, j) ∈ K(2,�D ,<) q . (4.1) Assuming that N∗n 6= I(i)∪I(j), let k ∈ N∗n\(I(i)∪I(j)). Otherwise, the lemma is trivial. From Proposition 4.1 and Lemma 4.2, which express that the functions ∆ij and ∆ji are 1-periodic in the variable λk, (4.1) implies that the functions Sij and Σij are 1-periodic in the variables λi, λj and λk. Then, the function d0 ij = Bij −∆ij is a particular solution of dijdji = Σij + ∆ij∆ji, being 1-periodic in the variable λk, for any k ∈ N∗n \ (I(i) ∪ I(j)), where the quantity Bij = Sij+ √ S2 ij+4Σij 2 is a root of the polynomial Pij(X) = X2 − SijX −Σij . Hence there exists a non- zero function gij such that dij = gijd 0 ij is the general solution of this equation, with the condition gijgji = 1. � 18 J. Avan, B. Billaud and G. Rollet Remark 4.1. Theorems 4.2, 4.3 and 4.4 will provide explicit expressions for the non-zero functions d0 ij and d0 ji, with (i, j) ∈ K(2,�D ,<) q and q ∈ N∗s, which will appear as the multiplicative invariant part of the d-coefficients. For (i, j) ∈ K(2,�D ,<) q , with q ∈ N∗s, there exists a second realization of the functions d0 ij and d0 ji, which also determine the functions gji, given by d0′ ij = d0′ ji = √ d0 ijd 0 ji and g′ij = gijg 0 ij , with g0 ij = √ d0 ij d0 ji , ∀ (i, j) ∈ K(2,�D ,<) q . Formally, this is the parametrization used in Lemma 4.3. In particular, both previous realizations of the functions g0 ij are also 1-periodic in the variable λk, for any k ∈ N∗n \ (I(i) ∪ I(j)). Moreover, let us point out that, extending the notation d0 ij to any pair of indices (i, j) ∈ Kq × Kq′ , with (q, q′) ∈ N∗(2,<) s , i.e. if we set d0 ij = √ Σqq′ in this case, the family of non-zero functions (d0 ij) (i,j)∈N∗(2,�D) n introduced by Lemmas 4.3 and 4.4 trivially satisfies (E1). Since we have introduced all the needed tools, we can now separately study (E1) in details. This particular treatment is justified by the fact that this equation, which is the only equation where three d-coefficients appear, is decoupled from other equations of system (S). It only constrains the functions gij , with (i, j) ∈ N∗(2,�D ,<) n . To this end, it is assumed that n ≥ 3, let i, j, k ∈ N∗n. Since d-coefficients are concerned, it is possible to consider that the triplet (i, j, k) is a ��D -triplet, i.e. (i, j, k) ∈ N∗(3,�D ) n . Otherwise, there exists {i′, j′} ⊆ {i, j, k} such that i′Dj′, and (E1) becomes trivial. We first establish that DQYBE shows another type of covariance, of which the twist covarian- ce is an example (cf. Proposition 2.1). This new symmetry of DQYBE is of great importance for characterizing the decoupled R-matrices. Let us now give the following definitions by analogy with [13]. Definition 4.1 (multiplicative 2-forms). Let I be a subset of the set of indices N∗n. i. A family of non-zero functions (αij)(i,j)∈I(2,�D) (resp. (αij)(i,j)∈I2) of the variable λ, such that αijαji = 1, for any (i, j) ∈ I(2,�D ) (resp. (i, j) ∈ I2), is called a ��D -multiplicative 2-form (resp. multiplicative 2-form). ii. A ��D -multiplicative 2-form (resp. multiplicative 2-form) (αij)(i,j)∈I(2,�D) (resp. (αij)(i,j)∈I2) is said to be ��D -closed (resp. closed), if it satisfies the cyclic relation αij(k) αij αjk(i) αjk αki(j) αki = 1, ∀ (i, j, k) ∈ I(3,�D ) (resp. ∀ (i, j, k) ∈ I3). iii. A ��D -multiplicative 2-form (resp. multiplicative 2-form) (αij)(i,j)∈I(2,�D) (resp. (αij)(i,j)∈I2) is said to be ��D -exact (resp. exact), if there exists a family of non-zero functions (αi)i∈I of the variable λ, such that αij = αi(j) αi αj αj(i) , ∀ (i, j) ∈ I(2,�D ) (resp. ∀ (i, j) ∈ I2). Classification of Non-Affine Non-Hecke Dynamical R-Matrices 19 Proposition 4.2. Let (αij)(i,j)∈I(2,�D) (resp. (αij)(i,j)∈I2) be a ��D -closed ��D -multiplicative 2-form (resp. closed multiplicative 2-form). If the matrix R is a solution of DQYBE, then the matrix R′ = n∑ i,j=1 ∆ije (n) ij ⊗ e (n) ji + n∑ i 6=j=1 αijdije (n) ii ⊗ e (n) jj . is also a solution of DQYBE. Proof. This is directly seen on system (S) and by remarking that the transformation defined above respects the D-classes, and then the ordered partition of the set of indices N∗n. (F1)–(F6) and (E2)–(E5) are factorized by d-coefficients, then either they are trivially verified (if iDj, when, for example, (F1)–(F6) are considered with indices ij) or the d-coefficients can be simplified (if i��D j, for the same example). (F7), (F8) and (E6) depend on d-coefficients only through the product dijdji, which is clearly invariant under the previous transformation, since (αij)(i,j)∈I(2,�D) (resp. (αij)(i,j)∈I2) is a ��D - multiplicative 2-form (resp. multiplicative 2-form). The same kind of argument applies to (E1), which is also invariant, since (αij)(i,j)∈I(2,�D) (resp. (αij)(i,j)∈I2) is in addition assumed to be ��D -closed (resp. closed). � Corollary 4.2. Let I be a subset of the set of indices N∗n of cardinality m, and (αij)(i,j)∈I(2,�D) (resp. (αij)(i,j)∈I2) be a ��D -closed ��D -multiplicative 2-form (resp. closed multiplicative 2-form). Following Proposition 2.3, the previous proposition implies that the contracted matrix (R′)I of the matrix R′ to the subset I is a solution of Glm(C)-DQYBE. Proposition 4.3. The family of non-zero functions (gij) (i,j)∈N∗(2,�D,<) n , introduced in Lemmas 4.3 and 4.4, is a ��D -closed ��D -multiplicative 2-form. Proof. This is an obvious corollary of the remark of Lemmas 4.3 and 4.4, since the family of non-zero functions (d0 ij) (i,j)∈N∗(2,�D,<) n introduced by Lemmas 4.3 and 4.4 trivially satisfies (E1). Then, using gijgji = 1, (E1) with indices ijk is simply the cyclic relation gij(k) gij gjk(i) gjk gki(j) gki = 1. � Corollary 4.3. The ��D -closed ��D -multiplicative 2-form (gij) (i,j)∈N∗(2,�D,<) n can be factorized out by ��D -multiplicative covariance, as described by Proposition 4.2. Proposition 4.4. Any (��D -)exact (��D -)multiplicative 2-form is (��D -)closed. Remark 4.2. Under the assumption that, for any (i, j) ∈ N∗(2,<) n , the non-zero function αij is a holomorphic function of the variable λ in a simply connected domain of Cn, there exists a multiplicative analog of the Poincaré lemma for differential forms, the so-called multiplicative Poincaré lemma. It enunciates that the reciprocal of Proposition 4.4 is also true, that is: a mul- tiplicative 2-form (αij)(i,j)∈N∗(2,<) n is exact, if and only if it is closed [13]. This directly implies that the multiplicative covariance of Proposition 4.2 coincides under this assumption with the twist covariance of Proposition 2.1. In particular, if I0 = N∗n, i.e. if the set of indices N∗n only contain free indices, a ��D -closed ��D -multiplicative 2-form is a closed multiplicative 2-form, and then is exact, which is the case e.g. for (weak) Hecke-type solutions of DQYBE5. In this case, by analogy with differential forms, the closed multiplicative 2-form (αij)(i,j)∈N∗(2,<) n will be refereed as a gauge 2-form, since it can 5For more details, see Subsection 5.5. 20 J. Avan, B. Billaud and G. Rollet universally be factorized out thanks to Proposition 2.1, in the sense that it is representation- independent. Considering the general problem, we do not succeed to solve whether or not any ��D -closed ��D - multiplicative 2-form is ��D-exact. However, as we will see in the proof of the following proposition, it does not really matter in practice, since the notion of ��D -multiplicative covariance is actually the minimal main tool allowing to achieve the characterization of the decoupled R-matrices. The issue which therefore remains is to get a general classification of ��D -closed ��D -multiplica- tive 2-forms, when the D-classes have a non-trivial structure. Note that if the set of indices N∗n is split into two D-classes any ��D -multiplicative 2-form is ��D -closed, since no cyclic relation exists. Proposition 4.5. Any R-matrix, solution of DQYBE, characterized by a block-upper triangu- lar reduced ∆-incidence matrix MR with two or more triangular blocks, is ��D -multiplicatively reducible to a decoupled R-matrix, and vice versa. Proof. This results from successive implementations of Proposition 2.2. Assuming s ≥ 2, let (q, q′) ∈ N∗(2,<) s and consider the solutions of system (S) restricted to the subsets Kq and Kq′ . According to Proposition 2.3, matrix elements of the R-matrix with both indices either in the subset Kq or in the subset Kq′ realize a contraction-type solution of a lower-dimensional, more precisely of a Nq-dimensional or Nq′-dimensional DQYBE. Due to the block-upper triangularity, the only remaining non-zero matrix elements are the d-coefficients dij , with (i, j) ∈ Kq ×Kq′ . Lemma 4.3 now solves this issue, for any pair of labels (q, q′) ∈ N∗(2,<) s . Indeed, Lem- mas 4.3 and 4.4, and Proposition 4.3 prove the existence of a ��D -closed ��D -multiplicative 2-form (gij) (i,j)∈N∗(2,�D,<) n , such that dij = √ Σqq′gij and dji = √ Σqq′gji, ∀ (i, j) ∈ Kq ×Kq′ , where Σqq′ is a non-zero constant, or dij = gijd 0 ij and dji = gjid 0 ji, ∀ (i, j) ∈ K(2,�D ,<) q , where the functions d0 ij and d0 ji are 1-periodic in the variable λk, for any k ∈ N∗n \Kq. Moreover, from Lemmas 4.1 and 4.2, the ∆-coefficients ∆ii′ , with i, i′ ∈ Kq, are also 1-periodic in the variable λk, for any k ∈ N∗n \Kq, with q ∈ N∗s. From Proposition 4.3 and its corollary, the ��D -multiplicative 2-form (gij) (i,j)∈N∗(2,�D,<) n is ��D - closed, and then can precisely be factorized out by the ��D -multiplicative covariance. This brings, on the one hand, the d-coefficients dij , with (i, j) ∈ Kq × Kq′ and (q, q′) ∈ N∗(2,<) s , to be equal to an overall block-pair dependent constant √ Σqq′ , and, on the other hand, the d-coefficient dij to be equal to d0 ij , for any (i, j) ∈ K(2,�D ) q . To summarize, any solution of DQYBE defined by its block-upper triangular reduced ∆- incidence matrix MR is necessarily ��D -multiplicatively covariant to a multiply decoupled R- matrix obtained from successive applications of Proposition 2.2. But this proposition shows that such decoupled R-matrices are also solutions of DQYBE. The reciprocal is obvious. � Corollary 4.4. It is therefore relevant to focus our discussion of solutions of system (S) to the cases ii and iii of Corollary 3.6, where the indices i, j ∈ Kq, with q ∈ N∗s and Nq ≥ 2. 4.3 Sum and determinant In this section, as stated in the following fundamental result, the functions Sij and Σij are shown to be actually constant independent of indices i and j, as soon as the pair of indices (i, j) is Classification of Non-Affine Non-Hecke Dynamical R-Matrices 21 a ��D -pair, depending then only on the subset Kq, i.e. only on the label q. Moreover, they actually parametrize the set of solutions of system (S) restricted to the subset Kq, to be specified later. Theorem 4.1 (inside a set Kq of ∆-classes). There exist a constant Sq and a non-zero con- stant Σq such that Sij = Sq and Σij = Σq, ∀ (i, j) ∈ K(2,�D ,<) q ; refereed as the “sum” and the “determinant” in the subset Kq. Moreover, denoting Dq = √ S2 q + 4Σq the “discriminant” in the subset Kq, there exists a fa- mily of ∑ p∈Pq (n (p) 0 + lp) signs (εI)I∈{I(i), i∈Kq} such that ∆I(i) = Sq + εI(i)Dq 2 , ∀ i ∈ Kq. (4.2) Proof. If n = 2, i.e. if N∗n = Kq = {i, j}, the theorem is a direct corollary of the proof of Lemma 4.4 and (4.1). Assuming that n ≥ 3, the proof of Lemma 4.4 and (4.1) shows that the functions Sij and Σij are 1-periodic in the variables λi, λj and λk, with k ∈ N∗n \ (I(i) ∪ I(j)). If Nq = 2, i.e. if Kq = {i, j}, the proof of the theorem ends. The proof goes now in three steps. 1. Periodicity. Assuming that Nq ≥ 3, it becomes possible to introduce a third index k ∈ Kq \ {i, j}. Two symmetrical possibilities k��D i or k��D j are to be considered. Indeed, ad absurdum, kDi and kDj leads to the contradiction iDj. • If k��D j, any d-coefficient involving one of the indices i, k and the index j is non-zero. Moreover, since i, k ∈ Kq, ∆ik 6= 0 or ∆ki 6= 0. Without loss of generality, it is possible to assume that ∆ik 6= 0, the case ∆ki 6= 0 being treated similarly by exchanging the indices i and k. Hence (E4) and (E5) both with indices ijk give ∆ij(k)∆jk + ∆ik∆ji(k) = ∆ik∆jk, (E′4) and ∆ij(k)∆jk + ∆ik(j)∆kj = ∆ik(j)∆ij(k). (E′5) From Proposition 4.1 (when kDi) and from Lemma 4.2 (when k��D i), the function ∆ik is 1-periodic in the variable λj . Then, by the substraction (E′4)–(E′5), we get that ∆ik{Sij(k)− Skj} = 0. From which we deduce that Sij(k) = Skj . However, we have seen that the function Skj is 1-periodic in the variable λk, since k��D j. The function Sij is thus also 1-periodic in the variable λk, for any k ∈ Kq \ ({i} ∪ I(j)). Moreover, we obtain that Skj = Sij and Σkj = ∆I(j){∆I(j) − Skj} = ∆I(j){∆I(j) − Sij} = Σij , ∀ (k, j) ∈ (Kq \ {i})(2,�D ). 22 J. Avan, B. Billaud and G. Rollet • If k��D i, the previous reasoning is symmetrically done, exchanging the indices i and j. This yields that the function Sij is 1-periodic in the variable λk, for any k ∈ Kq \ (I(i) ∪ {j}), and Sik = Sij and Σik = Σij , ∀ (i, k) ∈ (Kq \ {j})(2,�D ). If I(j) = {j}, or if I(i) = {i}, the functions Sij and Σij are thus in particular respectively 1-periodic in the variable λk, for any k ∈ I(i) \ {i} or k ∈ I(j) \ {j}, from application of the first or the second previous point, and then are constant from above. Assuming that |I(i)| ≥ 2 and |I(j)| ≥ 2, both previous points apply, implying that the functions Sij and Σij are actually 1-periodic in variable λk, with k ∈ (I(i) ∪ I(j)) \ {i, j}, and then are constant in this case too. This ends the proof of the periodicity property expressed in the theorem. 2. Existence of Sq and Σq. Always under the assumption that Nq ≥ 3, Step 1 has been seen to justify the existence of a ��D -pair of indices (i′, j′) ∈ K(2,�D ) q , distinct from the ��D -pair (i, j). Since these two ��D -pairs are distinct, it is always possible to impose that i′ ∈ Kq \ {i, j}. This suggests to rather adopt the notation (k, j′), where k ∈ Kq \ {i, j}. • If j′ ∈ {i, j}, then k��D j′��D i′, where we define the index i′ ∈ {i, j} so that {i′, j′} = {i, j}. From Step 1, if j′ = j, the case j′ = i being treated similarly by exchanging the indices i and j, we directly deduce that Skj′ = Sj′i′ = Sij and Σkj′ = Σij . If kDi′, the pairs of indices (k, j′) and (j′, i′) are the only ��D -pairs in {i′, j′, k}. On the contrary, if k��D i′, i.e. if (i′, j′, k) ∈ K(3,�D ) q , the pair of indices (k, i′) has to be also considered. The result we have just obtained applies to the indices k��D i′��D j′, leading to the second needed set of equations Ski′ = Sij and Σki′ = Σij . In particular, if Nq = 3, i.e. if Kq = {i, j, k}, the existence of the constant Sq is proved. • Assuming that Nq ≥ 4, since the first point of this reasoning already dealt with the case j′ ∈ {i, j}, we can consider here without loss of generality the case j′ /∈ {i, j}. However, since once more either k��D i or k��D j, and either j′��D i or j′��D j, there exists j1, j2 ∈ {i, j}, such that k��D j1 and j′��D j2. Defining the index i1 ∈ {i, j} so that {i1, j1} = {i, j}, then j2 ∈ {i1, j1}. If j2 = j1, then k��D j′��D j1��D i1, and the first point of this reasoning applies successively to the subsets {k, j′, j1} and {j′, j1, i1}, implying that Skj′ = Sj′j1 = Sj1i1 = Sij ⇒ Σkj′ = Σij . If j2 = i1, then k��D j′��D j2��D j1, and the first point of this reasoning applies successively to the subsets {k, j′, j2} and {j′, j2, j1}, implying that Skj′ = Sj′j2 = Sj2j1 = Sij ⇒ Σkj′ = Σij . This implies that there exists two constants Sq = Sij and Σq = Σij , such that Si′j′ = Sq and Σi′j′ = Σq, ∀ (i′, j′) ∈ K(2,�D ) q , ending the proof of the first part of the theorem. 3. Existence of (εI)I∈{I(i), i∈Kq}. The previous two steps now imply, from (4.1), that the family of constants (∆I)I∈{I(i), i∈Kq}, Sq and Σq satisfy the following quadratic equation ∆2 I(i) − Sq∆I(i) − Σq = 0. (4.2′) Classification of Non-Affine Non-Hecke Dynamical R-Matrices 23 From which we deduce the existence of a sign εI(i) ∈ {±} for each of the ∑ p∈Pq (n (p) 0 + lp) D-classes I(i) ⊆ Kq. � Remark 4.3. We have to insist on the fact that (4.2′) does not impose that the constant ∆I(i) and the sign εI(i) to be independent from the D-class I(i). Considering a ��D -pair of indices (i, j) ∈ K(2,�D ) q , ∆I(i) and ∆I(j) are solutions of (4.2′), which is equivalent to εI(i) = ±εI(j), and does not indeed constrain the family of signs (εI)I∈{I(i), i∈Kq}. This remark will be crucial later to distinguish between Hecke, weak Hecke and non-Hecke type solutions6. 4.4 Inside a ∆-class This section will present in details the explicit resolution of system (S) restricted to any ∆- class Jp, with p ∈ Pq, of the subset Kq. Case i of Corollary 3.6 being already solved in Theo- rem 4.1 thanks to Proposition 4.1, we have to focus on case ii, in which any pair of indices (i, j) under study is a ��D -pair. In general, the solution will be parametrized by the values of the sum Sq and the constant Tq = Dq−Sq Dq+Sq ∈ C∗ = C \ {0}, in addition with the family of signs (εI)I∈{I(i), i∈Kq}. More precisely we will see that three cases are to be distinguished. Remark 4.4. The quantity Tq is well defined and is non-zero, for any constants Sq and Σq 6= 0 (since, by construction, Dq = √ S2 q + 4Σq 6= ±Sq). Moreover, Sq = 0 if and only if Tq = 1. When Sq 6= 0, there exists tq = Dq Sq ∈ C \ {±1} such that Tq = −1−tq 1+tq = −1−|tq |2−2i=(tq) |1+tq |2 ∈ C∗, where the limit |tq| → ∞, or equivalently the limit Sq → 0, exists. We then deduce that Tq ∈ R∗−, if and only if tq ∈ ]−1, 1[. For later purpose, when Sq 6= 0 (when Tq 6= 1), we also introduce the non-zero constants Aq and Bq, viewed as functions of Tq, defined as Aq = { log(Tq) if Tq /∈ R∗−, log(−Tq)− iπ if Tq ∈ R∗−, and Bq = Sq +Dq 2 = Sq 1− eAq , where the principal value of the function log: C\R− −→ C is used when needed. When, Sq = 0, the constant Bq can also be defined, and is equal to √ Σq. Before beginning the resolution, we need to introduce the following technical lemma. Lemma 4.5 (multiplicative shift). Let A ∈ C, i ∈ Kq and a family of non-zero functions (βj)j∈I(i) of the variable λk, for any k ∈ I(i), and 1-periodic in any other variable, such that, for any j ∈ I(i) and j′ ∈ I(i) \ {j} βj(j) = eAεI(i)βj and βj(j ′) = eAεI(i)βj′ . (4.3) Then, there exists a non-zero constant fI(i) such that, for any j ∈ I(i) βj = βI(i) = eAεI(i)ΛI(i)fI(i), (4.3′) where we define the variable ΛI(i) = ∑ k∈I(i) λk. Proof. The case A = 0 being trivial, we will focus on the case A ∈ C∗. If |I(i)| = 1, i.e. I(i) = {i}, (4.3) reduces to βi(i) = eAεI(i)βi. Hence there exists a non-zero function fi, such that βi = eAεI(i)λifi, and the proof of the lemma ends. 6For more details, see Subsection 5.5, and particularly Propositions 5.6 and 5.7. 24 J. Avan, B. Billaud and G. Rollet Assuming that |I(i)| ≥ 2, from (4.3), for any j ∈ I(i) \ {i}, we deduce that βi(j) = eAεI(i)βj = βj(j), i.e. βi = βj . Hence there exists a non-zero function βI(i) of the variable λk, for any k ∈ I(i), and 1-periodic in any other variable, such that βI(i) = βj , ∀ j ∈ I(i). From (4.3), the function βI(i) satisfies βI(i)(j) = eAεI(i)βI(i), ∀ j ∈ I(i). (4.3′′) We now define the function fI(i) of the variable λk, for any k ∈ I(i), and 1-periodic in any other variable as fI(i) = e−AεI(i)βI(i)ΛI(i)βI(i). From (4.3′′), we directly deduce that fI(i) is now periodic in the variable λk, for any k ∈ I(i), and then is constant. � This result possesses an obvious linear limit. Lemma 4.6 (additive shift). Let i ∈ Kq and a family of functions (βj)j∈I(i) of the variab- le λk, for any k ∈ I(i), and 1-periodic in any other variable, such that, for any j ∈ I(i) and j′ ∈ I(i) \ {j} βj(j) = βj + εI(i) and βj(j ′) = βj′ + εI(i). (4.4) Then, there exists a constant fI(i) such that, for any j ∈ I(i) βj = βI(i) = εI(i)ΛI(i) + fI(i). (4.4′) Proof. Let a ∈ C and introduce, for any j ∈ I(i), the function of the variable λk, for any k ∈ I(i) βaj = eaβj . By construction, the family of functions (βaj )j∈I(i) satisfies the assumptions of Lemma 4.5. Hence, there exists a non-zero constant faI(i) such that βaj = βaI(i) = eaεI(i)ΛI(i)faI(i), ∀ j ∈ I(i). However, for any j ∈ I(i), the function a 7−→ βaj is holomorphic on C, and then as well as the functions a 7−→ βaI(i) and a 7−→ faI(i). These three functions admit a Taylor expansion in the neighboorhood of 0. In particular, there exists a constant fI(i) = d daf a I(i) ∣∣ a=0 such that βaj = f0 I(i) + aβj + o(a) = f0 I(i) + a(εI(i)ΛI(i) + fI(i)) + o(a) ⇒ βj = εI(i)ΛI(i) + fI(i) = βI(i). � We now enunciate the fundamental result of the resolution in any ∆-class Jp as well as in Theorems 4.4 and 5.1, as justified by the special treatment of Propositions 4.2, 4.3 and 4.4. Theorem 4.2 (trigonometric behavior). Assuming that Sq 6= 0, i.e. Tq ∈ C∗ \ {1}, there exist a family of n (p) 0 + lp non-zero constants (fI)I∈{I(i), i∈Jp} (with the convention that fI(min Jp) = 1), and a ��D -multiplicative 2-form (gij) (i,j)∈J(2,�D,<) p , such that the solution of system (S) restricted to the ∆-class Jp is given by the following expressions ∆I(i) = Sq 1− eAqεI(i) , ∀ i ∈ Jp; Classification of Non-Affine Non-Hecke Dynamical R-Matrices 25 ∆ij((λk)k∈I(i)∪I(j)) = Sq 1− eAq(εI(i)ΛI(i)−εI(j)ΛI(j)) fI(i) fI(j) = ∆I(i)I(j)(ΛI(i),ΛI(j)) and dij = gij{Bq −∆I(i)I(j)}, ∀ (i, j) ∈ J(2,�D ) p . Proof. 1. Diagonal ∆-coeff icients. From (4.2), if we adopt for a time the notations A±q = log(±Tq), where the exponent “±” means respectively that Tq /∈ R∗− or Tq ∈ R∗−, we deduce that 1− Sq ∆I(i) = εI(i)Dq − Sq εI(i)Dq + Sq = T εI(i) q = ±eA ± q εI(i) = eAqεI(i) ⇔ ∆I(i) = Sq 1∓ eA ± q εI(i) = Sq 1− eAqεI(i) . For the rest of the article, we will omit to make the explicit split between the cases Tq /∈ R∗− and Tq ∈ R∗−, unless otherwise stated. If np = 1 or if the ∆-class Jp is reduced to a single D-class (cf. Proposition 4.1), the proof of the theorem ends here (cf. case i of Corollary 3.6). 2. Off-diagonal ∆-coeff icients. Assuming that np ≥ 2 and that Jp is not reduced to a D-class, there exists a ��D -pair of indices (i, j) ∈ J(2,�D ) p . Let (i′, j′) ∈ J(2,�D ) p be a ��D -pair of indices such that (i′, j′) ∈ I(i) × I(j). From Lemmas 4.1 and 4.2, only the dependence in the variable λk, for any k ∈ I(i) ∪ I(j), of the function ∆i′j′ remains to be determined. To this end, (F ′1)⇔ (F ′2) with indices i′j′ is re-written, since ∆I(i)∆I(j)∆I(i)I(j)∆I(j)I(i) 6= 0, as ∆I(i) ∆i′j′(i′) − 1 = { 1− Sq ∆I(i) } ∆I(i) ∆i′j′ = eAqεI(i) ∆I(i) ∆i′j′ ⇔ 1 ∆i′j′(i) = eAqεI(i) ∆i′j′ + 1 ∆I(i) . Denoting βi′j′ = Sq ∆i′j′ − 1 = ∆j′i′ ∆i′j′ 6= 0, we deduce that βi′j′(i ′) = eAqεI(i) Sq ∆i′j′ + Sq ∆I(i) − 1 = eAqεI(i)βi′j′ , and by symmetry βi′j′(j ′) = e−AqεI(j)βi′j′ . • If |I(j)| = 1, i.e. if I(j) = {j}, the only ��D -pairs to consider are (i′, j), for any i′ ∈ I(i). From above, the function βi′j satisfies βi′j(i ′) = eAqεI(i)βi′j and βi′j(j) = e−AqεI(j)βi′j . • If |I(j)| ≥ 2, i.e. if I(j) is a D-class, let k ∈ I(j) and k′ ∈ I(j) \ {k}. Since k′��D i′��D k, for any i′ ∈ I(i), we have βi′k(k) = e−AqεI(j)βi′k and βi′k′(k ′) = e−AqεI(j)βi′k′ . Moreover, (E′4) can be used with indices i′kk′, and yields βi′k(k ′) = e−AqεI(j)βi′k′ . (4.5) In both cases, Lemma 4.5 is now applied to the family of functions (βi′k)k∈I(j) of the variable λk, for any k ∈ I(i) ∪ I(j). Hence, there exists a non-zero function βi′I(j) of the variable λk, for any k ∈ I(i), such that, for any i′ ∈ I(i) βi′j′ = e−AqεI(j)ΛI(j)βi′I(j) and βi′I(j)(i ′) = eAqεI(i)βi′I(j), ∀ j′ ∈ I(j). (4.6) 26 J. Avan, B. Billaud and G. Rollet • If |I(i)| = 1, i.e. if I(i) = {i}, the only ��D -pairs to consider are (i, j′), for any j′ ∈ I(j). From above, the function βiI(j) satisfies βiI(j)(i) = eAqεI(i)βiI(j). • If |I(i)| ≥ 2, let k ∈ I(i) and k′ ∈ I(i) \ {k}. The previous reasoning ensures the existence of a non-zero function βkI(j) of the variable λk, for any k ∈ I(i), which satisfies (4.6). Moreover, the exchange of the indices i and j, as well as the indices i′ and j′ in (4.5) yields βkI(j)(k ′) = eAqεI(i)βk′I(j), where the symmetry relation βi′j′ = 1 βj′i′ is used, for any (i′, j′) ∈ J(2,�D ) p | (i′, j′) ∈ I(i)× I(j). In both cases, Lemma 4.5 applies once more to the family of functions (βkI(j))k∈I(i) of the variable λk, for any k ∈ I(i), ensuring the existence of the non-zero constant fI(i)I(j) such that βi′j′(ΛI(i),ΛI(j)) = eAq(εI(i)ΛI(i)−εI(j)ΛI(j))fI(i)I(j), ∀ (i′, j′) ∈ I(i)× I(j). Finally, this implies that ∆ij((λk)k∈I(i)∪I(j)) = Sq 1− eAq(εI(i)ΛI(i)−εI(j)ΛI(j))fI(i)I(j) = ∆I(i)I(j)(ΛI(i),ΛI(j)), ∀ (i, j) ∈ J(2,�D ) p ; where, as in the second point above, by symmetry, fI(i)I(j)fI(j)I(i) = 1. 3. Existence of the functions (fI)I∈{I(i), i∈Jp}. If Jp = I(i) ∪ I(j), assuming that i < j, it is sufficient to set fI(i) = 1 and fI(j) = 1 fI(i)I(j) . Assuming that Jp 6= I(i) ∪ I(j), let k ∈ Jp \ (I(i) ∪ I(j)), i.e. (i, j, k) ∈ J(3,�D ) p . (E′4) can be used with indices ijk, and yields, by linear independence of the functions ΛI(i) 7−→ eAqεI(i)ΛI(i) , ΛI(j) 7−→ eAqεI(j)ΛI(j) and ΛI(k) 7−→ eAqεI(k)ΛI(k) 1 ∆I(j)I(i)∆I(i)I(k) + 1 ∆I(i)I(j)∆I(j)I(k) = 1 ∆I(i)I(j)∆I(i)I(k) ⇔ fI(i)I(k)fI(k)I(j) = fI(i)I(j). This second set of equations reduces the number of independent constants to the choice of a family of 2(lp + n (p) 0 ) non-zero constants (fIJ, fJI)I,J∈{I(i), i∈Jp}, for any fixed D-class I, for example I(imin), where imin = min Jp. The first set of equations fI(i)I(j)fI(j)I(i) = 1 reduces this number by half. The family of lp + n (p) 0 non-zero constants (fI)I∈{I(i), i∈Jp} (fixed by fI(imin) = 1 and fI(i) = 1 fI(i)I(imin) , for any i ∈ Jp \ {imin}) satisfies the expected properties. 4. d-coeff icients. Setting Bq = Sq 1−eAq = Bij , for any ��D -pair (i, j) ∈ J(2,�D ,<) p , Lemma 4.4 and Proposition 4.3 insure the existence of a ��D -closed ��D -multiplicative 2-form (gij) (i,j)∈J(2,�D,<) p such that dij = gijd 0 I(i)I(j), where d0 I(i)I(j) = Bq −∆I(i)I(j). Reciprocally, it is straightforward to check that the family of constants (∆I)I∈{I(i), i∈Jp} and the family of functions (∆I(i)I(j),∆I(j)I(i), dij , dji) (i,j)∈J(2,�D) p are indeed solutions of system (S) restricted to the ∆-class Jp. Note in particular that, as mentioned in the remark of Lemmas 4.3 and 4.4, the family of functions (dij) (i,j)∈J(2,�D,<) p obeys (E1). The set of solutions of system (S) restricted to any ∆-class Jp is exactly parametrized by the giving of the constants Sq, Σq and (εI, fI)I∈{I(i), i∈Jp} and the ��D -multiplicative 2-form (gij) (i,j)∈J(2,�D,<) p . � Classification of Non-Affine Non-Hecke Dynamical R-Matrices 27 Remark 4.5. Because we use the principal value of the logarithm function, if Tq ∈ C∗\{1}∪R∗− then =(Aq) = arg(Tq) ∈ ]−π, π[, and if Tq ∈ R∗−, then =(Aq) = −π. Therefore, if Tq ∈ C∗ \ {1} ∪ R∗−, the function ΛI(i) 7−→ eAqεI(i)ΛI(i) can be periodic of any period strictly greater than 2, but cannot be 2-periodic (since Aq 6= 0, the period is greater than 2). This happens if and only if Tq ∈ R∗−, in which case can arise 2-periodic trigonometric functions such as ΛI(i) 7−→ eiπεI(i)ΛI(i) . Moreover, Theorem 4.2 justifies the choosing of the quantity Aq, through the choosing of a particular complex logarithm. The expressions obtained for the solutions are indeed indepen- dent of this choice. Strictly speaking, =(Aq) may be a priori defined up to 2π. However, for any k ∈ Z, the function ΛI(i) 7−→ e2iπkεI(i)ΛI(i) is 1-periodic in any variable. Remembering that “constant quantity” means in fact “1-periodic function in any variable”, it can be re-absorbed in each constant of the family (fI)I∈{I(i), i∈Jp} by multiplying each one by the function ΛI(i) 7−→ e−2iπk(εI(i)ΛI(i)−εI(min Jp)ΛI(min Jp)). In particular, this preserves the convention fI(min Jp) = 1.7 Then, the universal covering of C∗ by the Riemann surface S = {(z, θ) ∈ C∗ × R | θ − arg(z) ∈ 2πZ}, associated with the logarithm function logS: (z, θ) ∈ S 7−→ log |z|+ iθ, allows to naturally con- tinuously extend expressions of Theorem 4.2, viewed as functions of Tq, to the surface S. This can be done as above by multiplying each constant fI(i), viewed as a function of Tq ∈ S, by the function Tq ∈ S 7−→ ei(θ−arg(Tq))(εI(i)ΛI(i)−εI(min Jp)ΛI(min Jp)), which is 1-periodic in any variable, for any Tq ∈ S. Theorem 4.3 (rational behavior). Assuming that Sq = 0, i.e. Tq = 1, there exist a family of n (p) 0 +lp constants (fI)I∈{I(i), i∈Jp} (with the convention that fI(min Jp) = 0), and a ��D -multiplicative 2-form (gij) (i,j)∈J(2,�D,<) p , such that the solution of system (S) restricted to the ∆-class Jp is given by the following expressions ∆I(i) = εI(i) √ Σq, ∀ i ∈ Jp; ∆ij((λk)k∈I(i)∪I(j)) = −∆ji({λk}k∈I(i)∪I(j)) = ∆I(i)I(j)(ΛI(i),ΛI(j)) = −∆I(j)I(i)(ΛI(j),ΛI(i)) = √ Σq εI(i)ΛI(i) − εI(j)ΛI(j) + fI(i) − fI(j) and dij = gij{Bq −∆I(i)I(j)}, ∀ (i, j) ∈ J(2,�D ) p . Proof. The proof of Theorem 4.2 to compute the functions ∆I(i)I(j) and dij , for any i, j ∈ Jp can be directly adapted here, since the family of constants (∆I)I∈{I(i), i∈Jp} are obviously obtained from (4.2). If np = 1 or if the ∆-class Jp is reduced to a D-class (cf. Proposition 4.1), the proof of the theorem ends here (cf. case i of Corollary 3.6). Assuming that np ≥ 2 and that Jp is not reduced to a D-class, there exists a ��D -pair of indices (i, j) ∈ J(2,�D ) p . Let (i′, j′) ∈ J(2,�D ) p be a ��D -pair of indices such that (i′, j′) ∈ I(i)× I(j). From Lemmas 4.1 and 4.2, only the dependence in the variable λk, for any k ∈ I(i)∪ I(j), of the function ∆i′j′ remains to be determined. To this end, (F ′1)⇔ (F ′2) with indices i′j′ is re-written, since ∆I(i)∆I(j)∆I(i)I(j) 6= 0 and denoting hij = √ Σq ∆ij = −βji, as hij(i) = hij + εI(i) and hij(j) = hij − εI(j), thanks to the equality Sq = ∆ij + ∆ji = 0. Moreover, for any k ∈ I(j) and any k′ ∈ I(j) \ {j}, (E′4) can be used with indices i′kk′, and yields βi′k(k ′) = βi′k′ − εI(j), 7Such manipulation will be also used in the proof of Propositions 5.2 and 5.3. 28 J. Avan, B. Billaud and G. Rollet the exact correspondent of (4.5), to which we apply Lemma 4.6. We deduce the existence of a constant fI(i)I(j), such that βi′j′(ΛI(i),ΛI(j)) = εI(i)ΛI(i) − εI(j)ΛI(j) + fI(i)I(j), ∀ (i′, j′) ∈ I(i)× I(j). Finally, this implies that ∆ij((λk)k∈I(i)∪I(j)) = √ Σq εI(i)ΛI(i)− εI(j)ΛI(j)+ fI(i)I(j) = ∆ij(ΛI(i),ΛI(j)), ∀ (i, j) ∈ J(2,�D ) p , where, as above, by symmetry, fI(i)I(j) = −fI(j)I(i). If Jp = I(i) ∪ I(j), assuming that i < j, it is sufficient to set fI(i) = 0 and fI(j) = −fI(i)I(j). Assuming that Jp 6= I(i) ∪ I(j), let k ∈ Jp \ (I(i) ∪ I(j)), i.e. (i, j, k) ∈ J(2,�D ) p . Hence (E′4) can be used with indices ijk, and yields, by linear independence of the functions ΛI(i) 7−→ ΛI(i), ΛI(j) 7−→ ΛI(j) and ΛI(k) 7−→ ΛI(k) 1 ∆I(j)I(i)∆I(i)I(k) + 1 ∆I(i)I(j)∆I(j)I(k) = 1 ∆I(i)I(j)∆I(i)I(k) ⇔ fI(i)I(k) + fI(k)I(j) = fI(i)I(j). This second set of equations reduces the number of independent constants to the choice of a family of 2(lp + n (p) 0 − 1) constants (fIJ, fJI)I,J∈{I(i), i∈Jp}, for any fixed D-class I, for example I(imin), where imin = min Jp. The family (fI)I∈{I(i), i∈Jp} defined as fI(imin) = 0 and fI(i) = −fI(i)I(imin), for any i ∈ Jp \ {imin}, yields the expected result, the rest of the proof being identical to the proof of Theorem 4.2. � 4.5 Inside a subset Kq We now end the resolution of system (S) by solving case iii of Corollary 3.6. Let q ∈ N∗s, and consider a subset Kq = ⋃ p∈Pq Jp such that rq ≥ 2, the case rq = 1 being already treated in Theorems 4.2 and 4.3. We now have to determine the cross-terms between two distinct ∆-classes. This is given by the following theorem. Theorem 4.4 (trigonometric behavior). Let q ∈ N∗s such that rq ≥ 2. Then, there exist two non-zero constant Sq and Σq, a family of signs (εI)I∈{I(i), i∈Kq}, a family of non-zero constants (fI)I∈{I(i), i∈Kq} (with the convention that fI(min Jp) = 1, for any p ∈ Pq), and a ��D -multiplicative 2-form (gij) (i,j)∈K(2,�D,<) q , such that the R-matrix, solution of system (S) restricted to the sub- set Kq, is given by R(q) = ∑ p∈Pq R(q) p + Sq ∑ (p,p′)∈P(2,<) q ∑ (i,j)∈Jp×Jp′ e·ij ⊗ e·ji + √ Σq ∑ (p,p′)∈P(2,<) q ∑ (i,j)∈Jp×Jp′ {gije·ii ⊗ e·jj + gjie · jj ⊗ e·ii}, where the family of matrices (R (p) q , (R (q,p) I(i) )i∈Jp) is defined, for any p ∈ Pq, as R(q) p = ∑ i∈Jp R (q,p) I(i) + ∑ (i,j)∈J(2,�D) p { ∆I(i)I(j)e · ij ⊗ e·ji + dije · ii ⊗ e·jj } = ∑ i∈Jp R (q,p) I(i) + ∑ (i,j)∈J(2,�D) p { ∆I(i)I(j)(e · ij ⊗ e·ji − gije·ii ⊗ e·jj) Classification of Non-Affine Non-Hecke Dynamical R-Matrices 29 + ∆I(j)I(i)(e · ji ⊗ e·ij − gjie·jj ⊗ e·ii) } +Bq ∑ (i,j)∈J(2,�D) p {gije·ii ⊗ e·jj + gjie · jj ⊗ e·ii}, and R (q,p) I(i) = ∆I(i) |I(i)| ∑ j,j′∈I(i) e·jj′ ⊗ e·j′j , ∀ i ∈ Jp. Remark 4.6. In this notation, the exponent “ · ” has to be chosen to appropriately specify the size of the matrices R (q) p or R(q). According to Proposition 2.3, there are three possibilities for the matrix R (q) p • if “ · = Jp”, then R (q) p is a solution of the Glnp(C)-DQYBE, which is the contraction to the ∆-class Jp of a R-matrix, solution of the Gln(C)-DQYBE or the GlNq(C)-DQYBE; • if “ · = Kq”, then R (q) p is the restriction to the ∆-class Jp of a solution of the GlNq(C)- DQYBE, which is the contraction to the subset Kq of a R-matrix, solution of the Gln(C)- DQYBE; • if “ · = (n)”, then R (q) p is the restriction to the ∆-class Jp of a R-matrix, solution of the Gln(C)-DQYBE. Similarly, there are two possibilities for the matrix R(q) • if “ · = Kq”, then R(q) is a solution of the GlNq(C)-DQYBE, which is the contraction to the subset Kq of a R-matrix, solution of the Gln(C)-DQYBE; • if “ · = (n)”, then R(q) is the restriction to the subset Kq of a R-matrix, solution of the Gln(C)-DQYBE. Proof. The existence of the constants Sq and Σq is insured by Theorem 4.1. Assuming that rq ≥ 2, for any ∆-class Jp, with p ∈ Pq, of the subset Kq, Theorems 4.2 and 4.3 give the expressions of the constants of the family (∆I)I∈{I(i), i∈Jp} and the functions of the family (∆I(i)I(j),∆I(j)I(i), dij , dji) (i,j)∈J(2,�D) p , depending on the signs of the family (εI)I∈{I(i), i∈Jp}, the constants of the family (fI)I∈{I(i), i∈Jp}, and the ��D -multiplicative 2-form (gij) (i,j)∈J(2,�D,<) p . It remains to determine the “crossed” functions (∆ij , dij , dji)(i,j)∈Jp×Jp′ |(p,p′)∈P (2,<) q of the solutions of system (S) restricted to the subset Kq. We also need to specify the dependence of the ��D - multiplicative 2-form (gij) (i,j)∈J(2,�D,<) p on the variable λk, for any k ∈ Kq \ Jp, which is done in Propositions 4.2, 4.3 and 4.4. From Theorem 4.1, we have ∆ij = Sq = ∆I(i)I(j), for any (i, j) ∈ Jp × Jp′ ∣∣ (p, p′) ∈ P(2,<) q . This implies also that the constant Sq is non-zero (otherwise ∆ij = Sq = 0, for any (i, j) ∈ Jp × Jp′ ∣∣ (p, p′) ∈ P(2,<) q ), which yields a contradiction with the construction of the subset Kq. The d-coefficients are deduced from the fact that d0 ij = d0 ji = √ Σq is a particular solution of dijdji = Σq, for any (i, j) ∈ Jp× Jp′ ∣∣ (p, p′) ∈ P(2,<) q , as in the proof of Theorem 4.2. This yields the expected result, the rest of the proof being almost identical. Reciprocally, it is straightforward to check that the family of functions (∆I(i)I(j),∆I(j)I(i), dij , dji) (i,j)∈K(2,�D) q are indeed solutions of system (S) restricted to the subset Kq. The set of solutions of system (S) restricted to the subset Kq is exactly parametrized by the giving of the constants Sq, Σq and (εI, fI)I∈{I(i), i∈Kq}, and the ��D -closed ��D -multiplicative 2-form (gij) (i,j)∈K(2,�D,<) q . This concludes the proof of the theorem. � 30 J. Avan, B. Billaud and G. Rollet Remark 4.7. Let us first insist on the fact that the proof above has been seen to justify that assuming that rq ≥ 2 implies that Sq 6= 0, and then forbids the rational behavior. In particular, if the rational behavior is assumed, i.e. if Sq = 0, then rq = 1, meaning that there exists a single label p ∈ N∗r such that Kq = Jp. A trigonometric R-matrix, solution of DQYBE restricted to the subset Kq, shows similitudes with a decoupled R-matrix presented in Proposition 2.2. As in this case, the d-coefficients dij and dji, for any (i, j) ∈ Jp × Jp′ , between two distinct subsets Jp and Jp′ , with (p, p′) ∈ P(2,�D ) q , are given by an overall Kq-dependent constant, up to ��D -multiplicative covariance (cf. Propositions 4.2, 4.3 and 4.4). The coupling (with ∆-coefficients) between the subsets Jp and Jp′ is minimal, in the sense that the coupling part ∑ (p,p′)∈P(2,<) q ∑ (i,j)∈Jp×Jp′ e·ij ⊗ e·ji of the R-matrix is non-dynamical. It is interesting to note a similarity of upper-triangularity structure of this non-dynamical part of the dynamical R-matrix with the Yangian R-matrix RYang = ∑ i<j eij⊗eji. structure also appears in the non-dynamical operators {R(q,p) I(i) }i∈Jp , describing the coupling inside any D-class. 5 General solution and structure of the set of solutions Given an ordered partition of the set of indices N∗n as described in Theorem 3.1, Proposition 4.5 and Theorems 4.1, 4.2, 4.3 and 4.4 allow to directly write the general form of any solution of DQYBE, compatible with this partition, where the only parameters of the R-matrix not explicitly constructed is the ��D -closed ��D -multiplicative 2-form (gij) (i,j)∈N∗(2,�D,<) n More precisely, we recall that Theorem 3.1 yields a first, set-theoretical, “parametrization” of the solutions of DQYBE, which is given in terms of an ordered partition of the indices set N∗n into the s ordered subsets {Kq, q ∈ N∗s}, being unions of the ∆-classes {Jp, p ∈ N∗r}, and the ordered partition of each ∆-class Jp into lp D-classes {I(i), i ∈ Jp}, either reduced to a single element (case of a D-class I(i) reduced to a free index i) or non-trivial (case of D-class I(i) generated by a non-free index i). Theorem 5.1 (general R-matrices). Let n ≥ 2 and an ordered partition of the set N∗n. Then, there exist a family of constants (Sq)q∈N∗s (with Sq 6= 0 if rq ≥ 2) two families of non-zero- constants (Σq)q∈N∗s and (Σqq′)(q,q′)∈N∗(2,<) s , a family of signs (εI)I∈{I(i), i∈N∗n}, a family of non- zero constants (fI)I∈{I(i), i∈N∗n} (with the convention that fI(min Jp) = 0, for any p ∈ Pq, with q ∈ N∗s | rq = 1, and fI(min Jp) = 1, for any p ∈ Pq, with q ∈ N∗s | rq ≥ 2), and a ��D -closed ��D -multiplicative 2-form (gij) (i,j)∈N∗(2,�D,<) n , such that the R-matrix, solution of DQYBE, is given by R = s∑ q=1 R(q) + ∑ (q,q′)∈N∗(2,<) s √ Σqq′ ∑ (i,j)∈Kq×Kq′ { gije (n) ii ⊗ e (n) jj + gjie (n) jj ⊗ e (n) ii } . Let us now characterize the structure of the moduli space of DQYBE. Putting aside the delicate issue of general ��D -closed ��D -multiplicative 2-forms, we see from Proposition 4.5 that the general solution of DQYBE is therefore built, up to the ��D -multiplicative covariance, in terms of solutions of DQYBE restricted to each subset Kq together with cross- terms (Σqq′)(q,q′)∈N∗(2,<) s between each pair of such subsets. This takes care of the interpretation of set-theoretical parameters {Kq, q ∈ N∗s} and c-number complex parameters (Σqq′)(q,q′)∈N∗(2,<) s as defining irreducible components of decoupled R-matrices according to Proposition 2.2. Classification of Non-Affine Non-Hecke Dynamical R-Matrices 31 5.1 Continuity properties of the solutions with respect to the constants (Sq,Σq)q∈N∗ s Consider the family of c-number complex parameters (Sq,Σq)q∈N∗s . As seen in Theorems 4.2, 4.3 and 4.4, the solution of system (S) restricted to any subset Kq is essentially characterized by the constants Sq and Σq 6= 0, through the quantity Tq 6= 0, except for the family of signs (εI)I∈{I(i), i∈Kq}, the set of constants (fI)I∈{I(i), i∈Kq} and the ��D -multiplicative 2-form (gij) (i,j)∈K(2,�D,<) q . We have particularly exhibited three cases to be distinguished i. Rational behavior: Sq = 0 and Tq = 1; ii. Trigonometric behavior (periodicity 2): Sq 6= 0 and Tq ∈ R∗−; iii. Trigonometric behavior (arbitrary periodicity strickly greater than 2): Sq 6= 0 and Tq ∈ C∗ \ {1} ∪ R∗−. This naturally rises the question whether these three types of solutions are distinct or if it is possible to connect them one to each other, typically in this situation by continuity arguments. Such connections exist and are described by the two following propositions. Let us immediately point out that the first two cases are clearly incompatible, and thus es- sentially different, in the sense that exploring the neighboorhood of the dimensionless variable Sq√ Σq = 0 imposes equivalently to explore the neighboorhood of Tq = 1, which cannot be asymp- totically reached by points in R∗−. In other words, for any fixed Σq 6= 0, the quantity Tq viewed as a function of Sq√ Σq defined on C is not continuous in 0, and solutions parameterized by Sq = 0 and Tq = 1 cannot be approached by solutions with Sq 6= 0 and Tq ∈ R∗−. The periodicity 2 of a trigonometric solution cannot then become infinite as required by the rational behavior. Remark 5.1. From Corollary 4.2, the ��D -multiplicative 2-form (gij) (i,j)∈K(2,�D,<) q can be fac- torized out independently in each subset Kq. In other words, the ��D -multiplicative 2-form (gij) (i,j)∈K(2,�D,<) q are non-relevant moduli to consider inside any fixed irreducible component. Hence, unless otherwise stated, we will assume in this section that gij = 1, for any (i, j) ∈ K(2,�D ,<) q . Let q ∈ N∗s such that rq = 1, a solution (parametrized by the constants S0 q = 0, Σq 6= 0 and (εI, f 0 I )I∈{I(i), i∈Kq}) of system (S) restricted to the subset Kq, and a parameter ξ ∈ C∗. Proposition 5.1 (from trigonometric to rational). There exists a solution (parametrized by the constants Sξq 6= 0, Σq and (εI, f ξ I )I∈{I(i), i∈Kq}) of system (S) restricted to the subset Kq, such that the solution with S0 q = 0 is the limit of the solution with Sξq 6= 0, when ξ → 0. In particular, if T ξq /∈ R∗−, for any ξ ∈ C∗, the piecewise solution (parametrized among others by the constants Sξq and (f ξI )I∈{I(i), i∈Kq}) is a continuous function of ξ on C. Proof. Since rq = 1, there exists an unique p ∈ N∗r such that Kq = Jp. The idea of the proof is to Taylor expand each considered quantity in the neighboorhood of ξ = 0, starting from Sξq . Such expansion exists as soon as the quantity under study is a sufficiently regular function of ξ ∈ C (or at least in a neighboorhood of ξ = 0). Let for example Sξq a non-zero holomorphic function of ξ on C∗ (or at least in a neighboorhood of ξ = 0), such that Sξq = Sqξ + o(ξ), with Sq = d dξS ξ q ∣∣ ξ=0 6= 0. Then, we deduce the following expansions Dξ q = 2 √ Σq{1 + o(ξ)} 6= 0 or T ξq = 1− Sq√ Σq ξ + o(ξ) /∈ R∗−, 32 J. Avan, B. Billaud and G. Rollet implying that Aξq = − Sq√ Σq ξ + o(ξ). Now, we reason similarly for the constants (f ξI )I∈{I(i), i∈Kq}. We assume that their Taylor expansion in the neighboorhood of ξ = 0 exists and stands, for any i ∈ Kq, f ξI(i) = 1− f0 I(i) Sq√ Σq ξ + o(ξ). (in particular, f ξI(i) = 1, for any ξ ∈ C∗, implies f0 I(i) = 0, and then we verify that f0 I(minKq) = 1, as required by Theorem 4.3). Then, we deduce that, for any i ∈ Kq ∆ξ I(i) = Sξq 1− eA ξ qεI(i) −−−→ ξ→0 εI(i) √ Σq = ∆0 I(i); and, for any (i, j) ∈ K(2,�D ) q ∆ξ I(i)I(j)(ΛI(i),ΛI(j)) = Sξq 1− eA ξ q(εI(i)ΛI(i)−εI(j)ΛI(j)) fξI(i) fξI(j) −−−→ ξ→0 √ Σq εI(i)ΛI(i) − εI(j)ΛI(j) + f0 I(i) − f 0 I(j) = ∆0 I(i)I(j)(ΛI(i),ΛI(j)). The limit for d-coefficients is trivially deduced from above, for any ��D -pair (i, j) ∈ K(2,�D ) q . � Corollary 5.1. Let s ≥ 2 and m ∈ N∗s. This proposition can be extended to any set m⋃ a=1 Kqa, where rqa = 1, for any a ∈ N∗m, providing that there exists a non-zero constant Σq, such that Σqa = Σqaqb = Σq, for any (a, b) ⊆ N∗(2,<) m . Proof. The case m = 1 already being treated, let assume that m ≥ 2. Since, from the remark of Theorem 4.4, there exists an unique pa ∈ N∗r such that Kqa = Jpa , the solution, parametrized (among others) by the constants Sξq 6= 0, Σq and (εI, f ξ I )I∈{I(i), i∈Kq}, of system (S) restricted to the subset Kq = n⋃ a=1 Jpa has the expected properties. The limit for coefficients in the family (∆ξ I(i)I(j), d ξ ij) (i,j)∈K(2, �D ) q are given by Proposition 5.1. The limit for the other ∆-coefficients is not problematic, since ∆ξ I(i)I(j) = Sq → 0 = ∆0 I(i)I(j), for any (i, j) ∈ Jpa × Jpb , with (a, b) ∈ N∗(2,<) m , and the limit for d-coefficients is trivially deduced from above. � Let q ∈ N∗s, a solution (parametrized by the constants S0 q 6= 0, Σ0 q , T 0 q = −1−t0q 1+t0q ∈ R∗−, i.e. t0q = D0 q S0 q ∈ ]−1, 1[, and (εI, f 0 I )I∈{I(i), i∈Kq}) of system (S) restricted to the subset Kq, and a parameter θ ∈ C \ R. Proposition 5.2 (from arbitrary periodicity to periodicity 2). There exists a solution (para- metrized by the constants Sθq 6= 0, Σθ q, T θq = −1−tθq 1+tθq /∈ R∗−, i.e. tθq = Dθq Sθq /∈ ]−1, 1[, and (εI, f θ I )I∈{I(i), i∈Kq}) of system (S) restricted to the subset Kq, such that the solution with T 0 q ∈ R∗− is the limit of the solution with T θq /∈ R∗−, when θ → 0. In particular, the piecewise solution (parametrized among others by the constants Sθq , Σθ q and (fθI )I∈{I(i), i∈Kq}) is a continuous function of θ on a neighboorhood of 0 in C \ R∗. Classification of Non-Affine Non-Hecke Dynamical R-Matrices 33 Proof. Let t0q ∈ ]−1, 1[. The proof consists in showing that the constants S0 q , Σ0 q and (f0 I )I∈{I(i), i∈Kq} parametrize a solution, characterized by t0q ∈ ]−1, 1[, which is the limit of a solution, to be precised, parametrized by the constants Sθ 6= 0, Σθ q and (fθI )I∈{I(i), i∈Kq}, and characterized by tθq /∈ ]−1, 1[, when θ → 0. Let for example Sθq and Σθ q two non-zero holomorphic functions of θ on a neighboorhood of θ = 0 in C \ R, such that Sθq = S0 q { 1− θ 2t0q + o(θ) } and Σθ q = Σ0 q { 1 + (D0 q)2+(S0 q )2 4Σ0 q θ t0q + o(θ) } . Then, using previous notation of the proof of Theorem 4.2, i.e.A0 q = log(−T 0 q )−iπ = log 1−θ0q 1+θ0q −iπ, we Taylor expand the quantity T θq = eA 0 q { 1− 2θ 1− (t0q) 2 + o(θ) } /∈ R∗−. From which, introducing the Heaviside function H, we deduce that Aθq = log(T θq ) = A0 q + 2iπH(=(θ)) + o(1), which does not converge in the general case when taking the limit θ → 0. However, it is possible to appropriately choose the functions (fθI )I∈{I(i), i∈Kq}, viewed as functions of θ, so that the constants (∆θ I )I∈{I(i), i∈Kq} and the functions (∆θ I(i)I(j))(i,j)∈K(2,�D) q converge to the constants (∆0 I )I∈{I(i), i∈Kq} and the functions (f0 I )I∈{I(i), i∈Kq} in this limit. More precisely, we have that, for any i ∈ Kq ∆θ I(i) = Sθq 1− eA θ qεI(i) = Sθq 1− eA 0 qεI(i)eo(1) −−−→ θ→0 S0 q 1− eA 0 qεI(i) = ∆0 I(i); and, for any (i, j) ∈ J(2,�D ) p , with p ∈ Pq, ∆θ I(i)I(j)(ΛI(i),ΛI(j)) = Sθq 1− eA θ q(εI(i)ΛI(i)−εI(j)ΛI(j)) fθI(i) fθI(j) = Sθq 1− eA 0 q(εI(i)ΛI(i)−εI(j)ΛI(j))e(2iπH(=(θ))+o(1))(εI(i)ΛI(i)−εI(j)ΛI(j)) fθI(i) fθI(j) . Remarking that the function ΛI(i) 7−→ e−2iπH(=(tq−θq))εI(i)ΛI(i) is 1-periodic in the variable ΛI(i), this prescribes a way to appropriately define the limit tq → θq. As in the remark following Theorem 4.2, it is actually sufficient to consider constants (fθI )I∈{I(i), i∈Kq}, from which it is possible to factorize the non-continuous part ΛI(i) 7−→ e−2iπH(=(θ))εI(i)ΛI(i) . Then, we will assume that there exist non-zero constants (fI)I∈{I(i), i∈Kq}, which verify, for any i ∈ Jp, with p ∈ Pq, fθI(i) = e−2iπH(=(θ))(εI(i)ΛI(i)−εI(min Jp)ΛI(min Jp))fI(i), where the constant fI(i) is a continuous function of θ ∈ C\R∗ (or at least in a neighboorhood of 0 in C \ R∗), with the additional condition that fI(i)(θ) −−−→ θ→0 f0 I(i), for any i ∈ Kq. In particular, we verify that fθI(min Jp) = fI(min Jp) = f0 I(min Jp) = 1, for any p ∈ Pq, as required by Theorem 4.2. This choice implies that, for any (i, j) ∈ J(2,�D ) p , with p ∈ Pq, ∆θ I(i)I(j)(ΛI(i),ΛI(j)) = Sθq 1− eA 0 q(εI(i)ΛI(i)−εI(j)ΛI(j))eo(1) f0I(i) f0I(j) 34 J. Avan, B. Billaud and G. Rollet −−−→ θ→0 S0 q 1− eA 0 q(εI(i)ΛI(i)−εI(j)ΛI(j)) f0I(i) f0I(j) = ∆0 I(i)I(j)(ΛI(i),ΛI(j)). If rq = 1, as far as ∆-coefficients are concerned, the proof of the theorem ends here. Assuming that rq ≥ 2, the limit for the other ∆-coefficients is not problematic, since ∆I(i)I(j) = Sq → Sθq = ∆θ I(i)I(j), for any (i, j) ∈ Jp × Jp′ , with (p, p′) ∈ P(2,<) q . The limit for d-coefficients is trivially deduced from above, for any (i, j) ∈ K(2,�D ) q . � It is therefore sufficient to consider the family of constants (Sq,Σq)q∈N∗s as general parameters characterized by the case iii above, cases i and ii being independent limits thereof. 5.2 Scaling of solutions Let us now examine the interpretation of the set-theoretical parameter identified with the specification of the ordered partition of any subset Kq into ∆-classes {Jp, p ∈ N∗r}. There exists another kind of continuity property of the solutions, which relies on the freedom allowed by the choice of the constants (fI)I∈{I(i), i∈Kq}. This result brings to light the fact that the solution built on a subset Kq such that rq = 1, i.e. Kq is reduced to a single ∆-class Jp, is the elementary solution, in the sense that it can generate any other solution by a limit process of an adequate re-scaling of the constants (fI)I∈{I(i), i∈Kq}. Let q ∈ N∗s such that rq ≥ 2, a solution (parametrized by the constants Sq 6= 0, Σq and (εI, fI)I∈{I(i), i∈Kq}) of system (S) restricted to the subset Kq and a parameter η > 0. Proposition 5.3. There exist a solution (parametrized by the constants Sq, Σq and (εI, f η I )I∈{I(i), i∈Kq}) of system (S) restricted to the subset Kq, such that rηq = 1 and that the solution with (fI)I∈{I(i), i∈Kq} is ��D -multiplicatively reducible to the limit of the solution with (fηI )I∈{I(i), i∈Kq}, up to a permutation of the indices in Kq, when η → 0+. Proof. We first construct a permutation σq: Kq −→ Kq, such that, after re-ordering, the subset Kq is an ordered partition as required by Theorem 3.1 in the case of rq = 1. We accordingly bring all free indices to the beginning of the subset. Introducing the ordered partition of the subset Kq in free subsets and ∆-classes Kq = ⋃ p∈Pq Jp = ⋃ p∈Pq lp⋃ l=0 I(p)l , this can be done by defining, for any p ∈ Pq, the permutation σq ∈ Sn, whose support is a subset of Kq, as σq(i) = i− ∑ p′∈Pq p′≤p np′ + np + ∑ p′∈Pq p′≤p n (p′) 0 − n(p) 0 , ∀ i ∈ I(p)0 ; and σq(i) = i+ ∑ p′∈Pq p′≥p n (p′) 0 − n(p) 0 , ∀ i ∈ Jp \ I(p)0 . Therefore, the subset Kq can be written as the following ordered partition Kσ q =  ⋃ p∈Pq σq(I (p) 0 ) ⋃ ⋃ p∈Pq σq(Jp \ I(p)0 )  , Classification of Non-Affine Non-Hecke Dynamical R-Matrices 35 where we recall that the exponent “σ” indicates that the permutation σq is applied. Moreover, the permutation σq respects the relation ��D , meaning that the pair of indices (i, j) ∈ K(2,<) q is a ��D -pair, if and only if the pair of indices (σq(i), σq(j)) ∈ (Kσ q )2 is a ��D -pair. Theorems 4.2 and 4.3 insure the existence of a full solution (parametrized by the constants Sq, Σq and (εI, f η I )I∈{I(i), i∈Kσq }, and the ��D -closed ��D -multiplicative 2-form (gηij)(i,j)∈(Kσq )(2,�D,<) ) of system (S) restricted to the subset Kq, where the constants (fηI )I∈{I(i), i∈Kσq } are defined, for any p ∈ Pq, as fηI(σq(i)) = η−p+pq+1fI(i), ∀ i ∈ Jp. In particular, we verify that fηI(minKq) = fηI(σq(minKq)) = fηI(σq(min Jpq+1)) = fI(min Jpq+1) = 1, as required by Theorem 4.2. Therefore, by construction, we have, for any p ∈ Pq fηI(σq(j)) fηI(σq(i)) = fI(j) fI(i) , ∀ (i, j) ∈ J(2,�D ) p ; and, for any (p, p′) ∈ P(2,<) q , fηI(σq(i)) fηI(σq(j)) = ηp ′−p fI(i) fI(j) −−−−→ η→0+ 0, ∀ (i, j) ∈ Jp × Jp′ . This directly implies the expected result, which is that, for any p ∈ Pq ∆η I(σq(i))I(σq(j)) = ∆I(i)I(j), ∀ (i, j) ∈ J(2,�D ) p ; and, for any (p, p′) ∈ P(2,<) q ∆η I(σq(i))I(σq(j)) −−−−→η→0+ Sq = ∆I(i)I(j) ⇔ ∆η I(σq(j))I(σq(i)) −−−−→η→0+ 0 = ∆I(j)I(i), ∀ (i, j) ∈ Jp × Jp′ . The limit for d-coefficients is trivially deduced from above, up to the multiplication by the ��D -closed ��D -multiplicative 2-form (gij) (i,j)∈K(2,�D,<) q , defined as, for any p ∈ Pq gij = 1, ∀ (i, j) ∈ J(2,�D ) p ; and, for any (p, p′) ∈ P(2,<) q gij = √ Σq Bq − Sq , gji √ Σq Bq , ∀ (i, j) ∈ Jp × Jp′ . � Remark 5.2. The construction of the scaled R-matrix in terms of the non-zero constants (fI)I∈{I(i), i∈Kq} guarantees that the scaled R-matrix with multiple ∆-classes still preserves the fundamental feature that a D-class is inside a single ∆-class, as shown in Proposition 3.5. We conclude that the solution built on any ordered partition of any subset Kq into ∆-classes {Jp, p ∈ N∗r} is obtained as the result of scaling procedure applied to a solution built on a single ∆-class Kq = Jp. The single ∆-class solution is then generic. 36 J. Avan, B. Billaud and G. Rollet 5.3 Re-parametrization of variables (λk)k∈N∗ n Here we propose an interpretation of the parameters (fI)I∈{I(i), i∈N∗n}. Proposition 5.4. Let a R-matrix, solution of DQYBE, be parametrized (among others) by the set of constants (fI)I∈{I(i), i∈N∗n}. Then, there exists a re-parametrization of the dynamical variable λk, for any k ∈ N∗n, which eliminates this dependence. Proof. From Theorem 4.2 (when there exists q ∈ N∗s such that the R-matrix, restricted to the subset Kq, has a trigonometric behavior, i.e. when Sq 6= 0) or Theorem 4.3 (when there exists q ∈ N∗s such that the R-matrix, restricted to the subset Kq, has a rational behavior, i.e. when Sq = 0), it is manifest that these parameters can be respectively re-absorbed in a re-definition of the dynamical variables (λk)k∈N∗n as λi −→ λi + εI(i) |I(i)| logS(fI(i)) Aq or λi −→ λi + εI(i) |I(i)| fI(i), ∀ i ∈ N∗n. Following the same argumentation as in the remark of Theorem 4.2 concerning the definition of the quantity Aq, this re-parametrization of the dynamical variables (λk)k∈N∗n (when the R- matrix has a trigonometric behavior) is indeed independent from the choice of the determination of the logarithm function, and of the choice of log(fI(i)), when fI(i) ∈ R∗−. This justifies that the constants (fI)I∈{I(i), i∈N∗n} should be advantageously seen as belonging to the Riemann surface S, as well as the use of function logS. � This re-parametrization of the dynamical variables (λk)k∈N∗n is the only one under which DQYBE is form-invariant, since it must preserve the translation λi −→ λi + 1, for any i ∈ N∗n. This is the reason why it will be refereed as the canonical parametrization of the dynamical variables (λk)k∈N∗n . Moreover, the family of signs (εI)I∈{I(i), i∈N∗n} cannot be re-absorbed in such a way, and represents a set of genuine relevant parameters of a generic solution of DQYBE, to be interpreted in the next subsection. To summarize, any R-matrix, solution of DQYBE, characterized by an ordered partition of the ∆-class Jp into D-classes {I(i), i ∈ Jp}, is built by juxtaposition following Theorem 3.1 and Proposition 4.5, and by ��D -multiplicative covariance following Propositions 4.2, 4.3 and 4.4, of solutions obtained by • limit following Propositions 5.1 and 5.2, • scaling following Proposition 5.3, • re-parametrization of the dynamical variables following Proposition 5.4; of a solution (parametrized by the non-zero constants Sp and Σp, the family of signs (εI)I∈{I(i), i∈Jp}, and the ��D -multiplicative 2-form (gij) (i,j)∈J(2,�D,<) p ) of DQYBE on a single ∆- class Jp. 5.4 Commuting operators The form of a generic R-matrix, solution of DQYBE given by Theorem 5.1, allows to immediately bring to light a set of operators which commute with the R-matrix. Proposition 5.5. For any R-matrix, solution of DQYBE, the operator R0 = s∑ q=1 ∑ p∈Pq ∑ i∈I(p)0 R (q,p) i = ∑ i∈I0 ∆iie (n) ii ⊗ e (n) ii Classification of Non-Affine Non-Hecke Dynamical R-Matrices 37 together with the family of operators (R (q,p) I(i) ) i∈Jp\I(p)0 , p∈Pq |q∈N∗s build a set of mutually commuting operators, commuting with the R-matrix. Proof. 1. Mutual commutation. For any q ∈ N∗s, for any p ∈ Pq, we recall the formula R (q,p) I(i) = ∆I(i) |I(i)| ∑ j,j′∈I(i) e (n) jj′ ⊗ e (n) j′j , ∀ i ∈ Jp. For any (i, j) ∈ (Jp \ I(p)0 )× (Jp′ \ I (p′) 0 ) | (i, j) ∈ (N∗n)2, with q, q′ ∈ N∗s and (p, p′) ∈ Pq × Pq′ , it is straightforward to check that[ R0, R (q,p) I(i) ] = [ R (q,p) I(i) , R (q′,p′) I(j) ] = 0. This relies on the fact that the indices appearing in the sum defining the operator R0 or R (q′,p′) I(j) do not appear in the sum defining the operator R (q,p) I(i) , because of the partitioning of the set of indices N∗n. 2. Commutation with the R-matrix. This implies in particular that any operators of the family (R0, (R (q,p) I(i) ) i∈Jp\I(p)0 , p∈Pq |q∈N∗s ) commute with their sum, being the operator R0 + s∑ q=1 ∑ p∈Pq ∑ i∈Jp\I(p)0 R (q,p) I(i) = s∑ q=1 ∑ p∈Pq ∑ i∈Jp R (q,p) I(i) = ∑ i∈N∗n R (q,p) I(i) . Remember now that a R-matrix, solution of DQYBE, is essentially the sum of three kinds of terms, which are the previous sum (containing all terms ∆jj′e (n) jj′ ⊗ e (n) j′j , when jDj′), terms such as ∆jj′e (n) jj′ ⊗ e (n) j′j and such as djj′e (n) jj ⊗ e (n) jj′ , when j��D j′. It remains to check the commutativity of any operators of the family (R0, (R (q,p) I(i) ) i∈Jp\I(p)0 , p∈Pq |q∈N∗s ) with the second and third kind of terms appearing in the expression of a generic solution of DQYBE. To this end, let i ∈ N∗n, j, j′ ∈ I(i) and (k, k′) ∈ N∗(2,�D ) n . Here, for convenience, we allow for once to have an equality between two indices distinctly labelled. However, there are compatibility conditions to fulfil. For example, we can have j = k, but not at the same time j′ 6= k′, otherwise kDk′. Such considerations give directly that( e (n) jj′ ⊗ e (n) j′j )( e (n) kk′ ⊗ e (n) k′k ) = ( e (n) kk′ ⊗ e (n) k′k )( e (n) jj′ ⊗ e (n) j′j ) = 0 and ( e (n) jj′ ⊗ e (n) j′j )( e (n) kk ⊗ e (n) k′k′ ) = ( e (n) kk ⊗ e (n) k′k′ )( e (n) jj′ ⊗ e (n) j′j ) = 0, which yields, for any i ∈ Jp \ I(p)0 , with q ∈ N∗s and p ∈ Pq [R,R0] = [ R,R (q,p) I(i) ] = 0. � 5.5 (weak) Hecke and non-Hecke R-matrices We recall that we have dropped any Hecke or weak Hecke condition in our derivation Nevertheless such conditions will be shown to arise in connection with the choice of the family of signs (εI)I∈{I(i), i∈N∗n}, parametrizing any R-matrix. Let us first give the definition for the (weak) Hecke condition following [14]. 38 J. Avan, B. Billaud and G. Rollet Definition 5.1 ((weak) Hecke condition). i. A R-matrix satisfies the Hecke condition with parameters %, κ ∈ C∗, such that % 6= −κ, if the eigenvalues of the permuted R-matrix Ř = PR, P being the permutation operator of spaces V1 and V2, are % on the one-dimensional vector space Vii = Ce(n) i ⊗ e (n) i , for any in- dex i ∈ N∗n, and %, −κ on the two-dimensional vector space Vij = Ce(n) i ⊗e (n) j ⊕Ce (n) j ⊗e (n) i , for any pair of indices (i, j) ∈ N∗(2,<) n , where (e (n) i )i∈N∗n is the canonical basis of the vector space V = Cn. ii. A R-matrix satisfies the weak Hecke condition with parameters %, κ ∈ C∗, such that % 6= −κ, if the minimal polynomial of the permuted matrix Ř is µŘ(X) = (X − %)(X +κ). Remark 5.3. A R-matrix, for which the set of indices N∗n is reduced to a single D-class is not strictly of (weak) Hecke-type with parameters %, κ ∈ C∗, such that % 6= −κ, since in this case µŘ(X) = X −∆I(1). By language abuse, it can be considered as a degenerate weak Hecke-type R-matrix with parameters %, κ ∈ C∗, such that % = −κ = ∆I(1). In the following, we will assume that the set of indices N∗n is not reduced to a single D-class, unless otherwise stated. The classification of Hecke-type solutions of DQYBE is well known [13], whereas the clas- sification of weak Hecke-type solutions of DQYBE remains unknown. Reference [14] presents two fundamental theorems treating separately trigonometric and rational cases. They stress that, according to the value of the parameter κ and up to gauge transformations and to addi- tional trivial transformations, such as scalar multiplication, or global linear re-parametrization of the dynamical variable λ, there essentially exist two distinct types of R-matrices satisfying the Hecke condition, the so-called basic trigonometric (% = 1 and κ 6= 1) and rational (% = κ = 1) Hecke-type solutions. Let r ∈ N∗n and {Jp, p ∈ N∗r} be an ordered partition of the set of indices N∗n into r ∆-classes. Using the canonical parametrization of the dynamical variables (λk)k∈N∗n of Proposition 5.4, by simple identification, these basic Hecke-type solutions are particular examples of solutions of DQYBE given by Theorems 4.2, 4.3, 4.4 and 5.1, as described in the following definitions. Definition 5.2 (basic trigonometric and rational Hecke-type behavior). i. The basic trigonometric Hecke-type R-matrix, solution of DQYBE, with parameter κ ∈ C∗ \ {1}, is, up to trivial transformations, a R-matrix, parametrized by • I0 = N∗n, i.e. N∗n does not contain any D-classes; • s = 1, i.e. N∗n = K1 = r⋃ p=1 Jp; • the non-zero constants S1 = κ− 1, Σ1 = κ and T1 = 1 κ ; • the signs εI(i) = εi = −1, for any i ∈ N∗n; • the multiplicative 2-form gij = −1, for any (i, j) ∈ N∗(2,�D ) n . ii. The basic rational Hecke-type R-matrix, solution of DQYBE, is, up to trivial transforma- tions, a R-matrix, parametrized by • I0 = N∗n, i.e. N∗n does not contain any D-classes; • s = r, i.e. Kq = Jq, for any q ∈ N∗s; • the constants Sq = 0 and the non-zero-constants Σq = 1, for any q ∈ N∗s; • the non-zero constants Σqq′ = 1, for any (q, q′) ∈ N∗(2,<) s ; Classification of Non-Affine Non-Hecke Dynamical R-Matrices 39 • the signs εI(i) = εi = −1, for any i ∈ N∗n; • the multiplicative 2-form gij = −1, for any (i, j) ∈ N∗(2,�D ) n . Remark 5.4. The trivial manipulations to get Definition 5.2 from [14] are a multiplication of the R-matrix by −1, and a global re-parametrization of the dynamical variables (λk)k∈N∗n as λ→ −λ, the rest relies on a simple identification. To be as exhaustive as possible, we propose an alternative formulation of the classification theorems of Hecke-type solutions of DQYBE [14], as well as a classification of the weak Hecke- type solutions of DQYBE. Proposition 5.6 (Hecke-type R-matrices). Any Hecke-type R-matrix, solution of DQYBE, with parameters %, κ ∈ C∗, such that % 6= −κ, is twist-gauge reducible to a basic trigonometric or rational Hecke-type R-matrix. Proof. This is an obvious corollary of the classification theorems of Hecke-type solutions of DQYBE [13], and of the two Definitions 5.2. However, we can explicit why no other R- matrix does satisfy the Hecke condition. Let R be a matrix, solution of DQYBE, built on the ordered partition {Jp, p ∈ N∗r} of the set of indices N∗n, which additionally satisfies the Hecke condition. The proof of the proposition relies on the fact that, by direct calculation, the permuted matrix Ř is expressed as follows Ř = n∑ i,j=1 ∆jie (n) ii ⊗ e (n) jj + n∑ i 6=j=1 djie (n) ij ⊗ e (n) ji , from which we deduce the characteristic polynomial of its restriction to the subspace Vij , for any pair of indices (i, j) ∈ N∗(2,<) n , as Pij(X) = ∣∣∣∣∆ji −X dji dij ∆ij −X ∣∣∣∣ = X2 − SijX − Σij . By projecting the permuted matrix Ř in the subspace Vii, we deduce immediately that ∆I(i) = ∆ii = %, for any i ∈ N∗n. If I0 ( N∗n, then there exists a pair of indices (i, j) ∈ N∗(2,<) n such that iDj. From Propo- sition 4.1, we have ∆ij = ∆ji = ∆I(i), i.e. Sij = 2∆I(i) and Σij = −∆2 I(i), yielding that Pij(X) = (X − ∆I(i)) 2. This leads to a contradiction. The matrix Ř does not indeed satis- fy the Hecke condition, since it has a single eigenvalue ∆I(i) on the subspace Vij . One must then assume that I0 = N∗n and the set of indices only contains free indices. For any pair of indices (i, j) ∈ N∗(2,<) n , Definition 5.1 directly implies that % and −κ are root of the polynomial Pij , i.e. Sij = %− κ and Σij = %κ. This means in particular that Sq = %− κ, Σq = %κ and Tq = κ % , for any q ∈ N∗s, as well as that Σqq′ = %κ, for any (q, q′) ∈ N∗(2,<) s , if s ≥ 2. But, these equalities are more constraining than that • If % 6= κ, then s = 1. Otherwise, if s ≥ 2, let (q, q′) ∈ N∗(2,<) s and a pair of indices (i, j) ∈ Kq ×Kq′ . Then, by construction, Sij = %− κ = 0, which leads to a contradiction. Hence any trigonometric solution has a single subset K1. • If % = κ, then rq = 1, for any q ∈ N∗s. Otherwise, let q ∈ N∗s, such that rq ≥ 2. By Theorem 4.4, this leads to a contradiction, since Sq = 0. Hence any rational solution is a decoupled R-matrix, for which any subset Kq is reduced to a single ∆-class. 40 J. Avan, B. Billaud and G. Rollet Finally, in both cases, from expressions for the non-zero constants (∆ii)i∈N∗n in Theorem 4.1, this additionally imposes that εI(i) = εi = 1, for any i ∈ N∗n. Since, for any R-matrix satisfying the Hecke condition with parameters %, κ ∈ C∗, the matrix 1 %R satisfies the Hecke condition with parameters 1 and κ′ = κ % , it is sufficient to assume that % = 1 without loss of generality. This gives the expected results up to the trivial manipulations presented in the remark of Definition 5.2, which particularly make the re-parametrizations Sq → S′q = −Sq = κ′ − 1 and Tq → T ′q = 1 Tq = 1 κ′ , for any q ∈ N∗s and εi → ε′i = −1, for any i ∈ N∗n. Finally let us expose why each non-zero function of the family (gij) (i,j)∈N∗(2,�D) n can be brought to −1. For any Hecke-type solution, I0 = N∗n, hence any pair of indices (i, j) ∈ (N∗n)2 is a ��D -pair, so that, according to Propositions 4.2, 4.3 and 4.4, the family of functions (−gij) (i,j)∈N∗(2,�D) n is a closed multiplicative 2-form, then exact under appropriate assumptions, multiplicative 2-form. This means that they are of gauge form, and can be factorized out thanks to Proposition 2.1, finally leaving the wanted factor −1 in front of any d-coefficient. It remains now to apply the trivial manipulations of the remark of Definition 5.2. � Remark 5.5. From Proposition 4.5, if s ≥ 2, a basic rational Hecke-type solution is by con- struction a decoupled R-matrix, whereas a basic trigonometric Hecke-type solution never is. These two different behaviors are unified as soon as the Hecke condition is dropped. In the same spirit, Proposition 5.1 generalizes the well-known property that basic rational Hecke-type solutions can be obtained as limits of basic trigonometric Hecke-type solutions of parameter κ ∈ C∗ \ {1}, when κ→ 1. By analogy with the terminology used for Hecke-type solutions, we will introduce the no- tion of basic trigonometric or rational weak Hecke-type solutions of DQYBE as follows. Let us particularly insist on the fact that, as we will see explicitly later, unlike the Hecke-type con- dition, the weak Hecke-type condition allow ∆-classes and do not constrain the choice of signs (εI)I∈{I(i), i∈N∗n}. Definition 5.3 (basic trigonometric and rational weak Hecke-type behavior). i. The basic trigonometric weak Hecke-type R-matrix, solution of DQYBE, with parameter κ ∈ C∗ \ {1}, is, up to trivial transformations, a R-matrix, parametrized by • s = 1, i.e. N∗n = K1 = r⋃ p=1 Jp; • the non-zero constants S1 = κ− 1, Σ1 = κ and T1 = 1 κ ; • the multiplicative 2-form gij = −1, for any (i, j) ∈ N∗(2,�D ) n . ii. The basic rational weak Hecke-type R-matrix, solution of DQYBE, is, up to trivial trans- formations, a R-matrix, parametrized by • s = r, i.e. Kq = Jq, for any q ∈ N∗s; • the constants Sq = 0 and the non-zero-constants Σq = 1, for any q ∈ N∗s; • the non-zero constants Σqq′ = 1, for any (q, q′) ∈ N∗(2,<) s ; • the multiplicative 2-form gij = −1, for any (i, j) ∈ N∗(2,�D ) n . In both cases, the ordered partition of set of indices N∗n remain free, and the family of signs (εI)I∈{I(i), i∈N∗n} is arbitrary. Proposition 5.7 (weak Hecke-type R-matrices). Any R-matrix, solution of DQYBE, which satisfies the weak Hecke condition with parameters %, κ ∈ C∗, such that % 6= −κ, is ��D - multiplicatively reducible to a basic trigonometric or rational weak Hecke-type R-matrix. Classification of Non-Affine Non-Hecke Dynamical R-Matrices 41 Proof. Thanks to the zero-weight condition, the permuted matrix Ř is block diagonal, up to a permutation in Sn2 , where the blocks are of the form ∆I(i) (in the subspace Vii) or ( ∆ji dji dij ∆ij ) . By projecting the permuted matrix Ř in the subspace Vii, we deduce immediately either ∆I(i) = % or ∆I(i) = −κ, for any i ∈ N∗n. Assuming that N∗n is not reduced to a single D-class, the restriction of the matrix Ř to the subspace Vij has to satisfy the weak Hecke condition, for any pair of indices (i, j) ∈ N∗(2,<) n , i.e. µij | µŘ, where µij is the minimal polynomial of the restriction of the matrix Ř. As above, this is trivially satisfied for any pair of indices (i, j) ∈ N∗(2,<) n , such that iDj. Let then now a ��D -pair of indices (i, j) ∈ N∗(2,�D ,<) n . The rest of the proof is almost identical to the proof of Proposition 5.6, since we have that µij = Pij = µŘ. The major differences are that the signs (εI)I∈{I(i), i∈N∗n} are no longer constrained to be all equal (from Theorems 4.2 and 4.3, we particularly deduce that εi = 1, if ∆ii = %, and εi = −1, if ∆ii = −κ), and that we have to use the ��D -multiplicative covariance instead of the twist covariance, when the set of indices N∗n contains D-classes, which was excluded in the case of Hecke-type R-matrices. Let us note that the ��D -multiplicative covariance does not affect the minimal polynomial µŘ, since we have just proved that µŘ = Pij , for any pair of indices (i, j) ∈ N∗(2,<) n , where the polynomial Pij is obviously invariant under such transformation. This occurrence only arises, when n ≥ 3, hence it could not arise in the classification of R-matrices, solution of Gl2(C), as seen for trigonometric R-matrices in [21]. � We now come to the main concluding statement. Theorem 5.2 (decoupling theorem). Any R-matrix, solution of DQYBE, parametrized (among others) by the family of constants (Sq,Σq)q∈N∗s , is ��D -multiplicatively reducible to a decoupled R-matrix, whose constituting blocks are either D-classes or satisfy the weak Hecke condition with the family of parameters (%q, κq)q∈N∗s , such that %q 6= −κq and Sq = %q − κq and Σq = %qκq, ∀ q ∈ N∗s. Proof. This is a corollary of Propositions 4.5 and 5.7 and Theorems 3.1, 4.1, 4.2, 4.3, 4.4 and 5.1. � Beforing concluding this article, Fig. 1 gives an example of a non-Hecke R-matrix, solution of Gl4(C)-DQYBE with rational behavior. We have chosen the set of indices N∗4 = K1 to be a single ∆-class J1, where the free subset is I0 = I(1) 0 = {1, 2} and where the remaining indices I(1) 1 = {3, 4} form a D-class. For any i ∈ I0 or for any i ∈ I(1) 1 , we have respectively denoted the quantities depending on I(i), such as the sign εI(i) or the variable ΛI(i), by “i” or by “ ” instead of the previous notation “I(i)”. We have also dropped the index “1” for the constants S1 and Σ1. We have chosen the signs ε1 = ε2 = 1 and ε = −1, and have fixed the non-zero constants f1, f2 and f to 0, and the ��D -multiplicative 2-form to 1. The R-matrix we present is non-Hecke thanks to the presence of different signs and the D-class I(1) 1 , but satisfies the weak Hecke condition, because it is not decoupled. 6 Conclusion We have proceeded to the exhaustive classification of the non-affine “non-Hecke”-type quan- tum Gln(C) dynamical R-matrices obeying DQYBE. In particular, we have succeeded to fully characterize its space of moduli, and prove that weak Hecke-type R-matrices are the elemen- tary constituting blocks of non-Hecke-type R-matrices. This classification then brings to light 42 J. Avan, B. Billaud and G. Rollet R =                                     1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 − 1 λ 1 − λ 2 0 0 1 λ 1 − λ 2 0 0 0 0 0 0 0 0 0 0 0 0 0 1 − 1 λ 1 + λ 0 0 0 0 0 1 λ 1 + λ 0 0 0 0 0 0 0 0 0 0 1 − 1 λ 1 + λ 0 0 0 0 0 0 0 0 1 λ 1 + λ 0 0 0 0 − 1 λ 1 − λ 2 0 0 1 + 1 λ 1 − λ 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 − 1 λ 2 + λ 0 0 1 λ 2 + λ 0 0 0 0 0 0 0 0 0 0 0 0 0 1 − 1 λ 2 + λ 0 0 0 0 0 1 λ 2 + λ 0 0 0 0 − 1 λ 1 + λ 0 0 0 0 0 1 + 1 λ 1 + λ 0 0 0 0 0 0 0 0 0 0 0 0 0 − 1 λ 2 + λ 0 0 1 + 1 λ 2 + λ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 − 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 − 1 0 0 0 0 − 1 λ 1 + λ 0 0 0 0 0 0 0 0 1 + 1 λ 1 + λ 0 0 0 0 0 0 0 0 0 0 − 1 λ 2 + λ 0 0 0 0 0 1 + 1 λ 2 + λ 0 0 0 0 0 0 0 0 0 0 0 0 0 − 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 − 1                                      Figure 1. Example of a non-Hecke R-matrix, solution of Gl4(C)-DQYBE. a wide range of new solutions to this equation, while the Hecke-type solutions appear as a very particular type of solutions. As a matter of fact, the parametrization of a general solution of DQYBE involves a large number of objects of different mathematical natures, which is drasti- cally restricted when the Hecke condition is considered. These results may pave the way for the classification of affine non-Hecke-type quantum dynamical R-matrices obeying DQYBE, or less ambitiously may be a first step for the under- Classification of Non-Affine Non-Hecke Dynamical R-Matrices 43 standing of the Baxterization of non-affine non-Hecke-type quantum dynamical R-matrices, whose general case still remains an open problem nowadays. Occurrence of weak Hecke building blocks, for which Baxterization procedure is known [2], at least in the non-dynamical case, allows to be quite hopeful in this respect. Moreover, the non-Hecke-type solutions of DQYBE are interesting by themselves. In recent developments of researches on the second Poisson structure of Calogero models emerge non- Hecke-type solutions of dynamical classical Yang–Baxter equation [3], such as r =  0 0 0 0 0 w1 λ1−λ2 2 λ1−λ2 0 0 2 λ2−λ1 w2 λ2−λ1 0 0 0 0 0  , where w1, w2 ∈ C∗. More precisely, this matrix is the solution for the matrix a, occurring in a general quadratic Poisson bracket algebra {l1, l2} = al1l2 + l1bl2 + l2cl1 + l1l2d, and can be obviously obtained as a semi-classical limit of a non-Hecke solution of DQYBE. Appendix Here we give an exhaustive lists of definitions, theorems, propositions and lemmas, along with their names or some keywords. List of definitions Definition 4.1 Multiplicative 2-forms . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18 Definition 5.1 (weak) Hecke condition. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .38 Definitions 5.2 and 5.3 Basic trigonometric and rational (weak) Hecke-type behavior . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38 & 40 List of theorems Theorem 3.1 ∆-incidence matrices . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13 Theorem 4.1 inside a set Kq of ∆-classes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21 Theorems 4.2 and 4.3 Trigonometric and rational behavior . . . . . . . . . . . . . . . . . . . . . 24 & 27 Theorem 4.4 Trigonometric behavior . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28 Theorem 5.1 General R-matrices . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30 Theorem 5.2 Decoupling theorem . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41 List of propositions Propositions 2.1, 2.2 and 2.3 Twist covariance, decoupled and contracted R-matrices . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .5 & 6 & 7 Proposition 3.1 and 3.2 d-indices . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7 Proposition 3.3, 3.4 and 3.5 ∆-indices . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8 Proposition 3.6 and 3.7 (Reduced) ∆-incidence matrix . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9 44 J. Avan, B. Billaud and G. Rollet Proposition 3.8 Triangularity . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10 Proposition 3.9, 3.10 and 3.11 Comparability of ∆-classes . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10 Proposition 3.12 Block upper-triangularity . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11 Proposition 4.1 Inside a D-class . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15 Propositions 4.2, 4.3 and 4.4 (��D -)multiplicative covariance . . . . . . . . . . . . . . . . . . . . . 18 & 19 Proposition 4.5 Decoupling proposition . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20 Propositions 5.1 and 5.2 Continuity propositions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31 & 32 Proposition 5.3 Scaling proposition . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34 Proposition 5.4 Re-parametrization of variables . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36 Proposition 5.5 Commuting operators . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .36 Proposition 5.6 and 5.7 (Weak) Hecke-type R-matrices . . . . . . . . . . . . . . . . . . . . . . . . . 39 & 40 List of lemmas Lemma 3.1 Order of ∆-classes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10 Lemmas 4.1 and 4.2 Dependences and forms of the ∆-coefficients . . . . . . . . . . . . . . . . . . . . 16 Lemmas 4.3 and 4.4 Dependences and forms of the d-coefficients . . . . . . . . . . . . . . . 16 & 17 Lemmas 4.5 and 4.6 Multiplicative and additive shifts . . . . . . . . . . . . . . . . . . . . . . . . . . 23 & 24 Acknowledgements This work was supported by CNRS, Université de Cergy Pontoise, and ANR Project DIADEMS (Programme Blanc ANR SIMI 1 2010-BLAN-0120-02). We wish to thank Dr. Thierry Gobron for interesting discussions and comments on the structure of ∆-incidence matrices, and Professor Laszlo Fehér for pointing out to us reference [5]. Finally we thank the anonymous referees for their careful study and their relevant contributions to improve the paper. References [1] Avan J., Babelon O., Billey E., The Gervais–Neveu–Felder equation and the quantum Calogero–Moser systems, Comm. Math. Phys. 178 (1996), 281–299, hep-th/9505091. [2] Avan J., Kulish P.P., Rollet G., Reflection K-matrices related to Temperley–Lieb R-matrices, Theoret. and Math. Phys. 169 (2011), 1530–1538, arXiv:1012.3012. [3] Avan J., Ragoucy E., Rational Calogero–Moser model: explicit forms and r-matrix structure of the second Poisson structure, arXiv:1207.5368. [4] Avan J., Talon M., Classical R-matrix structure for the Calogero model, Phys. Lett. B 303 (1993), 33–37, hep-th/9210128. [5] Balog J., Dabrowski L., Fehér L., Classical r-matrix and exchange algebra in WZNW and Toda theories, Phys. Lett. B 244 (1990), 227–234. [6] Balog J., Fehér L., Palla L., Chiral extensions of the WZNW phase space, Poisson–Lie symmetries and groupoids, Nuclear Phys. B 568 (2000), 503–542, hep-th/9910046. [7] Bernard D., On the Wess–Zumino–Witten models on Riemann surfaces, Nuclear Phys. B 309 (1988), 145– 174. [8] Bernard D., On the Wess–Zumino–Witten models on the torus, Nuclear Phys. B 303 (1988), 77–93. [9] Calogero F., Solution of a three-body problem in one dimension, J. Math. Phys. 10 (1969), 2191–2196. [10] Donin J., Mudrov A., Dynamical Yang–Baxter equation and quantum vector bundles, Comm. Math. Phys. 254 (2005), 719–760, math.QA/0306028. http://dx.doi.org/10.1007/BF02099449 http://arxiv.org/abs/hep-th/9505091 http://dx.doi.org/10.1007/s11232-011-0130-y http://dx.doi.org/10.1007/s11232-011-0130-y http://arxiv.org/abs/1012.3012 http://arxiv.org/abs/1207.5368 http://dx.doi.org/10.1016/0370-2693(93)90039-K http://arxiv.org/abs/hep-th/9210128 http://dx.doi.org/10.1016/0370-2693(90)90060-J http://dx.doi.org/10.1016/S0550-3213(99)00738-5 http://arxiv.org/abs/hep-th/9910046 http://dx.doi.org/10.1016/0550-3213(88)90236-2 http://dx.doi.org/10.1016/0550-3213(88)90217-9 http://dx.doi.org/10.1063/1.1664820 http://dx.doi.org/10.1007/s00220-004-1247-8 http://arxiv.org/abs/math.QA/0306028 Classification of Non-Affine Non-Hecke Dynamical R-Matrices 45 [11] Enriquez B., Etingof P., Quantization of classical dynamical r-matrices with nonabelian base, math.QA/0311224. [12] Etingof P., On the dynamical Yang–Baxter equation, in Proceedings of the International Congress of Mathe- maticians, Vol. II (Beijing, 2002), Higher Ed. Press, Beijing, 2002, 555–570, math.QA/0207008. [13] Etingof P., Schiffmann O., Lectures on the dynamical Yang–Baxter equations, in Quantum Groups and Lie Theory (Durham, 1999), London Math. Soc. Lecture Note Ser., Vol. 290, Cambridge Univ. Press, Cambridge, 2001, 89–129, math.QA/9908064. [14] Etingof P., Varchenko A., Solutions of the quantum dynamical Yang–Baxter equation and dynamical quan- tum groups, Comm. Math. Phys. 196 (1998), 591–640, math.QA/9801135. [15] Felder G., Conformal field theory and integrable systems associated to elliptic curves, in Proceedings of the International Congress of Mathematicians, Vol. 1, 2 (Zürich, 1994), Birkhäuser, Basel, 1995, 1247–1255, hep-th/9407154. [16] Felder G., Elliptic quantum groups, in XIth International Congress of Mathematical Physics (Paris, 1994), Int. Press, Cambridge, MA, 1995, 211–218, hep-th/9412207. [17] Furlan P., Hadjiivanov L.K., Isaev A.P., Ogievetsky O.V., Pyatov P.N., Todorov I.T., Quantum matrix algebra for the SU(n) WZNW model, J. Phys. A: Math. Gen. 36 (2003), 5497–5530, hep-th/0003210. [18] Gervais J.L., Neveu A., Novel triangle relation and absence of tachyons in Liouville string field theory, Nuclear Phys. B 238 (1984), 125–141. [19] Hadjiivanov L.K., Stanev Y.S., Todorov I.T., Regular basis and R-matrices for the ŝu(n)k Knizhnik– Zamolodchikov equation, Lett. Math. Phys. 54 (2000), 137–155, hep-th/0007187. [20] Isaev A.P., Twisted Yang–Baxter equations for linear quantum (super)groups, J. Phys. A: Math. Gen. 29 (1996), 6903–6910, q-alg/95110006. [21] Ju G., Luo X., Wang S., Wu K., A free-fermion type solution of quantum dynamical Yang–Baxter equation, Commun. Theor. Phys. 32 (1999), 557–562. [22] Knizhnik V.G., Zamolodchikov A.B., Current algebra and Wess–Zumino model in two dimensions, Nuclear Phys. B 247 (1984), 83–103. [23] Moser J., Three integrable Hamiltonian systems connected with isospectral deformations, Adv. Math. 16 (1975), 197–220. [24] Ruijsenaars S.N.M., Schneider H., A new class of integrable systems and its relation to solitons, Ann. Physics 170 (1986), 370–405. [25] Xu P., Quantum dynamical Yang–Baxter equation over a nonabelian base, Comm. Math. Phys. 226 (2002), 475–495, math.QA/0104071. http://arxiv.org/abs/math.QA/0311224 http://arxiv.org/abs/math.QA/0207008 http://arxiv.org/abs/math.QA/9908064 http://dx.doi.org/10.1007/s002200050437 http://arxiv.org/abs/math.QA/9801135 http://arxiv.org/abs/hep-th/9407154 http://arxiv.org/abs/hep-th/9412207 http://dx.doi.org/10.1088/0305-4470/36/20/310 http://arxiv.org/abs/hep-th/0003210 http://dx.doi.org/10.1016/0550-3213(84)90469-3 http://dx.doi.org/10.1023/A:1011025515626 http://arxiv.org/abs/hep-th/0007187 http://dx.doi.org/10.1088/0305-4470/29/21/021 http://arxiv.org/abs/q-alg/95110006 http://dx.doi.org/10.1016/0550-3213(84)90374-2 http://dx.doi.org/10.1016/0550-3213(84)90374-2 http://dx.doi.org/10.1016/0001-8708(75)90151-6 http://dx.doi.org/10.1016/0003-4916(86)90097-7 http://dx.doi.org/10.1007/s002200200621 http://arxiv.org/abs/math.QA/0104071 1 Introduction 2 Preparatory material 3 The -incidence matrix and equivalence classes 3.1 d-indices 3.2 -indices 3.3 (Reduced) -incidence matrix 3.4 Classification 4 Resolution 4.1 Preliminaries 4.2 Decoupling procedure 4.3 Sum and determinant 4.4 Inside a -class 4.5 Inside a subset Kq 5 General solution and structure of the set of solutions 5.1 Continuity properties of the solutions with respect to the constants (Sq,q)qN*s 5.2 Scaling of solutions 5.3 Re-parametrization of variables (k)kN*n 5.4 Commuting operators 5.5 (weak) Hecke and non-Hecke R-matrices 6 Conclusion References