Arxiv:2012.15778V2 [Math.RT] 18 Jan 2021 4.1

METAPLECTIC IWAHORI WHITTAKER FUNCTIONS AND SUPERSYMMETRIC LATTICE MODELS

BEN BRUBAKER, VALENTIN BUCIUMAS, DANIEL BUMP, AND HENRIK P. A. GUSTAFSSON

Abstract. In this paper we consider Iwahori Whittaker functions on n-fold metaplectic covers Ge of G(F ) with G a split reductive group over a non-archimedean local ﬁeld F . For every element φ of a basis of Iwahori Whittaker functions, and for every g ∈ Ge, we evaluate φ(g) by recurrence relations over the Weyl group using “vector Demazure-Whittaker operators.” Specializing to the case of G = GLr, we exhibit a solvable lattice model whose partition function equals φ(g). These models are of a new type associated with the quantum aﬃne super group Uq(glb(r|n)). The recurrence relations on the representation theory side then correspond to solutions to Yang-Baxter equations for the lattice models. Remarkably, there is a bijection between the boundary data specifying the partition function and the data determining all values of the Whittaker functions.

Contents 1. Introduction2 2. Supersymmetric lattice models7 2.1. The monochrome model7 2.2. Equivalent model via fusion 11 2.3. The partition function for the ground state system 14 2.4. The R-matrix and the Yang-Baxter equation 16 2.5. Recurrence relations for the partition functions 22 3. Metaplectic Iwahori Whittaker functions 23 3.1. Metaplectic preliminaries 23 3.2. Intertwining operators and recurrence relations 29 3.3. Metaplectic Demazure-Whittaker operators 32 3.4. Averaged metaplectic Demazure-Whittaker operators 35 4. Main theorem and conclusions 37 arXiv:2012.15778v2 [math.RT] 18 Jan 2021 4.1. Metaplectic Iwahori Whittaker functions for G = GLr 38 4.2. The main theorem 39 4.3. R-matrices and p-adic intertwiners 40 References 43

2020 Mathematics Subject Classiﬁcation. Primary: 22E50; Secondary: 82B23, 16T25, 05E05, 17B37, 11F70. 1 1. Introduction This work is the culmination of a series of papers [15, 9, 13, 10, 12] concerning a remarkable parallel between the theory of non-archimedean Whittaker functions and certain solvable lattice models. Our results are naturally divided into three parts represented by the subsequent sections of the paper:

· The construction of a novel supersymmetric solvable lattice model connected to the quantum ˆ group Uq gl(n|r) using a fusion construction similar to [36]. · The computation of all values of metaplectic Iwahori Whittaker functions for split reductive groups using vector Demazure–Whittaker operators. · Connecting the above results in the special case of a metaplectic n-cover of GLr through a natural bijection between the data that specify the values for the metaplectic Iwahori Whittaker function and the boundary data for the lattice model.

Taken together, these results show that quantum superalgebras and their associated lattice models are able to capture deep aspects of metaplectic group representations at the Iwahori level. Solvable lattice models arose in statistical physics. A special role in statistical mechanics is played by the partition function, a sum over states that controls the thermodynamic properties of the system. In solvable models, a remarkable tool, the Yang-Baxter equation, reveals unexpected symmetries of the partition function and connects to quantum groups. Solvable lattice models appear in a wide variety of different contexts, such as conformal and quantum field theory (e.g. [3, 6, 30, 29]), algebraic combinatorics ([51, 14, 20]), and integrable probability ([7, 28, 8, 1]). Surprisingly, there exist solvable lattice models whose partition functions match functions from p-adic representation theory. As the starting point of this discussion let us take the “Tokuyama model,” a kind of six-vertex model where a state of the model is an assignment of “spins” to the edges of a two-dimensional grid, with six allowed configurations of spins at each vertex. The partition functions are recognizable as spherical Whittaker functions on GLr(F ), with F a non-archimedean local field. These models appeared in [49, 32, 16] and their relation to GLr(F ) is explained in [17], Chapter 19. The spherical Whittaker functions represented by the Tokuyama model are unique in two ways: in this setting Whittaker functionals are unique, and the spherical vector in the representation is also unique. We will report on two generalizations, each removing a uniqueness. In [9], solvable lattice models called metaplectic ice were constructed whose partition functions are spherical Whittaker functions on a degree n metaplectic central extension of GLr(F ), provided F contains enough roots of unity. Due to non-commutativity of the metaplectic torus, there is no longer a unique Whittaker functional, leading to an nr- dimensional space of metaplectic spherical Whittaker functions. The setup in [9] exhibits a model for each Whittaker function in a chosen basis. A common feature for the Tokuyama model and metaplectic ice is the focus on spherical Whittaker functions, meaning that they arise from a spherical vector in the unramified principal series invariant under the maximal compact subgroup. Yet it is important to have models for other Whittaker functions, particularly those fixed by the smaller Iwahori subgroup, since these better capture the representation theory of the p-adic group. 2 For the nonmetaplectic Iwahori case treated in [12], inspired by ideas in [8], we exhibited lattice models whose partitions functions are Whittaker functions invariant under an Iwahori subgroup, thereby relaxing another condition for uniqueness. For GLr(F ) there is an r!- dimensional space of Iwahori fixed vectors. In this paper, we will simultaneously generalize both [9] and [12] by exhibiting solvable lattice models whose partition functions are Iwahori Whittaker functions on a metaplectic n-cover of GLr(F ). One way to tidily express these increasing levels of generality is via the quantum group R-matrices that imply solvability of the respective models. To be precise, the R-matrix will be a Drinfeld twist of a quantum group module as highlighted in Figure1. This will ultimately give us a connection between quantum and non-archimedean representation theories in Section 4.3 via the maps induced by intertwiners on the spaces of Whittaker functionals and Iwahori fixed vectors.

GLr metaplectic n-cover of GLr no scolors n scolors (or charges)

spherical vector Tokuyama ice [16] Metaplectic ice [15,9] n=1 no colors Uq(glb(1|1)) Uq(glb(1|n))

identify colors identify colors

Iwahori ﬁxed vector Iwahori ice [12] Metaplectic Iwahori ice n=1 r colors Uq(glb(r|1)) Uq(glb(r|n))

Figure 1. Summary of the different lattice models for non-archimedean Whittaker functions. The arrows indicate how results from this paper generalize and recover previous results of the indicated citations. The horizontal arrows reflect the fact that the models in the left columns are special cases (with n = 1) of the models in the right column. If the boundary spins of the colored system are restricted to a single color, we recover the Tokuyama or metaplectic ice models. To look at it another way, the vertical arrows reflect the fact that the spherical Whittaker function is the sum of Iwahori Whittaker functions.

The models in each of these four cases consist of a grid, whose edges are allowed to take on a limited number of values called spins. In the Tokuyama model, there are only two spins. In a state of the system, the spins can be organized into paths or lines through the grid, beginning on the top boundary and ending on the right. In the models of [12], the number of spins is increased, and the new spins are called colors. A state of the system may be described by the configuration of distinctly colored paths. Using such models with r colors we may represent a basis of Iwahori Whittaker functions on GLr. The models of [9] do not use colors, but a different scheme involving an attribute called charge, which is replaced by an equivalent scheme in this paper. If we reimagine [9] using the technology of this paper, there are again paths, but these move down and to the left to terminate on the left boundary. To distinguish these paths from the colored ones in [12], we will call their attribute super color, or simply scolor, instead of color. To model all nr Whittaker functionals, a given scolor may appear on multiple lines. 3 The colors and scolors are combined together in this paper, and roughly the colors are used to distinguish Iwahori fixed vectors while the scolors are used to distinguish different Whittaker functionals. The reason for the name “scolor” is to showcase the symmetry with colors, which is also manifested in the quantum super group Uq(glb(n|r)) associated to the Yang-Baxter equation for this model. With the quantum group interpretation, the r colors and the n scolors correspond to a basis of the n + r dimensional defining module for this quantum group. Thus we think of the scolors or super colors as supersymmetric counterparts of the colors. In order to achieve the results of this paper there is work to be done both in the theory of these lattice models and in the theory of metaplectic Iwahori Whittaker functions for the unramified principal series. Both sides of this parallel involve recurrences made with certain divided difference operators, called Demazure-Whittaker operators, in which there appear Gauss sums. In the representation theory side, such recurrences were developed in [18, 26, 44, 41]. This paper redevelops this theory by introducing new vector metaplectic Demazure-Whittaker operators. Previous literature [26, 44, 41, 42] has focused on a particular (averaged) metaplectic Whittaker functional and associated Demazure-Whittaker operators, while our operators treat the vector of all Whittaker functionals. We develop the metaplectic theory in the context of a general split reductive group G, with some new results. For example, while previous works have mainly focused on the values of the metaplectic Iwahori Whittaker functions restricted to the metaplectic torus, we compute all values for the metaplectic Iwahori Whittaker functions and specify a convenient set of combinatorial data to describe the associated group element arguments in G(F ). From our vector Demazure-Whittaker operators we also obtain a natural derivation of the intricate Chinta-Gunnells [25] Weyl group action and the averaged Demazure-Whittaker operators of [26, 44, 41] from the non-archimedean representation theory. Furthermore, with these operators we produce a representation of the metaplectic Iwahori Hecke algebra on the space of polynomial functions on the dual torus. This representation is isomorphic to the one appearing in [46]; the difference is that our representation arises naturally from the non-archimedean representation theory, while theirs follows solely from considering the affine Hecke algebra and the double affine Hecke algebra. The final part of the paper will connect the representation theory to the lattice models. Lattice models differ depending on the Cartan type, so in order to connect them to the representation theory we must specialize to the case of G = GLr. The strategy to show that the partition functions equal the Whittaker functions is to prove that they satisfy the same recurrence relations. On the representation theory side these recurrences come from the vector Demazure-Whittaker operators, and on the lattice model side, from Yang-Baxter equations. The partition function is determined by the boundary data for the lattice model. It is striking that there is a bijection between all possible boundary data and the data for the different values of the metaplectic Iwahori Whittaker function. This is schematically shown in Figure2. As we have explained, the r colors are associated with the data required to specify an Iwahori fixed vector and the n scolors are associated with the data required to specify a Whittaker functional. Hence the colors and scolors are of different origin, but remarkably the spin types required to accommodate both in a single theory — the colors and scolors — are

4 g ∈ G(F )  (w0, λ)     θ w   z1  zρφ (z; g) = Z   θ,w  z2     z3    + + + + +

Iwahori vector Whittaker functional

Figure 2. The metaplectic Iwahori Whittaker function φ is equal to the partition function Z up to the prefactor zρ where ρ = (r −1, r −2,..., 0) by TheoremA using the following dictionary. The r Langlands-Satake parameters z = (z1, . . . , zr) ∈ C determine the unramified principal series for φ and the weights for vertices at row i in the lattice model. The θ ∈ (Z/nZ)r labels the Whittaker functional for φ and the scolors on the left boundary of Z. The Weyl group element w labels the Iwahori fixed basis vector in the unramified principal series for φ and the order of the colors on the right boundary of Z. Lastly, φ is determined by its values on particular group elements g which are specified by certain pairs of weights λ and Weyl group elements w0. On the lattice model side λ and w0 determine the top boundary conditions.

not very different from each other. One apparent difference is that there are nr Whittaker functionals but only r! Iwahori fixed vectors; this is reflected in the boundary conditions, where scolors on the left boundary are allowed to repeat, whereas the colors on the right boundary are not allowed to repeat. Still the Yang-Baxter equations require the theory to accommodate multiple colors. The whole machinery therefore works for more general models in which both colors and scolors are allowed to repeat. One may ask, what is the meaning of such systems? One answer to this was given in [12], where such systems with repeating colors were used to represent parahoric Whittaker functions. It seems possible that further pursuing the relationship between colors and scolors, or in another optic to treat metaplectic Whittaker functionals on the same footing as parahoric vectors, may prove to be a useful exercise that is accessible by our methods. The R-matrices that appear in the Yang-Baxter equations are (up to Drinfeld twisting) known as “Perk-Schultz” or supersymmetric R-matrices with an already-substantial literature [43, 5, 53, 52, 35]. However the RTT Yang-Baxter equations presented in Section 2.4 are new and different from other known Yang-Baxter equations that involve the R-matrix alone. These models involve a new type of vertex, distinct from the R-matrix, whose Boltzmann weights appear in Figure4. Therefore the models in this paper are, as far as we are aware, completely new. A recent preprint of Aggarwal, Borodin and Wheeler [2] independently introduced

supersymmetric models that make use of modules for the quantum group Uq(slb(m|n). One may make a distinction (also explained in [12]) between models that are bosonic and fermionic. In fermionic systems, an edge may carry at most one instance of any given color or scolor, whereas in bosonic systems, an edge may carry a (s)color with a multiplicity. The models in [12] resembled those in [8], but were distinguished in that those of [12] are fermionic and those of [8] are bosonic. In contrast, the models in this paper and in the new paper [2] are 5 both fermionic. Their approach is quite diﬀerent from ours, and it will be interesting to explore connections between the two. Although our main purpose here is to expose a deep connection between Whittaker functions and solvable lattice models, it is possible that these methods will produce insights and applications to number theory and related representation theory. Realizations as partition functions are suitable for proving branching rules (related to Jacquet functors to Levi subgroups) and Cauchy-type identities related to the Rankin-Selberg method. Thus in [11] the key lemma (Statement B) of [17] was proved using the Yang-Baxter equation; this had applications to multiple Dirichlet series and was simpler than the original proof. See [21] for a demonstration of how Cauchy identities may be proved. Branching rules and Cauchy identities are also an important theme in the parallel theory in integrable probability [2, 8]. It seems an important question whether deeper level Whittaker functions and Hecke algebras can similarly be realized in terms of solvable lattice models.

The paper is outlined as follows. In Section2 we introduce our new class of supersymmetric lattice models which are constructed through a fusion process similar to that of [36]. The model that is used as an input in this fusion process is called monochrome and is detailed in Section 2.1, while the fusion process itself is explained in Section 2.2. In preparation for the proof of TheoremA we evaluate the partition function of a certain ground state system in Section 2.3. In Section 2.4 we prove that our lattice model is solvable; that is, we prove an associated Yang-Baxter equation in Corollary 2.9. To show this for every n and r, we first prove an auxiliary Yang-Baxter equation for the monochrome model in Theorem 2.6 which we reduce to the problem of verifying a finite number of equations. In the same section we also show that the R-matrix is a Drinfeld twist of the Uq(glb(r|n)) R-matrix in Proposition 2.11. From the Yang-Baxter equation we obtain a recursion relation for the partition function in Proposition 2.12 of Section 2.5 that we later compare with the recursion relations for the metaplectic Iwahori Whittaker functions. In Section3 we turn our attention to metaplectic Iwahori Whittaker functions for general reductive groups, with some introductory theory and notation introduced in Section 3.1. Here, we also compute a base case for the Whittaker functions in Proposition 3.8. The recursion relations for the Whittaker functions are obtained in Proposition 3.11 in Section 3.3 by studying intertwining operators. These Whittaker functions can be conveniently described using the (vector) metaplectic Demazure-Whittaker operators as shown in Theorem 3.12 in Section 3.3. These operators are shown to produce a representation of the affine metaplectic Hecke algebra in Theorem 3.15. An expression for all values of the metaplectic Iwahori Whittaker functions in terms of these operators is deduced in Corollary 3.13. In Section 3.4 we match our vector metaplectic Demazure-Whittaker operators with the averaged metaplectic Demazure-Whittaker operators of [26, 41] in Theorem 3.16, which also recovers the Chinta- Gunnells Weyl group action from non-archimedean representation theory. We conclude the paper with Section4 where we state and prove the main theorem, TheoremA in Section 4.2, after detailing the computation of the metaplectic Iwahori Whittaker functions for GLr and a particular convenient cocycle in Section 4.1. By comparing the partition function for the ground state system with the above base case for the Whittaker function and then using the matching recurrence relations, we can conclude the equality for all cases. Then in Section 4.3, we give an interpretation of the intertwining operators in the 6 non-archimedean representation theory in terms of different R-matrices on the lattice model side as illustrated by equation (4.9) and Figure 20. Acknowledgements. The authors thank Slava Naprienko, Manish Patnaik, Anna Puskás, Siddhartha Sahi, Jason Saied and Chenyang Zhong for useful discussions and comments. The authors would like to thank Amol Aggarwal, Alexei Borodin and Michael Wheeler for illuminating discussion regarding their work [2]. This work was supported by NSF grants DMS-1801527 (Brubaker) and DMS-1601026 (Bump). Buciumas was supported by the Australian Research Council DP180103150 and DP17010264, NSERC Discovery RGPIN-2019-06112 and the endowment of the M.V. Subbarao Professorship in Number Theory. Gustafsson was supported by the Swedish Research Council (Vetenskapsr˚adet)grant 2018-06774.

2. Supersymmetric lattice models 2.1. The monochrome model. We will describe a class of solvable lattice models, similar to the exactly solved models in [4]. Such a model consists of systems S, certain finite collections of states and associated Boltzmann weights for each state in the system, to be described more precisely in a moment. The sum Z(S) of these Boltzmann weights over all states in the system S is the partition function of the system. The solvability of the model amounts to the existence of Yang-Baxter equations that allow us to study these partition functions through the commutation properties of the row transfer matrices, following the method by which Baxter [4] studied such models. This paper treats three equivalent versions of the model, which we will call monochrome, color-fused and fully-fused. Collectively, we call these models metaplectic Iwahori ice in reference to six-vertex models often being called ice models. In this section we describe the monochrome model. Its admissible states are described by labelings of a rectangular grid with Nnr columns and r rows, where n and r are fixed integers (later, the cover degree and the rank of the associated metaplectic general linear group, respectively) and N is a sufficiently large integer. An example of such a grid with n = r = 3 and N = 2 is given in Figure3. In particular, there is a vertex at each crossing of a row and column in the grid, and each vertex has 4 adjacent edges. Edges that are adjacent to just one vertex are called boundary edges while edges adjacent to two vertices are called internal edges. A fixed labeling of all boundary edges is referred to as a boundary condition and, given a boundary condition, an (admissible) state is a labeling of internal edges satisfying certain given restrictions. We will dictate the restrictions for our monochrome models momentarily. Finally, a system S is defined as the set of admissible states for a fixed boundary condition, together with a set of Boltzmann weights β(s) for every admissible state s ∈ S. Before we describe monochrome systems formally, we give an intuitive description of its admissible states. An (admissible) state of the system assigns labels (or ‘attributes’) called color and scolor to the edges. Borrowing the naming convention from physics, the prefix s is used to imply that scolor is the supersymmetry partner of color. A color c is an attribute indexed by an integer modulo r. The set of possible colors will be denoted {c1, ··· , cr} and ordered so that ci < cj if and only if i < j.A scolor w is an attribute indexed by an integer modulo n. The possible scolors are denoted {w0, ··· , wn−1} and are ordered so that wi < wj if and only if i < j. 7 In any state, the edges that are labeled with a color c or scolor w will also be said to carry the color or scolor. According to the restrictions we will impose, the edges carrying any color c will form a connected path through our rectangular grid, traveling from the top boundary and exiting at the right boundary. Similarly the edges that carry the scolor w form a path beginning at the top boundary of the grid and end up at the left boundary. To separate the two, we will draw color-paths with thick solid lines and scolor-paths with dotted lines; see Figure7 for an example. More than one scolor-path may have the same scolor w, but for the current application we do not need systems in which different color-paths have the same color. Edges may be classified as horizontal (occurring in a row of the grid) or vertical (occuring in a column). In any (admissible) state s, every horizontal edge carries a unique attribute, either a color or scolor. Every vertical edge carries either no attribute, or two: a color and a scolor. In Figure7, the vertical edges that carry the color ci and the scolor wj are labeled ciwj, or simply ij when no confusion may arise. So roughly, admissible states are paths of a given color or scolor traveling through the grid as described above and the boundary conditions serve to dictate the beginning and ending location of each path. We now give a more formal description of the systems. First, we label each column of the model with a pair (c, w) consisting of a color c and scolor w according to the following cyclic ordering. Let (c, w) = (ci, wj) be such a pair with 1 6 i 6 r and 0 6 j 6 n − 1. Then the successor (c0, w0) to (c, w) in this cyclic ordering is defined by  (ci+1, wj) if i < r, 0 0  (2.1) (c , w ) = (c1, wj−1) if i = r, j > 0  (c1, wn−1) if i = r, j = 0.

Note that the colors ci first increase to r, then the accompanying scolor wj decreases by one and color is reset to c1, and then cycle repeats. The leftmost column is labeled (c1, wn−1) and each column to the right of any (c, w) is labeled by its successor (c0, w0) in the cyclic ordering. According to the size of our grid, each pair appears N times in the labeling. Next we label each row with a formal parameter zi where the index i increases from 1 to r from top row to bottom row; when we specialize these r parameters (z1, . . . , zr) to non-zero complex numbers, we call them spectral parameters. See Figure3 for an example of the labels on columns and rows. In the figure we have grouped each sequence of r columns of the same scolor into blocks enumerated from right to left. These blocks will become the columns of the color-fused system described in Section 2.2. Finally, we label each vertex T in the grid cw with the labels from its row and column, writing Tz for the vertex with column label (c, w) and row label z. In decorating the edges of our model with attributes color and scolor, we only allow for the configurations in Figure4 on the edges adjacent to any vertex. In particular, in the column labeled (c, w), the vertical edges are only allowed to carry both attributes c and w or neither. This explains the name monochrome for the system, since each column can only contain a specific color and scolor. A horizontal edge must carry a single color ci or a single scolor wj (independent of column labels) but not both. It is sometimes convenient to refer to the collection of all possible sets of attributes on a given edge. We call this collection the spinset for an edge. So the spinset for every horizontal edge is

(2.2) {c1, ··· , cr, w0, ··· wn−1} 8 block 5 block 4 block 3 block 2 block 1 block 0

c = c1 c2 c3 c1 c2 c3 c1 c2 c3 c1 c2 c3 c1 c2 c3 c1 c2 c3 w = w2 w2 w2 w1 w1 w1 w0 w0 w0 w2 w2 w2 w1 w1 w1 w0 w0 w0

z = z1

Figure 3. The monochrome grid. Every vertex (represented here by a dot) in each row is assigned the same spectral parameter zi, and every vertex in each column is assigned the same scolor and color; the colors c and scolors w rotate cyclicly through all the possibilities. The columns with the same scolors are combined into blocks numbered from right to left, and these blocks will later become the columns of the color-fused model. In this example, n = r = 3 and N = 2.

a1 a2 b1 b2 c1 c2 + cw cw + + cw

cw cw cw cw cw cw y Tz y a Tz a y Tz y a Tz a c Tz w w Tz c

+ cw cw + cw + v if c > a z if a = c, 1 z if c = a g(y − w) (1 − v)z 1 1 otherwise. 1 if c < a

Figure 4. The Boltzmann weight at a vertex of color c and scolor w labeled by cw. Here a is a color (which may be the same or different from c) and y is a scolor (which may be the same or different from w). We label each vertex configuration according to the first row of the table. If the vertex configuration of edges does not appear in this table, the Boltzmann weight is zero. For the Boltzmann weight g(y − w) we identify each scolor wj with the integer j.

cw while the spinset for a vertical edge adjacent to vertex Tz is the two-element set {+, (c, w)}, where + is a placeholder for no attribute. We will now describe the boundary conditions which assign to each boundary edge an element of its spinset. These conditions are parametrized by µ ∈ Zr with non-negative parts, r an element θ in (Z/nZ) and a Weyl group element w in Sr as summarized in Figure5. At row i, the left boundary edge is assigned the scolor θi and the right boundary edge the color cr+1−w−1(i), that is, the color (wP )i where P = (cr, cr−1, . . . , c1) is the palette of colors in decreasing order. Recall that vertical edges have either no attribute, which we denote by + and will call empty, or the color and scolor of the corresponding column, which we will called filled in. Thus we only need to specify which columns are filled in and which are empty. All edges at the bottom boundary are empty, that is, assigned +. For each j we fill in the column at the top boundary which is in block µj with color cr+1−j, that is, color (P )j. We recall the following definition from [12] for a root datum with weight lattice Λ, root system ∆ with positive and negative roots ∆± and Weyl group W as usual.

9 block number µj

color (P )j

z1 row i z2 row i . scolor θi . color (wP )i

r Figure 5. Boundary conditions as determined by µ ∈ Z with non-negative parts, w ∈ Sr and θ ∈ (Z/nZ)r. The block number refers to the blocks shown in Figure3 for the monochrome system, or equivalently to the column numbers of the color-fused system discussed in Section 2.2. Here P = (cr, cr−1, . . . , c1) is a tuple of the colors ordered from largest to smallest.

0 ∨ Deﬁnition 2.1. Let w ∈ W , let αi be a simple root and αi the corresponding coroot. A weight λ ∈ Λ is w0-almost dominant if

( 0 −1 + ∨ 0 if (w ) αi ∈ ∆ , (2.3) hαi , λi > 0 −1 − for all simple roots αi. −1 if (w ) αi ∈ ∆ . r In this section we use the GLr weight lattice Λ which we identify with Z . Remark 2.2. We showed in Proposition 7.2 of [12] that for each weight µ there exists a unique pair (w0, λ) with λ ∈ Λ being w0-almost dominant and w0µ = λ+ρ where ρ = (r−1, r−2,..., 0). Concretely, this means that w0µ is a dominant weight, but if µ contains repeated entries, there are several permutations w0 that would make w0µ dominant. The weight λ = w0µ − ρ may not be dominant, but the Weyl group element w0 that makes λ w0-almost dominant is unique in the sense that the order of the colors at the top boundary read from left to 0 0 0 right is w P . Indeed, for each j we fill in the column in block (w µ)j with color (w P )j, and 0 0 0 0 since (w µ)j > (w µ)j+1 we only need to check the case when (w µ)j = (w µ)j+1 which means 0−1 0−1 0−1 that hαi, λi = −1 and w αi < 0. In that case w (j) > w (j + 1) which is equivalent to 0 0 (w P )j = cr+1−(w0)−1(j) < cr+1−(w0)−1(j+1) = (w P )j+1. With the boundary spins fixed, a state s of the system assigns a spin — that is, an element of the spinset — to every edge of the grid. See Figure7 for an example. We call the set of states with given boundary spins a system, denoted by Sµ,θ,w. Given a state s ∈ Sµ,θ,w, cw each vertex Tz is assigned a Boltzmann weight, according to the values in Figure4 where we identify each scolor wj with the integer j. These depend on c, w and z for the vertex, as well as a parameter v and a function g of the integers modulo n that satisfies the following condition. Assumption 2.3. The function g satisfies g(0) = −v, and if a is not congruent to 0 modulo n, then g(a)g(−a) = v. Note that this assumption is satisfied when g is a certain non-archimedean Gauss sum and the parameter v is associated to the cardinality of the residue field. This will be important when connecting these lattice models to metaplectic Whittaker functions in a later section. The Boltzmann weight β(s) of the state s is the product of the Boltzmann weights at all vertices of s. The partition function Z(S) of a system S is the sum of the Boltzmann 10 weight β(s) over all states s in the system S and will be considered as a function of the row parameters z = (z1, . . . , zr). 2.2. Equivalent model via fusion. In [12] we employed a procedure we termed fusion to introduce additional equivalent versions of a lattice model. The terminology was chosen owing to the resemblance with lattice model incarnations of the fusion procedure for tensor products of quantum group modules and their subquotients as described by Kulish, Reshetikhin, and Sklyanin [36], though the precise algebraic explanation of the fusion process in [12] remains a mystery. Nevertheless the foundations presented in [12] allow us to conclude Yang-Baxter equations for “fused” models from similar equations for “unfused” models. Here, we will employ a similar procedure to the monochrome systems of the previous section.

Block of unfused (monochrome) vertices Equivalent fused vertex B B B 1 2 ··· r B

E E Er−1 A 1 2 ··· C A C c1w c2w crw w Tz Tz Tz Tz ··· D1 D2 Dr D

Figure 6. Fusion. This procedure replaces a sequence of vertices by a single vertex. Left: vertices in the unfused system. Right: a single replacement vertex in the fused system, with replacement edges B and D. The spinset for the fused vertical edge B is the Cartesian product of the spinsets of the edges Bi, and similarly for D.

The fusion procedure replaces a sequence of adjacent vertices in a model by a single vertex. We may preserve the weight of the state in the process by taking the Boltzmann weight of the fused vertex to be the partition function of the adjacent vertices that were fused together. c1w crw In our case, we fuse a block of adjacent vertices Tz ,...,Tz in a given row of the grid w over a complete set of r colors {ci} into a single fused vertex which is labeled Tz . The fusion process is illustrated in Figure6. In the notation of the figure, the spinset for the fused vertical edge B is the Cartesian product of the spinsets of the edges Bi, and similarly for D. Thus specifying spins for the edges A, B, C and D determines the spins for Bi and Di for the (unfused) monochrome configuration. In particular, the spin for Bi is either of the form ciw or + so the spinset for B is either cw, where c is some subset of colors, or +. If there are spins for the edges Ei that make the Boltzmann weights nonzero using the monochrome vertex configurations of Figure4, it is easy to see that they are uniquely determined. So the partition function for these adjacent vertices is just the product of the Boltzmann weights at w the unfused vertices. Hence the Boltzmann weight of the fused vertex Tz is just defined to c1w crw be the product of these monochrome weights at Tz ,...,Tz . In Figure7 we have a state of a monochrome system, and in Figure8 we have the corresponding state of the corresponding color-fused system. Since the (unfused) monochrome grid has Nnr columns, the corresponding color-fused grid has Nn columns. Again the rows are numbered from 1 to r in increasing order from top to bottom, and the columns are labeled in descending order from Nn − 1 to 0 moving left to right. Another example of fusion is shown in Figure9. 11 block 5 block 4 block 3 block 2 block 1 block 0

c1 c2 c3 c1 c2 c3 c1 c2 c3 c1 c2 c3 c1 c2 c3 c1 c2 c3 w2 w2 w2 w1 w1 w1 w0 w0 w0 w2 w2 w2 w1 w1 w1 w0 w0 w0

+ + + + + 3+1 + + + + 2+2 + + + + 1+0 + +

z1 1 1 1 1 1 1 3 3 3 0 0 0 0 0 0 0 1 1 1

+ + + + + + + + 3+0 + 2+2 + + + + + + +

z2 0 0 0 0 0 0 0 0 0 3 3 3 3 3 3 3 3 3 3

+ + + + + + + + + + 2+2 + + + + + + +

z3 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2

+ + + + + + + + + + + + + + + + + +

Figure 7. A state of a system Sµ,w,θ using the expanded grid and monochrome weights with µ = (4, 2, 0), θ = (1, 0, 2), w = s1s2, and n = r = 3. For readability we have labeled a vertical edge with the attribute ckwl by kl where the scolor label is underlined by a dotted line in recognition of how the scolor-paths are styled. The vertex in the column labeled ciwj and row labeled zk is ciwj of type Tzk , with Boltzmann weights as in Figure4. Note how the color and scolor paths are mirrored.

5 4 3 2 1 0

w2 w1 w0 w2 w1 w0

+ 3+1 + 2+2 + 1+0

z1 1 1 3 0 0 0 1

+ + 3+0 2+2 + +

z2 0 0 0 3 3 3 3

+ + + 2+2 + +

z3 2 2 2 2 2 2 2

+ + + + + +

Figure 8. The color-fused system equivalent to Figure7. For reference we have added row and column numbers. Note that column numbers decrease from left to right, and row numbers increase from top to bottom. This system is Sµ,θ,w where µ = (4, 2, 0), θ = (1, 0, 2), and w = s1s2.

In Figure 10 we write the Boltzmann weights of the color-fused model where a vertical w edge is described by a set Σ of colors. At a color-fused vertex Tz the set determines which columns (c, w) with c ∈ Σ that are occupied in the corresponding monochrome block. If the set Σ is empty, we denote the vertical edge by a plus sign. Note that these weights precisely match the weights in Figure 12 of [12], when one allows the existence of only one scolor w. 12 block 4 column 4

c1 c2 c3 w1 w1 w1 w1

+ + 3+1 31

z1 2 2 1 1 z1 2 1

+ 2+1 3+1 = (2 ⊗ 1) ∧ (3 ⊗ 1)

z2 1 1 1 3 z2 1 3

+ 2+1 + 21

Figure 9. An example of color-fused block where a resulting vertical state contains more than one color. We denote by (2 ⊗ 1) ∧ (3 ⊗ 1) the fused edge carrying the colors c2 and c3 together rn with the scolor w1. There are 2 possible states for the vertical fused edges, and these are in bijection with the exterior algebra over the free vector space space spanned by pairs consisting of one color and one supercolor, hence the notation. Note that the number of color paths and the number of scolor paths going through a vertical edge have to be equal.

The appearance of the term g(w − y)|Σ| in the bottom right entry of the table is surprising, but the exponent on the Gauss sum indeed checks out in the matching with Figure 12 of [12] because g(w − y) = g(0) = −v in that degenerate setting.

−+ Σ Σi j Σ

w w w ci Tz ci ci Tz cj cj Tz ci

Σ Σ Σ zv|Σ[i+1,r]| (1 − v)z(−v)|Σ[i+1,j−1]|v|Σ[j+1,r]| 0 + − Σi Σi Σ

w w w w Tz ci ci Tz w y Tz y

Σ Σ Σ (−v)|Σ[1,i−1]|v|Σ[i+1,r]| (1 − v)z(−v)|Σ[i+1,r]| g(y − w)|Σ[1,r]|

Figure 10. The Boltzmann weights for the color-fused model with scolor w, with vertex labeled w Tz . Here Σ is an arbitrary set of colors (which may be empty), and |Σ| denotes the cardinality of Σ. Moreover, Σ[i,j] denotes the colors in Σ with indices in the interval [i, j]. We introduce + − the operation Σi := Σ ∪ {ci} to be used only if ci 6∈ Σ and similarly for Σi ∪ {ci} = Σ and −+ + Σi j ∪ {ci} = Σj . We denote by ci, cj colors with i < j and y is a scolor (which may be the same or diﬀerent from w).

We will now translate the boundary conditions from the monochrome system to the color-fused system where a block of columns in the former becomes a single column in the 13 r latter. Since they are equivalent, the system will still be denoted Sµ,θ,w where µ ∈ Z with r non-negative parts, θ is an element in (Z/nZ) and w a Weyl group element in Sr. Here r is the number of rows in the grid with Nn color-fused columns. As before, along the bottom boundary of the color-fused model, all column vertices are assigned +, meaning the empty attribute. Along the right-hand border, the horizontal boundary edge in row i receives color cr+1−w−1(i), and along the left-hand border, the horizontal boundary edge in row i receives scolor θi. This leaves the top boundary. As before, the vertical boundary edges receive attributes according to the parts of µ so that each color appears exactly once along the top boundary edges. However, more than one color may appear in a single edge by gathering the colors in a block of columns in the monochrome system into a single column in the color-fused system. Speciﬁcally, column i gets attributes cwj where i ≡ j (mod n) and c is the subset of colors

{cr+1−j1 , . . . , cr+1−j` } such that the parts µj1 = ··· = µj` = i, and if this subset is empty we denote the edge by +. Recall that according to Remark 2.2, we may either parametrize the top boundary by r 0 0 0 µ ∈ Z or by a pair of elements (w , λ) where w is in Sr and λ is a w -almost-dominant weight (Deﬁnition 2.1) such that w0(µ) = λ + ρ with ρ = (r − 1,..., 1, 0). We have seen that color-fused systems are equivalent to monochrome systems by using a fusion process to collapse r consecutive columns over a full set of colors. It is natural to make a further collapse of the color-fused systems by means of a fusion process over a full set of n scolors, resulting in a grid with N columns and r rows. We refer to this third, equivalent collection of systems as the fully-fused model. See Figure 11 for the result of performing this collapse on our running example. Again, none of the partition functions is altered by the fusion process, as the Boltzmann weights are taken to be the product of Boltzmann weights in the un-fused model. The fully-fused model is best for making connections to associated quantum group modules.

Remark 2.4. Once a state of the system is fixed, colored and scolored lines emerge. The colored lines begin at the top edge and end at the right edge; the scolored lines begin at the top edge and end at the left edge. In the monochrome and color-fused systems, if the colored paths are specified, then the scolored paths are determined, so (if an allowable configuration results) the state is uniquely determined. This follows from the facts (apparent from the Boltzmann weights) that every vertical edge carries a color and a scolor, and every horizontal edge carries a color or a scolor (but not both); that the number of edges at a vertex carrying a given scolor (or color) is even; and that (in the monochrome and color-fused systems) the scolor carried by any vertical edge is predetermined. From these properties of the Boltzmann weights, it is obvious that each scolored line beginning at the top edge (and specified by boundary conditions) can be traced uniquely if the colored lines are known. It is similarly true in the monochrome systems that the scolored lines determine the state.

2.3. The partition function for the ground state system.

Lemma 2.5. Let θ ∈ ( /n )r and let µ ∈ r . Then there is a unique w0 ∈ W such that Z Z Z>0 the colors for the top boundary edges for the system Sµ,θ,w0 in the monochrome picture is 0 0 w P from left to right where P = (cr, cr−1, . . . , c1) which is equivalent to λ := w µ − ρ being 14 31 22 ∧ 10

1 0 1

30 22

0 3 3

+ 22

2 2 2

+ +

Figure 11. The state of the fully-fused system that corresponds to the states in Figures7 and8.

w0-almost dominant (Deﬁnition 2.1). Furthermore,

( `(w0) λ+ρ v z if θ ≡ λ + ρ mod n (2.4) Z S 0 (z) = µ,θ,w 0 otherwise.

Proof. We will use the monochrome version of the system and prove that the system Sµ,θ,w0 contains only at most a single state, which we call a ground state. Since the colors in P are unique and the colors of the top boundary are determined by µ there is a unique w0 such that the order of colors is w0P . That this is equivalent to λ being w0-almost dominant follows from Remark 2.2 which is based on Proposition 7.2 of [12]. Note that λ + ρ = w0µ is a dominant weight which means that we fill in top boundary columns at 0 the blocks (λ + ρ)i with color (w P )i in the order from left to right increasing with i. From Remark 2.4 we have that a state is determined by its colored paths (or, alternatively its scolored paths). By the assumptions, the colors on the top boundary from left to right is the same as the colors on the right boundary from top to bottom. We will now follow the possible colored paths from the top boundary to the right boundary and argue that there is at most one admissible state for Sµ,θ,w0 . From the admissible vertex configurations in Figure4 we see that the colored paths can only go rightwards and downwards. This means that the path starting at the left-most color on the top boundary must immediately go straight to the right boundary. Now the path starting at the second left-most color needs to go down and right to the right boundary on the second row. However, all the horizontal edges on the first row to the right of the second left-most color on the top boundary are occupied by the previous colored path. Thus it must first go straight down to the second row and then straight right to the right boundary. Continuing in this fashion we see that all colored paths are L-shaped in this way, and that there is at most one admissible state since an admissible state is determined by its colored paths. We will now argue that we get an admissible state if and only if θi ≡ (λ + ρ)i mod n at all rows i. An example of such a state, which we call a ground state, is shown in Figure 12. As seen in the vertex configurations of Figure5, the scolor paths, starting from the top boundary and ending at the left boundary, can only go leftwards and downwards, and the latter only when together with a color path. If we consider the left-most scolor path starting from the top boundary we see the it has to follow the first color path to the first row and immediately go straight to the left boundary. The second scolor path needs to follow the second color path 15 until the second row and then go straight to the left boundary. Continuing in this fashion we see that the scolored paths all have a vertically mirrored L-shape and that we clearly get an admissible state of the system provided it satisfies the boundary conditions. The left boundary scolor at row i matches the scolor of the i-th left-most scolor on the top boundary which is (λ + ρ)i mod n by construction and thus the boundary conditions are satisfied if and only if θ ≡ λ + ρ mod n. It remains to compute the weight of this admissible state. The involved vertices are of types a1, c2, a2 and b2 shown in Figure4 where only the latter two make any non-trivial contributions. For both type a2 and b2 we obtain, at row i, a factor of zi whenever the vertex color equals the horizontal edge color. This happens once for every color block of r columns (λ+ρ)i in the fully expanded system, which means that we get a contribution of zi . The only remaining non-trivial contributions come from the c > a case in type a2 each with a factor of v. These occur whenever we have a crossing of two colored paths where the color a of the horizontal edge is less than the color c of the vertical edge. The ordering of the colors on the top and right boundaries is given by w0P , and thus the number of such crossings is `(w0). Hence, the weight of the unique admissible state, called the ground state, is v`(w0)zλ+ρ which concludes the proof.

block 3 block 2 block 1 block 0

c1 c2 c3 c1 c2 c3 c1 c2 c3 c1 c2 c3 w1 w1 w1 w0 w0 w0 w1 w1 w1 w0 w0 w0

+ 2+1 + + + 3+0 + + + 1+0 + +

z1 1 1 2 2 2 2 2 2 2 2 2 2 2

+ + + + + 3+0 + + + 1+0 + +

z2 0 0 0 0 0 0 3 3 3 3 3 3 3

+ + + + + + + + + 1+0 + +

z3 0 0 0 0 0 0 0 0 0 0 1 1 1

+ + + + + + + + + + + +

Figure 12. The unique state, the so called ground state, for the system Sµ,θ,w0 with µ = (2, 3, 0), 0 0 w = s1, θ = (1, 0, 0) with n = 2 and r = 3. Note that w is such that the colors of the top boundary is of the order w0P from left to right.

2.4. The R-matrix and the Yang-Baxter equation. The Yang-Baxter equation may be expressed as an identity of partition functions consisting of three vertices, two of which are of the kind discussed in the previous section denoted T , and one of a diﬀerent kind denoted R which involves only horizontal edges as illustrated in Figure 13. A solution to the (RTT) Yang-Baxter equation given sets of weights for the T -vertices, is a third set of weights for the R-vertices so that the partition functions agree. The Boltzmann weights at each vertex encode an endomorphism of a tensor product of vector spaces, so we may speak of a solution 16 to the Yang-Baxter equation as a set of Boltzmann weights or equivalently as encoded in an “R-matrix” describing the associated endomorphism, and we use these terms interchangeably. Our ultimate goal in this section is a solution to the Yang-Baxter equation for the fully-fused

models of the previous section using the weights at vertices Tzi and Tzj . We accomplish this by describing solutions to auxiliary Yang-Baxter-like equations for the monochrome vertices labeled T cw and T cw, described in the following theorem. By a zi zj similar argument as in Lemma 5.4 of [12], such auxiliary solutions guarantee a solution of the Yang-Baxter equation for the color-fused and fully-fused models.

β β

T cw T cw τ ν zi θ τ zj ψ θ

Rcw γ Rc0w0 zi,zj δ zi,zj

µ T cw ρ T cw ρ σ zj σ zi φ

α α

Figure 13. The auxiliary Yang-Baxter equations. These are miniature systems with ﬁxed boundary spins α,σ,τ,β,θ,ρ. The partition function is the sum over the interior spins µ, ν, γ, δ, φ, ψ which are allowed to vary. Here (c0, w0) is the successor to (c, w) deﬁned in (2.1). Theorem 2.6 asserts that the partition functions of these two systems are equal.

Theorem 2.6. Let c and w be a color and scolor, respectively. Let σ, τ, θ, ρ be elements of the horizontal spinset (2.2), and let α, β be spins of vertical type cw, that is, elements of the two-element set {+, cw}. Let (c0, w0) be the successor to (c, w) according to (2.1). Then using the Boltzmann weights for the vertices T cw from Figure4 and Rcw from Figure 14, zi zi,zj the partition functions of the two systems in Figure 13 are equal. The proof is postponed to later in this section. Note that the R-matrix is realized as a vertex that is attached to four edges of horizontal type, with the corresponding admissible spinsets consisting of a single attribute (one color or one scolor). Ultimately the proof of Theorem 2.6 may be reduced to a computation that can be checked using a computer. But it is not a priori clear how this can be checked for all n and r and for all choices of the function g in a finite manner. Thus our main task is to show how the proof can be reduced to a finite number of equations that need to be checked. We have taken g to be an arbitrary function of the integers modulo n, subject to Assumption 2.3. Now let us define a particular choice for g:

−v if a ≡ 0 mod n, (2.5) g (a) = √ 0 − v otherwise. √ Here v is a ﬁxed square root of v.

Lemma 2.7. Suppose Theorem 2.6 is true for g = g0. Then it is true for an arbitrary g satisfying Assumption 2.3. 17 a a b a b b x x Rcw Rcw Rcw Rcw zi,zj zi,zj zi,zj zi,zj a a a b a a x x n (1 − v)zi a < b < c n (1 − v)zi c 6 a < b n n n−1 n n v(zi − zj ) a < b (1 − v)zi zj a < c 6 b n n zi − vzj n n n−1 zj − vzi zi − zj a > b (1 − v)zj zi a > c > b n (1 − v)zj a > b > c n (1 − v)zj c > a > b y x y y a x x a Rcw Rcw Rcw Rcw zi,zj zi,zj zi,zj zi,zj x y x x x a a x

n zj y−x n n (1 − v)zi ( z ) y > x n n n n g(y − x)(z − z ) zi v(z − z ) z − z i j (1 − v)zn( j )y−x y < x i j i j j zi x x a a Rcw Rcw zi,zj zi,zj a a x x

n zj x−w n zi x−w (1 − v)z ( ) a < c, w x (1 − v)zj ( z ) a < c, w 6 x i zi 6 j n zj x−w n zi x−w (1 − v)z ( ) a < c, w > x (1 − v)zi ( ) a < c, w > x j zi zj n zj x−w−1 n zi x−w−1 (1 − v)zi ( z ) a > c, w < x (1 − v)zj ( z ) a > c, w < x zi j n j x−w−1 n zi x−w−1 (1 − v)zj ( ) a c, w x (1 − v)z ( ) a c, w x zi > > i zj > >

Figure 14. The R-matrix depending on a color c and charge w, as well as a pair of spectral parameters zi and zj. Here a and b are distinct colors and x and y are distinct scolors. We allow a or b to equal c, and w to equal x or y.

Proof. Write g(a) = γ(a)g0(a) where γ(a) = 1 if a ≡ 0 mod n and γ(a)γ(−a) = 1. Consider the scolored lines in both ﬁgures of Figure 13 that begin in one of the three vertices β, θ and ρ, and terminate in one of τ, σ, α. We ask what possible contributions of γ(a) with a 6≡ 0 mod n can occur. From the Boltzmann weights in Figures4 and 14 we see that such contributions occur when two scolored lines with distinct scolors cross. For both systems, such crossings can be predicted just from the order of the scolors of the boundary spins β, θ, ρ and τ, σ, α, with the following caveat: we cannot determine from this information whether two scolored lines cross twice or not at all. But if they cross twice, they produce factors γ(a) and γ(−a) for some a, and hence these factors cancel. Since modulo this caveat, the crossings of the scolored lines (with distinct scolors) is determined by the order of the boundary spins, the product of the factors γ that can occur are the same for the left and right systems in Figure 13.

Remark 2.8. We made the choice g = g0 in equation (2.5) in order to simplify computations. Ultimately, we will specialize g to certain quantities that have number theoretic meaning called Gauss sums. The R-matrix under this specialization (or any other that satisﬁes Assumption 2.3) will correspond to a Drinfeld twist of the Uq glb(r|n) R-matrix (see Proposition 2.11). Drinfeld twisting is a procedure that changes the comultiplication, antipode 18 and universal R-matrix of a given quantum group, and is particularly useful in constructing new solutions of the Yang-Baxter equation. It was introduced by Drinfeld [31] and further studied by Reshetikhin [45]. Proof of Theorem 2.6. To make the auxiliary Yang-Baxter equation of Figure 13 amenable to proof by checking a ﬁnite number of equations independent of n and r we make a change of basis for the vector spaces of horizontal edges used to the describe the various weights for T cw and Rcw in Figures4 and 14. z zi,zj Starting from a basis {vA} where A is an element of the horizontal spinset (2.2) this change of basis will perform a rescaling of each basis element that will depend not only on the spin A but also on the data zi, c and w of a vertex to the left of which the horizontal edge may w−x+1−n be attached. By choosing this rescaling to be zi if A = x is a scolor and 1 otherwise we obtain the Boltzmann weights of Figure 15 and the R-matrix remains unchanged except n for the weights shown in Figure 16. Here we have replaced zi with ζi and similarly for ζj, and we will be able to prove the equivalent corresponding Yang-Baxter equation for these weights while keeping ζi and ζj as independent parameters. Note that all horizontal edges for the R-matrix Rcw all attach to the left of a vertex (c, w) with color c and scolor w, while zi,zj for T cw the left and right horizontal edges attach to the left of vertices (c, w) and (c0, w0) zi respectively where (c0, w0) is the successor of (c, w) according to (2.1).

a1 a2 b1 + cw cw

cw cw cw y Tz y a Tz a y Tz y

+ cw cw  1 c < c 1 if c < cr v c > a  r z if c = cr, w > w0 z c = a g(y − w) × z c = cr, w > w0 z/ζ if c = c , w = w 1 c < a  r 0 z/ζ c = cr, w = w0 b2 c1 c2 + + cw

cw cw cw a Tz a c Tz w w Tz c

+ cw +  1 c < cr z if a = c,  (1 − v)ζ × z c = cr, w > w0 z/ζ 1 otherwise.  z/ζ c = cr, w = w0

Figure 15. The Boltzmann weights after a change of basis from Figure4. Here we have replaced n z with ζ. In the Yang-Baxter equation (z, ζ) will be replaced by either (zi, ζi) or (zj, ζj) where ζi n n and ζj will be regarded as a independent parameters in the proof separate from zi and zj .

Now there are two key observations. The ﬁrst observation is that at most three colors and at most three scolors can appear on the boundary and interior edges of the systems in Figure 13. Indeed, based on the conservation of colors and scolors in T and R we get that the colors of (β, σ, τ) are the same as the colors of (α, ρ, θ) in some permutation (also 19 y y x x a a Rcw Rcw Rcw zi,zj zi,zj zi,zj x x a a x x

2 −1 −1 2 −1 −1 (1 − v)ζi ζj zi zj a < c, w 6 x (1 − v)ζj ζi zj zi a < c, w 6 x −1 −1 (1 − v)ζi y > x (1 − v)ζizi zj a < c, w > x (1 − v)ζj zj zi a < c, w > x 2 −1 2 −1 (1 − v)ζj y < x (1 − v)ζi ζj a > c, w < x (1 − v)ζj ζi a > c, w < x (1 − v)ζi a > c, w > x (1 − v)ζj a > c, w > x

Figure 16. The modiﬁed weights for the R-matrix after a change of basis. Remaining weights n are unchanged from Figure 14 with the exception of replacing each occurrence of zi with ζi and ±1 ±1 similarly for ζj (such that the weights are all polynomial in ζi and ζj ). The parameters zi and ζi, as well as zj and ζj, will be treated as independent parameters in the proof.

counting non-colored edges). Similarly, we have that the scolors of (β, ρ, θ) are the same as those of (α, σ, τ) in some permutation. Besides the colors and scolors of the boundary edges in Figure 13 we also have the colors and scolors of the vertices which are (c, w) and (c0, w0). Thus, we can have at most five different colors and at most five different scolors appearing in the Yang-Baxter equation. The second observation is that (with g = g ) the Boltzmann weights for Rcw and T cw in 0 zi,zj zi Figures 15 and 16 depend only on the inequalities between the colors and scolors and not their actual values. This is the reason for the particular change of basis we made. Therefore if we check the auxiliary Yang-Baxter equation displayed in Figure 13 for these weights with r = 5 and n = 5, every possible relative ordering of color and scolor has occurred in these calculations. Note that for this argument, it is crucial that ζi and ζj are not linked n to zi and zj by the equation ζi = zi , even though this is what we later use to relate these weights to those of Figures4 and 14. Thus if the auxiliary Yang-Baxter equation with the change of basis is valid for n = r = 5, then it is true for all n and r, and so is the auxiliary Yang-Baxter equation in the statement of the theorem. The case n = r = 5 with the weights in Figures 15 and 16 was checked using a computer program (in Sage using symbolic manipulations), and so the theorem is proved. After repeated use of these auxiliary Yang-Baxter equations, passing through a complete cycle of pairs of colors and scolors (c, w), we immediately arrive at the following: Corollary 2.9. The fully-fused lattice model satisfies the Yang-Baxter equation, written as an identity of endomorphisms on V ⊗ V ⊗ W :

(Rzi,zj )12(Tzi )13(Tzj )23 = (Tzj )23(Tzi )13(Rzi,zj )12 where V is a vector space with basis indexed by the set of colors and scolors and W is a vector space with basis indexed by the power set of the set of pairs (c, w) of colors and scolors. We use the R-matrix weights in Figure 17 and the notation (A)`k means applying the endomorphism A to the `- and k-th factors in the tensor product. Remark 2.10. The set of color-scolor pairs (c, w) is in bijection with the set of odd negative √ roots of U1/ v(gl(r|n)). Therefore the power set that is the basis of W in Corollary 2.9 is a 20 basis of the exterior algebra on the space spanned by the odd negative roots. This observation may be the key to a quantum group interpretation of W .

a a b a b b x x

Rzi,zj Rzi,zj Rzi,zj Rzi,zj a a a b a a x x n n n n n v(zi − zj ) a < b (1 − v)zi a < b n n zi − vzj n n n zj − vzi zi − zj a > b (1 − v)zj a > b y x y y a x x a

Rzi,zj Rzi,zj Rzi,zj Rzi,zj x y x x x a a x

n zj y−x n n (1 − v)zi ( z ) y > x n n n n g(y − x)(z − z ) zi v(z − z ) z − z i j (1 − v)zn( j )y−x y < x i j i j j zi x x a a

Rzi,zj Rzi,zj a a x x z (1 − v)zn( j )x (1 − v)zn( zi )x i zi j zj

Figure 17. The R-matrix for the fully-fused model. Here a and b are distinct colors and x and y are distinct scolors. The weights are obtained from the special case of Figure 14 where c = c1 and w = wn−1.

The fully-fused R-matrix in Figure 17 is obtained from the unfused R-matrix in Figure 14 by taking the vertex color and scolor to be c = c1 and w = wn−1. We may similarly obtain a fully-fused R-matrix from the one in Figure 16, which is related to the one in Figure 17 by 1−n the change of basis introduced in the proof of Theorem 2.6. The extra factor of zi in this change of basis was made to facilitate the comparison with known R-matrices from quantum groups resulting in the following proposition.

Proposition 2.11. The fully-fused R-matrix in Figure 17 for the metaplectic Iwahori ice √ model is a Drinfeld twist of the U1/ v glb(r|n) R-matrix.

Proof. After the change of basis made in the proof of Theorem 2.6 we consider the fully-fused √ version of the R-matrix in Figure 16 and compare it to Kojima’s U1/ v glb(r|n) R-matrix in [35] acting on a tensor product of evaluation representations. In fact, Kojima considers √ the U1/ v slb(r|n) R-matrix, but the two quantum groups have the same R-matrix up to some normalization factor.√ We start from the R-matrix in [35] parametrized by q and z which we equate to 1/ v and ζi/ζj respectively, and then rescale each weight by the factor ζi − vζj. Note that Kojima considers a graded Yang-Baxter equation, while we are not, which means that weights where all edge conﬁgurations are colors are multiplied by −1. The resulting R-matrix then diﬀers from the one in Figure 16 by a Drinfeld twist parametrized by 21 a function φ(a, b) where a and b are horizontal edge states rescaling weights of the form

b a

Rzi,zj a b

The procedure of using such a function φ to perform a Drinfeld twist is explained in [9, Section 4] and [10, Section 4]. The function satisfies φ(a, b)φ(b, a) = 1 and is here given by √ −1/ v if a and b are colors and a > b,  √ − v if a and b are colors and a < b,  √ g(b − a)/ v if a and b are scolors and a 6= b, (2.6) φ(a, b) = √ v if a is a scolor and b is a color,  √ 1/ v if a is a color and b is a scolor,  1 otherwise. 2.5. Recurrence relations for the partition functions. The results of the previous section allow us to demonstrate that the partition functions satisfy the following recurrence. For x ∈ Z/nZ, let bxcn denote the least nonnegative residue of x mod n, that is, the integer such that 0 6 bxcn < n and x ≡ bxcn (mod n). Similarly, let dxen denote the integer such that 0 < dxen 6 n and x ≡ dxen (mod n). Then for x, y ∈ {0, 1, . . . , n − 1}, y − x if y > x (2.7) dy − xen = y − x + n if y 6 x. Proposition 2.12. Let µ ∈ r , θ ∈ ( /n )r, and w ∈ S defining a system S . Then, Z>0 Z Z r µ,θ,w for the simple reflection si, the following equality of partition functions hold

−dθi−θi+1enαi −nαi (2.8) (1 − v)z Z(Sµ,θ,w)(z) + g(θi − θi+1)(1 − z )Z(Sµ,siθ,w)(z) =

( −nαi (1 − v)Z(Sµ,θ,w)(siz) + (1 − z )Z(Sµ,θ,s w)(siz) if `(siw) > `(w) = i −nαi −nαi (1 − v)z Z(Sµ,θ,w)(siz) + v(1 − z )Z(Sµ,θ,siw)(siz) if `(siw) < `(w)

Proof. Let y = bθicn and x = bθi+1cn. Applying the Yang-Baxter equation in the fully-fused model in rows i and i + 1, we have the following identity (as pictured in Figure 18)

 y x   y y 

(2.9) wt  Rzi,zi+1  Z(S )(z) + wt  Rzi,zi+1  Z(S )(z) =   µ,siθ,w   µ,θ,w x y x x

 c d   d d 

wt  Rzi,zi+1  Z(S )(s z) + wt  Rzi,zi+1  Z(S )(s z).   µ,θ,siw i   µ,θ,w i d c c c Here wt(R) denotes the Boltzmann weight of the vertex R, and we assume x =6 y and c =6 d. If x = y the two terms on the left-hand side agree, but by the Yang-Baxter equation, only one of them should be included. Note that the colors in our model are distinct, meaning that c 6= d. 22 Inserting the Boltzmann weights from Figure 17, being careful to distinguish between the cases when x 6= y versus x = y, and using (2.7) proves the statement.

y y Tzi ··· Tzi d Tzi+1 ··· Tzi+1 d

Rzi,zi+1 = Rzi,zi+1

x Tzi+1 ··· Tzi+1 c x Tzi ··· Tzi c

Figure 18. Left: the system Sµ,θ,w with the R-matrix attached. Right: after using the Yang- Baxter equation.

Note that from the recursion relations from Proposition 2.12 and the base case from Lemma 2.5 we can compute the partition function for all systems Sµ,θ,w.

3. Metaplectic Iwahori Whittaker functions 3.1. Metaplectic preliminaries. We begin with a brief review of the construction of metaplectic groups and their unramified principal series representations. The discussion follows [10] whose treatment of metaplectic groups is based on the work of Matsumoto [37], Kazhdan and Patterson [34], Brylinski and Deligne [19], and McNamara [38, 39]. Because they may be stated uniformly, the constructions in this section will be described for covers associated to any split, reductive algebraic group G. In Section4, where we perform explicit calculations for comparison with the lattice models of the prior section, we will restrict to a particular cover of the general linear group. Let F be a non-archimedean local field and let o be its ring of integers with uniformizer $ and maximal prime ideal p = $o. Denote the cardinality of the residue field o/p by v−1. Fix a positive integer n. We shall work under the assumption that v−1 ≡ 1 mod 2n which implies that F contains the 2n-th roots of unity. Let µn denote the group of n-th roots of × unity in F and fix an embedding ι : µn → C . Let G be a connected reductive group defined and split over F . Let ∆∨ be the root system of G, which we assume to be irreducible. Remark 3.1. The root system ∆∨ is dual to the root system ∆ of the Langlands dual group Gb . The reason for this notational choice is that in this paper it is the representation theory of the dual group that will be more prevalent.

Let T be a fixed split maximal torus of G, and let Tb be the corresponding maximal torus of Gb . Then ∆∨ is a subset of the character group X∗(T). The dual root system ∆ ∗ is a subset of the cocharacter group X∗(T), which is identified with X (Tb); we will denote ∗ ∨ ∨ Λ = X∗(T) = X (Tb). If α ∈ ∆ , the corresponding element of ∆ is denoted α. We fix a pinning. This means that for every root α∨ ∈ ∆∨ there is a morphism 1 u iα∨ : SL2 → G, whose image we will denote Gα∨ . Let xα∨ (u) = iα∨ 1 . Then the image 23 of iα∨ is the subgroup generated by one-parameter unipotent subgroups Nα∨ on which T acts (through conjugation) by the character α∨ ∈ X∗(T), and

Gα∨ = hNα∨ , N−α∨ i.

We require the maps iα∨ to satisfy commutation relations as in Chevalley [24, Section III]. ∨ ∨ ∨ The root system ∆ may be partitioned into positive and negative roots, ∆+ and ∆−. These choices correspond to positive and negative Borel subgroups of G, denoted B and B− − − Q respectively. Then N(F )T(F ) = B(F ) and N (F )T(F ) = B (F ) where N := ∨ ∨ Nα∨ α ∈∆+ − Q and N := ∨ ∨ Nα∨ is its opposite. α ∈∆− Let K := G(o) be the group generated by the xα∨ (o), which is a maximal compact subgroup of G(F ). Let W = NG(T)/T denote the Weyl group of G with simple reflections ∨ si corresponding to each simple root αi . We may choose Weyl group representatives in K. A metaplectic n-fold cover of G(F ) is a central extension of G(F ) by µn: p (3.1) 1 −→ µn −→ Ge −−→ G(F ) −→ 1. 2 As a set Ge ≡ G(F ) × µn, with the group multiplication given by a cocycle σ ∈ H (G(F ), µn). The group K splits in Ge, so we fix a splitting and identify K with its image in Ge. Let −1 −1 Be = p B(F ) and Te = p T(F ) . As in [19, 38, 39], the cover (3.1) is associated with a W -invariant symmetric bilinear form B :Λ × Λ → Z whose restriction to the root lattice is required to be even. Thus the × ∼ quadratic form Q(µ) := B(µ, µ)/2 is integer-valued. If x ∈ F and λ ∈ Λ = X∗(T) let xλ be a representative of x under the cocharacter of T corresponding to λ. By abuse of notation, we will denote a representative in Te by the same notation xλ. The bilinear form B is characterized by λ µ B(λ,µ) × (3.2) [x , y ] = (x, y)n , x, y ∈ F , λ, µ ∈ Λ, λ µ where [x , y ] is the group commutator and (x, y)n ∈ µn is the n-th power Hilbert symbol. This identity holds regardless of the representatives xλ and yµ. Since we are assuming that λ µ F contains the 2n-th roots of unity, ($, $)n = 1. Thus, by (3.2), $ and $ commute and so we may choose the representatives $λ such that $λ$µ = $λ+µ. See Section 5 in [10] for a discussion on the bilinear form in relation to the metaplectic group. Let T (o) = Te ∩ K. Since the cover splits over K, the group T (o) may (like K) be regarded as a subgroup of either G(F ) or G. Denote by H := C T (o). It is a maximal abelian e Te subgroup of Te by Lemma 1 of [38]. λ The map λ 7→ $ induces an isomorphism Λ → Te /µnT (o). Let (3.3) Λ(n) := {µ ∈ Λ | B(y, µ) ≡ 0 mod n for all y ∈ Λ}. It follows from Lemma 1 in [38] (and its proof) that the subgroup H corresponds to Λ(n) ∼ under this map. In other words the isomorphism Λ −−→ Te /µnT (o) induces an isomorphism (n) ∼ Λ = H/µnT (o). Therefore (3.4) Λ/Λ(n) ∼= Te /H. We will denote by Γ a fixed choice of representatives of Λ/Λ(n). We next turn to the definition of the metaplectic L-group given in McNamara [38]. In that paper, like this one, it was assumed that G is split and that the ground field F contained the 24 group µ2n of 2n-th roots of unity. These assumptions simplify many technical issues. See Weissman [50] for a theory of the metaplectic L-group without these limitations. Identifying Λ with X∗(Tb), if z ∈ Tb(C) and λ ∈ Λ, denote by zλ ∈ C the image of z under the rational character λ. Let us denote by C[Λ] the Laurent polynomial ring spanned over C by the polynomial functions z 7→ zλ (λ ∈ Λ), or by abuse of notation, by zλ. We will denote by C(Λ) its field of fractions. The rings C[Λ] and C(Λ) may be identified with the ring of regular and rational functions, respectively, on Tb(C). Similarly we define C[Λ(n)] and C(Λ(n)). Let Tb 0 be an algebraic torus that is the quotient of Tb by the finite subgroup of z such that zλ = 1 for λ ∈ Λ(n). Thus C[Λ(n)] and C(Λ(n)) may be identified with the rings of regular and rational functions on Tb 0(C). We will denote by pˆ : Tb 0(C) → Tb(C) the canonical surjection. If α ∈ ∆, let n (3.5) n := n := n α gcdn, Q(α) i αi as in [39, 38]. Then let

(3.6) ∆n = {nαα | α ∈ ∆}. (n) ∗ 0 McNamara [38, Theorem 12] proves that ∆n is a subset of Λ = X (Tb ), and that it is a ∨ ∼ root system. Moreover let α ∈ X∗(Tb) = Hom(Λ, Z) be the coroot corresponding to α. Then ∨ 0 ∼ (n) McNamara proves that the image of α in X∗(Tb ) = Hom(Λ , Z) is a multiple of nα. Let ∨ −1 ∨ 0 ∆n = {nα α |α ∈ ∆} ⊂ X∗(Tb ). ∗ 0 0 ∨ Then (X (Tb ), ∆n,X∗(Tb ), ∆n) form a root datum (Springer [48]) and by [48, Theorem 2.9] they are associated with a reductive group Gb 0. This is the metaplectic L-group.

Remark 3.2. If ∆ is simply-laced, or if it is doubly-laced and n is odd, or if G = G2 and n is not a multiple of 3, then the Cartan types of Gb and Gb 0 are the same. Otherwise, they are dual, since long roots of ∆ correspond to short roots of ∆n. This dichotomy is reﬂected in the local Shimura correspondence [47].

Let us define the unramified principal series of Ge. A function f on H, Te or Ge is called genuine if f(εg) = ι(ε)f(g) for ε ∈ µn. A quasicharacter of H is called unramified if it is trivial on T (o). Unramified genuine quasicharacters of H may be parametrized by elements 0 0 of Tb (C). If ez ∈ Tb (C), let us pick z ∈ Tb(C) such that pˆ(z) = ez under the surjection 0 λ λ pˆ : Tb(C) −→ Tb (C). There is a unique unramified character χ = χz such that χ($ ) = z . (The character χz actually only depends on ez.) Given an unramified genuine quasicharacter χ on H let

Te i(χ) = IndH χ = {f : Te → C, f(ht) = χ(h)f(t) for all t ∈ Te , h ∈ H}. The group Te acts on i(χ) by right translation. It is a finite-dimensional Te-module that may be shown to be irreducible. We inflate i(χ) from T to B and induce to G to obtain I(χ) := IndGe i(χ). Explicitly, we e e e Be have 1/2 I(χ) = {smooth f : Ge → i(χ) | f(bg) = δ χ (b)f(g) for all b ∈ B,e g ∈ Ge}, 25 where δ is the modular quasicharacter. The group Ge acts by right translation on I(χ) denoted 0 by π. For χ = χz defined above we denote I(χz) by I(z). Although it only depends on ez ∈ Te we will not put the tilde into the notation.

Remark 3.3. To avoid complications we shall impose two conditions on the character χz. We assume that znαα 6= v±1 and znαα 6= 1 for α ∈ ∆. The first condition assures that the intertwining integrals (3.18) do not vanish, and that the unramified principal series is irreducible. Although we impose this assumption, the vanishing of the intertwining integral is precisely the tool needed to understand the reducibility of the principal series representations; see [34, Theorem I.2.9] and its proof. The second condition assures that the intertwining integrals do not have poles. This condition can be removed later using Proposition 3.14, where it is shown that the Demazure-Whittaker operators that define the recurrence formulas satisfied by Whittaker functions do not have poles, due to cancellation of denominators.

Te Lemma 3.4. The module i(χ) = IndH (χ) has a T (o)-ﬁxed vector v0 that is unique up to λ scalar multiple. If t ∈ H then t · v0 = χ(t)v0. The vectors $ · v0 as λ runs through a set of coset representatives for Λ/Λ(n) are a basis of i(χ).

Proof. A complete set of double coset representatives for H\Te /T (o) consists of $λ as λ runs through a set of coset representatives for Λ/Λ(n). By Mackey theory, M i(χ)|T (o) = ξλ, λ∈Λ/Λ(n)

λ −λ where ξλ is the character of T (o) defined by ξλ(t) = χ($ t$ ). Since χ is genuine, ξλ is the trivial character only if $λ lies in the centralizer H of T (o). Hence there is, up to scalar multiple, only one T (o)-fixed vector v0. Now another application of Mackey theory shows that i(χ)|H contains a copy of χ. Since χ|T (o) is trivial, this must be the span of v0, proving that t · v0 = χ(t)v0 for t ∈ H. z z Define Φ◦ to be the unique K-fixed vector in I(z) normalized so that Φ◦(k) = v0 when k ∈ K. It may be defined by z 1/2 Φ◦(bk) = δ (b)b · v0,

according to the metaplectic version of the Iwasawa decomposition Ge = BKe . Using v0 here guarantees that this is well-defined. Denote by J + and J − the positive and negative Iwahori subgroups of G(o). Thus if k = o/$o is the residue field, these are the preimages of B(k) and B−(k) respectively under reduction modulo $. We will work mostly with J − and so we will denote it as J. Let I(z)J be the space of Iwahori fixed vectors in I(z). The dimension of this space is |W | according to [39, Proposition 5.1]. Moreover, we have the decomposition Ge = tw∈W BwJe , where we have chosen a representative for w in K, which (by abuse of notation) we will also z J denote as w. This allows us to define a basis Φw ∈ I(z) for w ∈ W as follows: ( Φz(b) if w0 = w (3.7) Φz (bw0l) := ◦ b ∈ B,e w0 ∈ W, l ∈ J. w 0 otherwise,

P z z From this definition it is easy to see that w∈W Φw = Φ◦. 26 We next describe basics about Whittaker functionals on metaplectic principal series. Recall that N−(F ) may be embedded in Ge using the trivial section. In fact, later results in this section will require the following slightly stronger condition given as Proposition (b) of Appendix 1 of [40], which is attributed to Deligne. We refer the reader to [40] for the proof, which is stated for global fields but works similarly for local fields as noted in the final sentence of the appendix.

− Proposition 3.5 (Deligne). There exists a unique T(F )-equivariant section sF of N (F ) into Ge lifting the natural projection map p. Let ψ be a non-degenerate character of N−(F ) viewed as a subgroup of Ge. We assume ∨ that if α is a simple root, then the character ψ ◦ x−α∨ of F is trivial on o but no larger fractional ideal. The space of complex valued Whittaker functionals − (3.8) Wz := {f : I(z) → C, f(nv) = ψ(n)f(v) for all n ∈ N (F ), v ∈ I(z)} z µ has a natural basis Ωµ indexed by coset representatives of the form $ for Te /H. We will now deﬁne this basis explicitly. We begin by considering the map Λ → i(χ) deﬁned by −λ λ (3.9) λ 7→ z $ · v0. If λ ∈ Λ(n) then $λ ∈ H and so by Lemma 3.4 this map is constant on the cosets of Λ(n). The image of our set of representatives Γ of Λ/Λ(n) under the map (3.9) is a basis of i(χ). z ∗ Let Lγ = Lγ enumerated by γ ∈ Γ be the dual basis of i(χ) , so ( zλ if λ − γ ∈ Λ(n), (3.10) Lz($λ · v ) = γ 0 0 otherwise.

z Note that Lγ does not depend on the choice of coset representative γ. Thus we may enumerate z (n) the basis by Lµ with any µ ∈ Λ/Λ . Let Wχ : I(χ) → i(χ) be the i(χ)-valued Whittaker function Z (3.11) Wχ(f) = f(n)ψ(n)−1 dn . N−(F ) z We then define the basis {Ωµ} of Wz as z z χ (n) (3.12) Ωµ := Lµ ◦ W , µ ∈ Λ/Λ . The metaplectic Iwahori-Whittaker functions we consider are with respect to the contragredient representation I(z−1) and defined as 1/2 z−1 z−1 (n) (3.13) φµ,w(z; g) := δ (g)Ωµ π(g)Φw g ∈ G,e µ ∈ Λ/Λ . Remark 3.6. There are different normalizations of the Whittaker function in the literature. We use the same normalization for the Whittaker function as in [10]. However, the normalization for Lµ is different in [10], which compensates for the extra factor in equation (48) of the same paper. The normalization used in both [9] and [39] is different from ours. These Whittaker functions are invariant with respect to J acting on the right and are − (N (F ), ψ)-invariant on the left. They are also genuine. Since Ge = tw∈W BwJe and Be = 27 −λ 0 N(F )Te we then have that φµ,w(z; g) is determined by its values on g = $ w where λ ∈ Λ and w0 ∈ W . Our main goal in this section is the determination of these metaplectic Iwahori-Whittaker −λ 0 0 functions φµ,w(z; g) on all such elements g = $ w with λ ∈ Λ and w ∈ W . In the next subsection, we provide this answer in terms of recursively defined “metaplectic Demazure operators” — certain divided difference operators made with data from our metaplectic cover. In the remainder of this subsection, we prove two preliminary results on the vanishing of these functions at certain elements of the above form and evaluate φµ,w at such g in the special case where w0 = w. This latter result may be viewed as the base case for the recursive relations to follow. The proofs of these initial results are similar to the special case where the cover is trivial, that is, the linear algebraic group case. The proofs in this special case may be found in [12], but require the section on N−(F ) to be T(F )-equivariant as in Proposition 3.5.

0 −λ 0 0 Lemma 3.7. Let λ ∈ Λ and w ∈ W . Then φµ,w(z; $ w ) = 0 unless λ is w -almost dominant (as in Deﬁnition 2.1).

Proof. The proof follows [12, Lemma 3.5] and [23, Lemma 5.1]. Denote φµ,w(z; g) by φ(g). 0 Assume that λ is not w -almost dominant and let αi be a simple root which does not 1 satisfy (2.3). By the construction of the character ψ there exists a t ∈ $ o such ψ ◦x−αi (t) 6= 1. − Then, since u := x−αi (t) ∈ N (F ), (3.14) ψ(u)φ($−λw0) = φ(u$−λw0) = φ($−λw0j)

0−1 λ −λ 0 −hα ,λi 0−1 i where j = w $ u$ w = x−w αi ($ t). No cocycles appear in this multiplication because of the T(F )-equivariant section on N−(F ) and the splitting of the cover over K. We must now analyze two possible cases:

0−1 + −hαi,λi · w αi ∈ ∆ and by our assumption this means that hαi, λi < 0 which implies $ t ∈ o, 0−1 − −hαi,λi · w αi ∈ ∆ which by our assumption implies that hαi, λi < −1 and $ t ∈ p. In either case j ∈ J, and since φ is J-invariant, (3.14) gives that ψ(u)φ($−λw0) = φ($−λw0) which must therefore vanish. Proposition 3.8. Let µ ∈ Λ/Λ(n), w ∈ W and λ ∈ Λ a w-almost dominant weight. Then ( v`(w)zλ if µ + λ ∈ Λ(n) (3.15) φ (z; $−λw) = µ,w 0 otherwise, where `(w) denotes the length of a reduced expression for w. Proof. This proof proceeds similarly to Proposition 3.6 in [12], where the cover is trivial. By deﬁnition −λ 1/2 −λ z−1 −λ z−1 φµ,w(z; $ w) = δ ($ )Ωµ π($ w)Φw Z 1/2 −λ z−1 z−1 −λ −1 = δ ($ )Lµ Φw (n$ w)ψ(n) dn . N−(F ) We may apply the variable change n 7→ $−λn$λ to the above integral without introducing cocycles, according to the T(F )-equivariant section on N−(F ) from Proposition 3.5. This 28 λ z−1 −λ change of variables does produce a change in measure by δ($ ). Since Φw ($ g) = 1/2 −λ λ z−1 δ ($ )z Φw (g) we may write Z Z 1/2 −λ z−1 −λ −1 λ z−1 −λ λ −1 (3.16) δ ($ ) Φw (n$ w)ψ(n) dn = z Φw (nw)ψ($ n$ ) dn. N−(F ) N−(F ) The proof now proceeds as in Proposition 3.6 of [12] by use of the Iwahori factorization for −1 − Jw = wJw , which allows one to conclude that the integrand is supported on Nw where −1 + − Y x−α(o) if w α ∈ ∆ , (3.17) Nw = −1 − x−α(p) if w α ∈ ∆ . α∈∆+ z−1 To complete the proof, we show that the value of the integrand is Φ◦ (1) = v0 so (3.16) is λ − z−1 z−1 just z v0 times the volume of Nw . From (3.7) it follows that Φw (nw) = Φ◦ (1) since the argument of the former is in wJ. It remains to show that $−λn$λ is in the kernel of ψ. For this it is suﬃcient to show that if α = αi is a simple positive root then −λ λ − $ x−αi (t)$ ∈ N (o) −1 + where, using (3.17), we may assume that t ∈ o if w (αi) ∈ ∆ and t ∈ p otherwise. Now

−λ λ hλ,αii $ x−αi (t)$ = x−αi ($ t), and because λ is w-almost dominant, $hλ,αiit is in o. Again no cocycles appear from our − `(w) choice of splitting. The volume of Nw is v . From the deﬁnition (3.10) we then have ( `(w) λ (n) −1 −1 v z if µ + λ ∈ Λ φ (z; $−λw) = v`(w)Lz (zλv ) = v`(w)Lz ($−λ · v ) = µ,w µ 0 µ 0 0 otherwise.

3.2. Intertwining operators and recurrence relations. Define the intertwining operator z Aw : I(z) → I(wz) as Z z −1 (3.18) Aw(f)(g) := f(w ng)dn N(F )∩wN−(F )w−1 when the integral is absolutely convergent, and by analytic continuation otherwise. Because χz is unramified, this definition does not depend on the representative of the Weyl group element. The intertwining operator can be restricted to the space of Iwahori fixed vectors, and we can explicitly compute the action of this operator.

Lemma 3.9. Let si ∈ W be a simple reﬂection and w ∈ W . Then ( (1 − c (s z))Φsiz + Φsiz if `(s w) > `(w), z z αi i w siw i (3.19) As (Φw) = i (v − c (s z))Φsiz + vΦsiz if `(s w) < `(w), αi i w siw i where 1 − vzniαi (3.20) cα (z) = . i 1 − zniαi Moreover (3.21) AsizAz = c (z)c (z). si si αi −αi 29 This is a generalization of Casselman [22, Theorem 3.4]. It is worth noting that

1 − cαi (siz) = cαi (z) − v, so there are diﬀerent ways of writing (3.19).

Proof of the Lemma. The case `(siw) > `(w) of (3.19) follows from [41, Lemma 6.3] using 1−v 1 − cαi (siz) = 1 − c−αi (z) = 1−zniαi . The second case of (3.19) follows from the ﬁrst case by interchanging the roles of w and siw, after some algebra. For (3.21), see [39, Proposition 4.4].

The intertwining integrals induce maps Wz−1 → Wwz−1 that can be expanded as

wz z X z (3.22) Ωµ ◦ Aw = τµ,ν(w; z)Ων. ν∈Λ/Λ(n) where τµ,ν is called the Kazhdan-Patterson scattering matrix that will be given explicitly later. We will often apply it for a simple reflection w = si and will then drop the Weyl group element from the notation τµ,ν(si; z) = τµ,ν(z). Applying this twice to the composition Ωz AsizAz and using (3.21) gives µ si si X c (z)c (z) if λ ≡ µ mod Λ(n), (3.23) τ (s z)τ (z) = αi −αi µ,ν i ν,λ 0 otherwise. ν∈Λ/Λ(n) The Whittaker functions we study are with respect to the contragredient representation I(z−1), for which equation (3.19) becomes (3.24) −1 −1 ( siz −niαi siz −1 −1 1 (1 − v)Φ + (1 − z )Φ if `(siw) > `(w), Az (Φz ) = w siw si w −n α −1 −1 1 − z i i (1 − v)z−niαi Φsiz + v(1 − z−niαi )Φsiz if `(s w) < `(w). w siw i Let g(a), a ∈ Z be the Gauss sum as defined in [39, Section 13], although with a different normalization. If g˜(a) are the Gauss sums in [39], then g(a) = vg˜(a). For our purposes, it is important to note that the Gauss sums are periodic modulo n and satisfy (3.25) g(0) = −v, g(a)g(n − a) = v for a 6≡ 0 mod n which is Assumption 2.3. The following result was proved for GLr by Kazhdan and Patterson [34, Lemma I.3.3], applied to the metaplectic Casselman-Shalika formula in [27], and generalized by [39] to a much larger family of covers for reductive groups, including those we work with here.

Proposition 3.10 (Kazhdan-Patterson, McNamara). For a simple reﬂection si and µ, ν ∈ Λ/Λ(n),

−1 −1 X −1 (3.26) Ωsiz ◦ Az = τ (z−1)Ωz µ si µ,ν ν ν∈Λ/Λ(n) where the coeﬃcients τµ,ν can be written as 1 2 τµ,ν = τµ,ν + τµ,ν 30 1 (n) 2 (n) and where τ vanishes unless ν ∼ µ mod Λ and τ vanishes unless ν ∼ si(µ) + αi mod Λ . Moreover,

j B(αi,µ) k 1 − v − − αi 1 −1 Q(αi) n (3.27) τµ,µ(z ) = z i 1 − z−niαi where bxcn denotes the least positive residue of x mod n, and

(3.28) τ 2 (z−1) = gB(α , µ) − Q(α )zαi . siµ+αi,µ i i Proof. We refer the reader to McNamara [39, Theorem 13.1], but make several comments about necessary modifications in the present set-up. Our normalization for the Whittaker z−1 functional Ωµ is slightly different accounting for the different z-exponents in (3.27) and z−1 (3.28), and we also use a different normalization for Aw explaining an overall factor of −1 −1 −αi cαi (z ). Furthermore, we are applying this formula at z , so McNamara’s xα is our z . Thus l m 1 − v B(αi,µ) 1 −1 −ni n Q(α ) αi−si(µ)+µ (3.29) τµ,µ(z ) = z i i 1 − z−niαi where dxe denotes the smallest integer at least x. From the bilinearity and W -invariance B(αi,µ) B(αi,µ) of B we get that si(µ) − µ = − αi which also tells us that is an integer. Then, Q(αi) Q(αi) x using nd n e − x = b−xcn for x ∈ Z we obtain (3.27). Note that (3.28) may also be written as

(3.30) τ 2 (z−1) = gQ(α ) − B(α , µ)zαi . µ,siµ+αi i i

As seen from (3.26) the coeﬃcients τµ,ν are independent of the coset representatives for µ −1 −1 siz z and ν since Ωµ and Ων are, which is an advantage of our normalization. Furthermore, we see directly that the exponent in (3.27) is invariant under shifts µ → µ + λ with λ ∈ Λ(n) since B(αi, λ) ≡ 0 mod n which one can show is equivalent to B(αi, λ)/Q(αi) ≡ 0 mod ni. Combining the two formulas (3.24) and (3.26) we can prove the following functional equation for metaplectic Iwahori Whittaker functions.

Proposition 3.11. The following functional equation holds:

(3.31) τ 1 (z−1)φ (z; g) + τ 2 (z−1)φ (z; g) = µ,µ µ,w µ,si(µ)+αi si(µ)+αi,w

( −niαi 1 (1 − v)φµ,w(siz; g) + (1 − z )φµ,siw(siz) if `(siw) > `(w),

−niαi −niαi −niαi 1 − z (1 − v)z φµ,w(siz) + v(1 − z )φµ,siw(siz) if `(siw) < `(w).

Proof. The intertwining integral is a Ge-homomorphism and therefore

−1 −1 −1 −1 −1 −1 Ωsiz ◦ Az π(g)Φz = Ωsiz π(g)Az Φz . µ si w µ si w By the definition of φ in (3.13) the statement now follows by expanding the left- and right-hand sides using (3.24) and (3.22). 31 3.3. Metaplectic Demazure-Whittaker operators. Recall that C[Λ] is the ring of Laurent polynomials in zλ with λ ∈ Λ over C and that C(Λ) is the corresponding field of fractions. The results of the previous section allow us to compute all Iwahori Whittaker functions recursively. Let us fix g ∈ Ge and z ∈ Tb(C). For fixed w ∈ W , assemble the values γ φw := φγ,w(z; g)(γ ∈ Γ) into a vector γ |Γ| φw = φw(z; g) = (φw)γ∈Γ ∈ C(Λ) . |Γ| Theorem 3.12. Define operators Ti on C(Λ) by

−1 −1 −1 1 − v (3.32) Ti = Dn (z ) + siτ (z ) Dn (z ) = i i zniαi − 1 where, for f ∈ C(Λ)|Γ|, τ (z−1)f ∈ C(Λ)|Γ| with components µ (3.33) τ (z−1)f = τ 1 (z−1)f µ + τ 2 (z−1)f si(µ)+αi . µ,µ µ,si(µ)+αi

The operator Ti is invertible and 1 (3.34) T−1 = zniαi D (z−1) + s τ (z−1) . i v ni i Furthermore, ( Tiφw if `(siw) > `(w) (3.35) φsiw = −1 Ti φw if `(siw) < `(w). Proof. Both cases of (3.35) follow from Proposition 3.11. That the operator in (3.34) is −1 actually the inverse of Ti can be checked by direct computation using (3.23). However we will instead make use of a fact that we will prove below in Theorem 3.15, namely the Hecke 1 relation (Ti − v)(Ti + 1) = 0. This implies that Ti is invertible with inverse v (Ti + 1 − v). niα Together with the fact that Dni (z) + 1 − v = z Dni (z), this proves (3.34).

We will call the operators Ti vector metaplectic Demazure-Whittaker operators. Now let us consider how these operators allow us to compute arbitrary Iwahori Whittaker functions − on Ge. Since a complete set of coset representatives for N (F )\G/Jµe n consists of elements of −λ 0 0 the form g = $ w with λ ∈ Λ and w ∈ W we may restrict to such g. Then φw = 0 unless λ is w0-almost dominant. Assuming this, we have the following explicit algorithm.

0 0 −λ Corollary 3.13. Let w, w ∈ W and λ a w -almost dominant weight, and let g = $ w2. |Γ| γ Let φw = φw(z; g) ∈ C(Λ) as above whose components are given by φw = φγ,w(z; g). 0 Consider any path in the Weyl group from w to w given by simple reﬂections sik , . . . , si1 as 0 0 0 w → si1 w → · · · → sik ··· si1 w = w and let eij = 1 if sij is an ascent in this path in the

Bruhat order, and otherwise let eij = −1. Then, 0 e e (3.36) φ (z; $−λw0) = v`(w )T ik ··· T i1 L(zλ), w ik i1 λ |Γ| λ γ z−1 λ z−1 −λ where L(z ) ∈ C(Λ) with components L(z ) = Lγ (z v0) = Lγ ($ · v0).

−λ 0 `(w0) λ Proof. The value of φw0 (z; $ w ) is given explicitly by Proposition 3.8 as v L(z ). Then −λ 0 −λ 0 Theorem 3.12 shows that φw(z; $ w ) can be obtained recursively from φw0 (z; $ w ) by −1 applying operators Ti and Ti in the stated order. 32 Demazure operators, as the term is commonly understood, are “divided difference” operators −αi −1 −αi such as f 7→ (1−z ) f −z si(f) that at first glance appear to introduce denominators, but actually do not since the numerator is automatically divisible by the denominator. In our L −γ (n) next proposition we will show that the space γ∈Γ z C[Λ ] is Ti-invariant meaning that on this space of polynomials Ti does not introduce denominators. In particular, since the −λ 0 `(w0) λ L −γ (n) −λ 0 vector φw0 (z; $ w ) = v L(z ) is in γ∈Γ z C[Λ ], this means that φw(z; $ w ) = 0 e e v`(w )T ik ··· T i1 L(zλ) from (3.36) also is, as expected from (3.10) and the definition of ik i1 φµ,w(z; g) in (3.13).

(n) (n) −1 |Γ| Proposition 3.14. Let R be C[Λ ] or C(Λ ). The operators Ti and Ti on C(Λ) = L L −γ γ∈Γ C(Λ) descend to operators on γ∈Γ z R.

Even though our Whittaker functions are polynomial we will still need to consider Ti as L −γ (n) an operator on γ∈Γ z R with R = C(Λ ) in Theorem 3.16. Note that Ti is not R-linear because of the si that appears in the deﬁnition (3.32). The integrality proved in this result allows us to remove the assumption, imposed in Remark 3.3, that znαα 6= 1 for α ∈ ∆.

−1 −1 Proof. We will prove the statement for Ti; the proof for Ti follows similarly. Let D(z ) := −1 L −γ Dni (z ) and f ∈ γ∈Γ z R . By linearity it is enough to consider f with a single non- (n) zero component γ with f γ = z−γh where h ∈ C(Λ(n)) for R = C(Λ(n)) and h = zλ with (n) (n) (n) λ ∈ Λ for R = C[Λ ]. Then the only non-zero components of (Tif) are γ −1 γ 1 −1 γ (Tif) = D(z )f + siτγ,γ(z )f (3.37) (T f)siγ+αi = s τ 2 (z−1)f γ = g(B(α , γ) − Q(α ))z−(siγ+αi)s h. i i siγ+αi,γ i i i γ −γ Since sih ∈ R, it remains to show that (Tif) ∈ z R. Using (3.29), we get that 1 − v 1 − zkniαi s γ kniαi+si(γ)−γ −γ −γ i (3.38) (Tif) = 1 − z si z h = z (v − 1) h zniαi − 1 1 − zniαi l m B(αi,γ) (n) γ −γ where k = ∈ . For R = (Λ ) we see that (Tif) is clearly in z R. For R = niQ(αi) Z C (n) λ(n) (n) (n) (n) C[Λ ] we could assume that h = z . We have that siλ − λ = − B(αi, λ )/Q(αi) αi. (n) (n) Since B(αi, λ ) ≡ 0 mod n is equivalent to B(αi, λ )/Q(αi) ≡ 0 mod ni we get that (n) (n) (n) (n) siλ − λ = k niαi with k ∈ Z. Thus, in this case (n) kniαi (k+k )niαi 1 − z s (n) 1 − z (3.39) i h = zλ 1 − zniαi 1 − zniαi (n) which is a geometric series in C[Λ ].

The ﬁnite Iwahori Hecke algebra Hv is the C-algebra with generators Ti subject to the 2 quadratic relations Ti = (v − 1)Ti + v and the braid relations. The (extended) aﬃne Iwahori Hecke algebra Hev adds an abelian subalgebra spanned by generators ϑλ with λ ∈ Λ subject to ϑλϑµ = ϑλ+µ together with the Bernstein relations

ϑλ − ϑsiλ (3.40) ϑλTi − Tiϑsiλ = (v − 1) . 1 − ϑ−αi The term extended is used here since, unless Λ is the root lattice, this Hecke algebra is slightly larger than the Hecke algebra of the usual aﬃne Weyl group, which is a Coxeter group. 33 As a variant, we may replace Λ here by the lattice Λ(n). Then the Bernstein relation must be modiﬁed to read

ϑλ − ϑsiλ (3.41) ϑλTi − Tiϑsiλ = (v − 1) . 1 − ϑ−niαi (n) Let Hev be the corresponding (extended) affine metaplectic Hecke algebra. Recall that v−1 is the cardinality of the residue field. For the algebraic group G(F ), Iwahori J and Matsumoto [33] showed that Hev−1 acts on I(χ) by convolution. For metaplectic forms and representations, Savin [47] computed the corresponding genuine Hecke algebra which acts (n) by convolution in a similar manner, and showed that it is isomorphic to the algebra Hev−1 . (n) (See Remark 3.2.) This convolution action is complementary to the action of Hev that we (n) are considering here. For example, both the intertwining operator action of Hev and the (n) L J convolution action Hev−1 appear together in an action on w I(wz) in [10], and these actions commute with each other. The next result gives an action of this Hecke algebra on functions. It may be extended to an action of the affine Hecke algebra as in [10, Theorem 2.2].

(n) |Γ| Theorem 3.15. We have an action of the affine metaplectic Hecke algebra Hev on C(Λ) |Γ| in which the generator Ti acts by the operator Ti on C(Λ) and the generator ϑλ multiplies |Γ| −λ f ∈ C(Λ) by the constant z . Thus the Ti satisfy the braid relations and the quadratic Hecke relation (Ti − v)(Ti + 1) = 0. Proof. We will deduce this from [10, Theorems 1.1 and 1.2]. These depend on the choice of the finite Weyl group of a root system acting on a weight lattice Λ. We will take the root (n) system to be ∆n and the weight lattice to be Λ . To avoid confusion, we will write roots in 0 ∆n as αe. Thus αe = nαα which means more precisely that if ez ∈ Tb (C) and if z is chosen so 0 α nαα thatp ˆ(z) = ez in terms of the canonical surjectionp ˆ : Tb(C) −→ Tb (C), then eze = z . The theorems also require, for each ez in a W -invariant open subset of the torus (which for 0 us will be Tb ) a vector space M(ez). We also need, for each simple reflection si ∈ W , a linear z map Ai : M(ez) → M(siez). A general schema to produce Hecke algebra representations is given in [10], subject to a couple of assumptions.

siz z Assumption 1 of [10]. The composition Ai ◦ Ai is the

αi −αi 1 − veze 1 − vez e (3.42) α −α . 1 − ezei 1 − ez ei Assumption 2 of [10]. Suppose that mij is the order of sisj. Let wij = ··· sisjsi = ··· sjsisj be the longest element of the dihedral group generated by si and sj. (There are mij factors z on both sides.) Then the Ai braid in the sense that if mij is the order of sisj then

sj siz siz z sisj z sj z z ··· Ai Aj Ai = ··· Aj Ai Aj with mij factors on both sides as maps M(z) → M(wijz).

On these two assumptions, [10, Theorems 1.1 and 1.2] give a representation of Hev on L −1 −1 w M(wz ). (We have replaced z in [10] by z .) The operator Ti sends f ∈ M(wz) into −1 (wz)−1 −1 −1 Di(wz )f + Ai f ∈ M(wz ) ⊕ M(siwz ), 34 −1 1−v where Di(z ) = α . To apply this, we will use the Hecke algebra corresponding to the zi −1 metaplectic L-group, for which z is replaced by ez and αi by αei, so 1 − v 1 − v D (z−1) = = , i e α n α ezei − 1 z i i − 1 which is D (z−1) in the notation (3.32). We take M(z) = (Λ)|Γ|. The map Aez : (Λ)|Γ| → ni e C i C |Γ| −1 C(Λ) is given by the matrix τ(z ), where z ∈ Tb(C) such that p(z) = ez. We must verify that two Assumptions are satisﬁed. αi niαi Assumption 1 is satisﬁed by (3.23), on substituting eze = z into (3.42); Assumption 2 follows from the property of the intertwining operators, that Aw2z ◦ Az = Az provided w1 w2 w1Lw2 `(w1w2) = `(w1) + `(w2). Theorem 1.1 of [10] then gives an action on M(z) = w M(wz). Following the proof of [10, Theorem 2.2], we may extract the action (3.32) on functions f in a single copy of C(Λ)|Γ|. To this end, consider the map

|Γ| M |Γ| w C(Λ) −→ C(Λ) , f 7→ ( f)w∈W . w∈W

|Γ| This copy of C(Λ) is stable under the action of Hv, and the Ti induce the operators Ti given by (3.32). The extension to Hev follows from [10, Theorem 1.2].

3.4. Averaged metaplectic Demazure-Whittaker operators. In Theorem 3.12 we define metaplectic Demazure-Whittaker operators acting on a vector of metaplectic Whittaker |Γ| functions φw ∈ C(Λ) similar to those in [10]. Recall that the Whittaker model of I(z) has dimension equal to the cardinality of Γ in the metaplectic setting. Therefore it is natural to consider a vector valued Whittaker function, as this allows us to study all Whittaker functions simultaneously. Other literature, such as [26, 41, 39, 46] define metaplectic Demazure- Whittaker operators acting on a particular averaged metaplectic Whittaker function, which we |Γ| will show is the sum of the components of φw ∈ C(Λ) . These averaged Demazure-Whittaker operators are defined in the above literature by (implicitly) making use of the Chinta-Gunnells action introduced in [25], while our operators are derived purely from the representation theory of the non-archimedean group. Indeed, in our case the operators are defined from the recursion relation in Proposition 3.11, which was in turn obtained from the intertwining operators and the Kazhdan-Patterson scattering matrix τµ,ν. We will in this section relate the two different types of metaplectic Demazure-Whittaker operators and show that the Chinta-Gunnells action arises naturally from our vector Demazure-Whittaker operators. (n) Recall that C(Λ(n)) denotes the field of rational functions in zλ with λ(n) ∈ Λ(n) over C and denote by C(Λ(n))[Λ] the ring of polynomials in zλ with λ ∈ Λ over this field. From our chosen representatives Γ of Λ/Λ(n) we have an isomorphism

M −γ (n) ∼ (n) (3.43) ϕ : z C(Λ ) −−→ C(Λ )[Λ] γ∈Γ by taking the sum of components. 35 (n) Then Ti induces an operator Ti on C(Λ )[Λ] such that the following diagram commutes

L −γ (n) Ti L −γ (n) γ∈Γ z C(Λ ) γ∈Γ z C(Λ )

(3.44) ϕ ∼ ϕ ∼

Ti C(Λ(n))[Λ] C(Λ(n))[Λ].

(n) Clearly, the Ti also satisfy the braid and the quadratic Hecke relations of Hev . Recall that L −γ (n) Proposition 3.14 shows that Ti descends to an operator on γ∈Γ z C[Λ ]. Furthermore, ϕ restricts to an isomorphism

M −γ (n) ∼ z C[Λ ] −−→ C[Λ] γ∈Γ

Thus, Ti descends to an operator on C[Λ]. P Let φΣ,w(z; g) := ϕ(φw) = γ∈Γ φγ,w(z; g). This particular metaplectic Iwahori Whittaker function has been studied extensively in the literature including for example [39, 41, 46]. By Theorem 3.12 and (3.44) it follows that, for a simple reflection si, ( TiφΣ,w(z; g) if `(siw) > `(w) (3.45) φΣ,siw(z; g) = −1 Ti φΣ,w(z; g) if `(siw) < `(w). 0 0 0 With the same setup as in Corollary 3.13, with a path w → si1 w → · · · → sik ··· si1 w = w in the Weyl group, we obtain that 0 e e (3.46) φ (z; $−λw0) = v`(w ) ik ··· i1 zλ ∈ [Λ]. Σ,w Tik Ti1 C Chinta and Gunnells [25] defined an action of the Weyl group on (in our notation) C(Λ(n))[Λ] that has appeared in subsequent works such as [26, 41, 44, 27, 39, 46]. Furthermore, Chinta, Gunnells, Puskásand Patnaik in several papers [26, 44, 41] study metaplectic Iwahori Whittaker functions through operators that they call Demazure-Lusztig operators, though in this paper we use the term Demazure-Whittaker operator. The Chinta-Gunnells Weyl group action, and related Hecke actions are studied in Sahi, Stokman and Venkateswaran [46], where they are given a conceptual basis related to the double affine Hecke algebra. We will now relate these works to the Hecke algebra action Ti. Specifically, we will show in our next Theorem 3.16 that the operator Ti can be expressed in terms of the Chinta-Gunnells action of [25] and that it matches the metaplectic Demazure operators defined by Chinta- Gunnells-Puskásin [26] and further studied by Patnaik-Puskásin [41]. Comparing (3.46) −λ 0 with [41, Corollary 5.4], this means that φΣ,w(z; $ ), with w = 1, equals the particular metaplectic Iwahori Whittaker function studied in Patnaik-Puskás[41] up to normalization.

From the isomorphism ϕ it is possible to recover the vector φw from φΣ,w. Note however, that our Corollary 3.13, or indeed (3.46), allows us to determine φΣ,w(z; g) for all values of g, and not only for g = $−λ. (n) For a simple reflection si ∈ W the Chinta-Gunnells action on C(Λ )[Λ] is linearly extended λ (n) from the action si ? z h where h ∈ C(Λ ) and λ ∈ Λ which is defined as (3.47) siλ j B(λ,αi) k (s h)z ) αi λ i Q(αi) n (ni−1)αi −nα αi si ? z h := (1 − v)z i − g(−Q(αi) − B(λ, αi))z (1 − z i ) . 1 − vz−niαi 36 To compare with equation (4.7) of [41], note that our Gauss sums are normalized differently: if gm is the Gauss sum in [41], then vgm = g(−m). Recall from (3.32) and (3.20) that

−niαi −1 −1 1 − v −1 1 − vz (3.48) D(z ) := Dn (z ) = cα (z ) = . i zniαi − 1 i 1 − z−niαi (n) Theorem 3.16. The operators Ti on C(Λ )[Λ] deﬁned from Ti by (3.44) are given by −1 −niαi −1 (3.49) Ti : f 7→ D(z )f − z cαi (z )si ? f, and match the operators in [26, 41], there called Demazure-Lusztig operators.

Proof. By linearity it is enough to consider f = zλh with λ ∈ Λ and h ∈ C(Λ(n)). Let (n) −1 L −γ (n) γ ∈ Γ such that −γ ≡ λ mod Λ and let f := ϕ (f) ∈ γ∈Γ z C(Λ ) . Then the only non-zero component of f is f γ = zλh. Similar to the computation in (3.37), but here 1 using (3.27) for τγ,γ instead of (3.29) the only non-zero components of (Tif) are

j B(αi,γ) k (1 − v) − αi γ −1 γ 1 −1 γ −1 Q(αi) n λ (Tif) = D(z )f + siτγ,γ(z )f = D(z ) + z i si z h (3.50) 1 − zniαi (T f)siγ+αi = s τ 2 (z−1)f γ = g(B(α , γ) − Q(α ))z−αi s zλh. i i siγ+αi,γ i i i Thus, since b− B(αi,γ) c = b B(αi,λ) c and g(B(α , γ) − Q(α )) = g(−B(α , λ) − Q(α )) we have Q(αi) ni Q(αi) ni i i i i that Tif = ϕ(Tif) equals ! j B(αi,λ) k 1 αi −1 Q(αi) n niαi −αi λ D(z ) + (1 − v)z i + g(−B(αi, λ) − Q(αi))(1 − z )z si z h. 1 − zniαi

−n α 1−vz i i −niαi −1 niαi −αi −niαi (ni−1)αi Since 1−zniαi = −z cαi (z ) and (1 − z )z = −(1 − z )z we conclude that

−1 −niαi −1 (3.51) Tif = D(z )f − z cαi (z )si ? f. This expression agrees with (4.10) of [41].

We have thus shown that the averaged Demazure-Whittaker operators Ti that arise from the component sum of the vector Demazure-Whittaker operators Ti, as well as the particular metaplectic Iwahori Whittaker function φΣ,w agree with those studied in [39, 41, 46]. Our equation (3.46) is [41, Corollary 5.4] and the fact that Ti descends to an operator on C[Λ] is [41, Claim 4.4]. The braid relations for Ti follow from those for Ti established in Theorem 3.15 which is based on [10]. Independent proofs of this fact can be found in [26] and [41, Appendix B], with a more conceptual proof in [46] based on a general construction of Weyl group and Hecke algebra representations.

4. Main theorem and conclusions The main theorem of this paper shows a relation between partition functions for the metaplectic Iwahori ice model described in Section2 and the metaplectic Iwahori Whittaker functions for G = GLr described for general G in Section3. Before we can state the main theorem, TheoremA, in Section 4.2 we need to specify the particular metaplectic cover we consider as well as some notation for G = GLr in the upcoming section. 37 4.1. Metaplectic Iwahori Whittaker functions for G = GLr. In this section we apply ∼ r the theory from the previous sections to the case G = GLr. We identify Λ = X∗(T ) = Z and use the standard elementary basis ei, such that αi = ei − ei+1. We choose a cover Ge such that the bilinear form B, which determines the group multiplication on Te by (3.2), becomes

B(ei, ej) := δij. which implies that Q(αi) = 1 and ni = n. Such a cover can be obtained by a particular cocycle as explained in [10, Section 5] and involves the 2n-power Hilbert symbol, which is one reason we assume that F contains the 2n-th roots of unity µ2n. (n) ∼ (n) ∼ r With this bilinear form Λ = nΛ and thus Λ/Λ = (Z/nZ) . We may choose the set Γ of coset representatives to be the set of weights µ ∈ Λ such that ρ − µ ∈ {0, 1, . . . , n − 1}×r where ρ = (r − 1, r − 2,..., 0). Letting the least nonnegative integer residue b·cn modulo n act elementwise we may write this as

r (4.1) µ = ρ − bθcn, θ ∈ (Z/nZ) .

Note that then B(αi, µ) = bθi+1cn − bθicn − 1.

We will denote the Whittaker function φρ−bθcn,w(z; g) simply by φθ,w(z; g). A similar z convention will be used for Ωθ and τθ,θ0 . In this setting, Proposition 3.8, which computes −λ 0 0 φθ,w(z; $ w ) when w = w , becomes the following

Corollary 4.1 (Proposition 3.8). For G = GLr, where we choose the representatives Γ of (n) r Λ/Λ as µ = ρ − bθcn with θ ∈ (Z/nZ) we get that ( v`(w)zλ if θ ≡ λ + ρ mod n (4.2) φ (z; $−λw) = θ,w 0 otherwise.

We will now rewrite the expressions for τ more explicitly. With µ = ρ − bθcn we note for (3.28) that si(ρ − bθcn) + αi = ρ − sibθcn. Furthermore, we have that x ∈ Z/nZ, 1 b−xcn = d1 − xen − 1, where dxen is the least strictly positive integer modulo n. Hence, τ and τ 2 coeﬃcients can be reduced to the following.

r Corollary 4.2 (Proposition 3.10). For G = GLr and θ ∈ (Z/nZ) with components θi the τ 1 and τ 2 coeﬃcients of (3.27) and (3.28) become 1 − v 1 αi −dθi−θi+1enαi τθ,θ = z z (4.3) 1 − z−nαi τ 2 = g(θ − θ )zαi . θ,siθ i i+1 r Corollary 4.3 (Proposition 3.11). For G = GLr and θ ∈ (Z/nZ) with components θi and w ∈ W the following functional equation holds: αi −dθi−θi+1enαi −nαi (4.4) z (1 − v)z φθ,w(z; g) + g(θi − θi+1)(1 − z )φsiθ,w(z; g) =

( −nαi (1 − v)φθ,w(siz; g) + (1 − z )φθ,s w(siz) if `(siw) > `(w), = i −nαi −nαi (1 − v)z φθ,w(siz) + v(1 − z )φθ,siw(siz) if `(siw) < `(w). 38 nr Lastly, let f ∈ (C(nΛ)[Λ]) with components enumerated by θ. The operator Ti defined in (3.32) reduces to 1 θ (dθi−θi+1en−1)αi θ −αi nαi siθ (4.5) (Tif) = (1 − v) z si − 1 f + g(θi − θi+1)z (1 − z )sif . 1 − znαi 4.2. The main theorem. We will now combine the results of Sections2 and 4.1 for the main theorem of this paper. Theorem A. Let w ∈ W , θ ∈ ( /n )r and let µ ∈ r . Let w0 ∈ W such that the order of Z Z Z>0 0 the colors for the top boundary edges for the system Sµ,θ,w in the monochrome picture is w P 0 from left to right where P = (cr, cr−1, . . . , c1). Then, the condition on µ and w is equivalent to λ := w0µ − ρ being w0-almost dominant and, with the conventions of Section 4.1, ρ −λ 0 (4.6) Z Sµ,θ,w (z) = z φθ,w(z, $ w ). For the definition of w0-almost dominant see Definition 2.1. Note that the partition functions determine all values for the metaplectic Whittaker function at any g ∈ Ge and that the set of boundary conditions maps bijectively to the set of data that determine all values of the metaplectic Iwahori Whittaker functions as pictured in Figure 19. This combinatorial r data consists of θ ∈ (Z/nZ) that determines the Whittaker functional Ωθ, w ∈ W which 0 0 determines the Iwahori fixed vector Φw, and the pairs (w , λ) where λ ∈ Λ and w ∈ W such that λ is w0-almost dominant which determines the argument g = $−λw0. Indeed, recall that φθ,w(z; g) is a genuine function on Ge, invariant under the Iwahori subgroup J on the right, − (N (F ), ψ)-invariant on the left and that the center Z(Ge) acts by scalars. Thus φθ,w(z; g) − is determined by its values on a set of double coset representatives for N (F )Z(Ge)\G/Jµe n which we may choose as g = $−λw0 with λ ∈ Λ and w0 ∈ W . That λ is w0-almost dominant follows from Lemma 3.7. Using the center, we may assume that λ + ρ, and therefore also µ has only non-negative parts.

 (w0, λ)     θ w  z1 ρ ρ 1/2 z−1 z−1   z φ (z; g) = z δ (g)Ω π(g)Φ = Z   = Z S 0−1 θ,w θ w  z2  w (λ+ρ),θ,w    z3    + + + + +

Figure 19. Dictionary between boundary data for lattice model and data determining the values of the metaplectic Iwahori Whittaker function illustrating the equality of TheoremA.

Proof. That the condition for µ and w0 is equivalent to λ := w0µ−ρ being w0-almost dominant follows from Remark 2.2. We prove the remaining statement by induction on the length of ww0−1. The base case w = w0 follows immediately from Lemma 2.5 and Corollary 4.1.

0 Using Theorem 3.12 we may determine φθ,siw from φθ ,w, and thus it is enough to verify that −ρ Zθ,w(z) := z Z Sµ,θ,w (z) satisﬁes the same recursion relations. In fact, it is enough to 39 compare with the recursion relations for φ in Corollary 4.3 which determine the operator Ti in (4.5) for GLr. Indeed, from the recursion relations for Z Sµ,θ,w (z) in Proposition 2.12 obtained from the Yang-Baxter equation, we get that

−dθi−θi+1enαi −nαi (4.7) (1 − v)z Zθ,w(z) + g(θi − θi+1)(1 − z )Zsiθ,w(z) =

( −nαi (1 − v)Zθ,w(siz) + (1 − z )Zθ,siw(siz) if `(siw) > `(w) = z−αi −nαi −nαi (1 − v)z Zθ,w(siz) + v(1 − z )Zθ,siw(siz) if `(siw) < `(w) which matches those for φ in Corollary 4.3. The factor z−αi comes from the fact that ρ −αi ρ ρ ρ−αi Z Sµ,θ,w (siz) = si z Zθ,w(z) = z z Zθ,w(siz) since siz = z . 4.3. R-matrices and p-adic intertwiners. In this section we will present a dictionary relating the following topics: p-adic intertwiners ←→ quantum group R-matrices ←→ lattice models R-matrices. In Proposition 2.11 we connected the R-matrix for the lattice model introduced in this paper to the Uq(glb(r|n)) R-matrix. We shall now refine this connection and show how R- matrices coming from quantum affine groups relate to non-archimedean local field intertwining integrals. √ × Let q = 1/ v. For fixed z ∈ C , let Vr,n(z) be an r + n dimensional complex vector space + − with basis given by elements vi (z) and vj (z), where 1 6 i 6 r indexes the set of colors and 1 6 j 6 n indexes the set of scolors. The weights of the R-matrix in Figure 16 can be used to create a map

Rq(z1, z2): Vr,n(z1) ⊗ Vr,n(z2) → Vr,n(z1) ⊗ Vr,n(z2) which satisﬁes the ungraded Yang-Baxter equation:

Rq(z1, z2)1,2Rq(z1, z3)1,3Rq(z2, z3)2,3 = Rq(z2, z3)2,3Rq(z1, z3)1,3Rq(z1, z2)1,2. The above should be thought of as an equality in End(Vr,n(z1) ⊗ Vr,n(z2) ⊗ Vr,n(z3)), and the subscripts i,j mean that we apply the map to the (i, j) entry of the tensor product. In equations (2.3)-(2.8) of [35], Kojima computes the Uq(glb(r|n)) R-matrix : Vr,n(z1) ⊗ Vr,n(z2) → Vr,n(z1) ⊗ Vr,n(z2), where Vr,n(z) is the evaluation module of Uq(glb(r|n)) of dimension (r|n), i.e. its even part has dimension r and its odd part has dimension n. This R-matrix satisfies the graded Yang-Baxter equation (see equation (2.11) in [35]). We shall denote by (r|n) −1 Rq (z1z2 ): Vr,n(z1) ⊗ Vr,n(z2) → Vr,n(z1) ⊗ Vr,n(z2) the (ungraded) Drinfeld twist of Kojima’s R-matrix described in the proof of Proposition 2.11. ∼ It follows by the same proposition that, after identifying Vr,n(z) = Vr,n(z), the two R-matrices are the same. A Drinfeld twist is a process which introduces extra complex factors in the (r|n) −1 R-matrix. Here we use a particular Drinfeld twist to introduce Gauss sums in Rq (z1z2 ) which otherwise do not appear naturally in the quantization of a universal enveloping algebra. + − Define Vr(z) and Vn(z) to be the subspaces of Vr,n(z) with basis vi (z) and vj (z), respectively. We define Vr(z) and Vn(z) in a similar way, and let Vr(z) := Vr(z1)⊗· · ·⊗Vr(zr), (r|n) −1 and Vn(z) := Vn(z1) ⊗ · · · ⊗ Vn(zr). Denote the restriction of Rq (z1z2 ) to Vr(z1) ⊗ Vr(z2) r −1 r −1 by Rq(z1z2 ). It is a well known fact in quantum group theory that Rq(z1z2 ) is then the 40 Uq(glb(r)) R-matrix (up to a Drinfeld twist), where one can now think of Vr(z) as the evaluation (r|n) −1 module for Uq(glb(r)). Moreover, restricting Rq (z1z2 ) to Vn(z1) ⊗ Vn(z2), produces the n −1 Uq−1 (glb(n)) R-matrix (again up to a Drinfeld twist), which we will denote by Rq−1 (z1z2 ). Recall the space of Whittaker functionals Wz defined in (3.8) and the space of Iwahori fixed vectors in the unramified principal series I(z)J from (3.7). Intertwiners induce maps on these spaces. It was shown in [9, Theorem 1.1] and in [12, Theorem 10.5, Remark 10.6] that there is a vector space embedding ξz and a vector space isomorphism ηz J n ∼ n ξz : I(z) ,−→ Vr(z ), ηz : Wz −−→ Vn(z ) ∼ = ∼ (4.8) = n n Vr(z ) Vn(z ) such that the following diagrams commute: non-archimedean quantum group lattice model

J ξz n ∼ n I(z) Vr(z ) Vr(z )

Az r nαi r nαi si τRq (z )i,i+1 τRq(z )i,i+1 ξ J siz n ∼ n (4.9) I(siz) Vr(siz ) Vr(siz )

ηz n ∼ n Wz Vn(z ) Vn(z ) Az τRn (znαi ) τ n (znαi ) si q−1 i,i+1 Rq−1 i,i+1 η siz n ∼ n Wsiz Vn(siz ) Vn(siz ) Here τ is the flip map, that is τ : v ⊗ w 7→ w ⊗ v, not to be confused with the scattering matrix of Section 3.2. Recall from the pictorial description of our main result in Figure 19 the natural correspondence between the boundary data for the lattice model and the data determining the values of the Whittaker function. We will now explain, by using this pictorial description, the natural correspondence between the two methods of deducing functional equations for the partition functions on the one side (as in Section 2.5) and for the Whittaker functions on the other (as in Section 3.2). On the lattice model side, we use the Yang-Baxter equation which allows us to move an R-vertex (essentially the R-matrix) from right to left. On the non-archimedean representation theory side one uses the fact that the intertwiner is a G˜-homomorphism, which allows us to move it from interacting with the Iwahori fixed vector z−1 z−1 to interacting with the Whittaker functional. Let AΩθ := Ωθ ◦ A. The two processes are pictured and compared in Figure 20. Note that in the non-archimedean representation theory side we only see the Uq(glb(r)) and Uq−1 (glb(n)) R-matrices, which appear as scattering matrices of intertwining operators on Iwahori fixed vectors and Whittaker functionals, respectively. On the lattice model side, one needs the full Uq(glb(r|n)) super R-matrix to solve the Yang-Baxter equations for the model. Yet this super R-matrix does not have any obvious role in the p-adic theory; only the gl(r) and gl(n) pieces do. It would be interesting to understand if the full super quantum group also has a direct role in the non-archimedean representation theory. 41 −1 −1 −1 −1 −1 −1 Az Ωsiz π(g)Φz = Ωsiz π(g) Az Φz si θ w θ si w

n r Rq−1 (r|n) Rq Rq

= ... = = ... =

+ + + + + + + + + + + +

Figure 20. Comparison between the two process of obtaining functional equations for the partition functions of the lattice model and the metaplectic Iwahori Whittaker functions. Different r restrictions of the R-matrices are identified with different actions of the intertwiners where Rq is an n (r|n) R-matrix for Uq(glb(r)), Rq−1 an R-matrix for Uq−1 (glb(n)), and Rq an R-matrix for Uq(glb(r|n)).

The Uq−1 (glb(n)) R-matrix also makes an appearance if we consider an alternative way of proving the braid relations for our vector metaplectic Demazure-Whittaker operators Ti. Recall that in the proof of Theorem 3.15 we use that τ of (3.33) satisﬁes similar braid relations which is Assumption 2 of [10]. As illustrated above, this Kazhdan-Patterson scattering matrix

appears as the Uq−1 (glb(n)) R-matrix on the lattice model side for G = GLr. The braid relations are then translated to the Yang-Baxter equation illustrated in Figure 21 involving three R-vertices (which we shall call the RRR Yang-Baxter equation). Furthermore, while the proof of Theorem 3.15 used the fact that the function g(a) appearing in τ is actually a Gauss sum coming from the metaplectic representation theory, the RRR Yang-Baxter equation holds for any function g(a) satisfying Assumption 2.3, as the R-matrix

comes from the Drinfeld twist of Uq−1 (glb(n)). Thus, at least for GLr, the Demazure-Whittaker operators Ti may be deﬁned for these more general parameters g(a) as well, in accord with [46] where they work with the scalar Demazure-Whittaker operators in similar generality.

zk zi zk zi R R R zj zj = zj zj R R R zi zk zi zk

Figure 21. The RRR Yang-Baxter equation with three R vertices that can be used to prove the braid relations for the vector metaplectic Demazure-Whittaker operators Ti when G = GLr.

Remark 4.4. Sahi, Stokman and Venkateswaran [46] construct the metaplectic Demazure- Whittaker operators Ti in (3.49) (and therefore the Chinta-Gunnells action) without making use of the representation theory of the p-adic group and for general parameters g(a). They consider quotients of a representation of affine Hecke algebras induced from the finite Hecke algebra. The latter is built by using (finite) R-matrices (see [46, Lemma 3.1]) which give rise to affine R-matrices via Baxterization. 42 The operators Ti that are studied in this paper and [26, 41] using non-archimedean methods can be seen to agree with the metaplectic Demazure-Whittaker operators of [46] via direct comparison in terms of the Chinta-Gunnells action as in Theorem 3.16. However, it would be desirable to give a conceptual explanation why the construction of [46] (which does not use the p-adic group) produce the same operators. Given the appearances of affine R-matrices in both [46] and the current paper, it would be worthwile to understand if their method of constructing quotients of induced Hecke representations has an analogue in the non-archimedean world.

References [1] A. Aggarwal. Current fluctuations of the stationary ASEP and six-vertex model. Duke Math. J., 167(2):269–384, 2018. [2] A. Aggarwal, A. Borodin, and M. Wheeler. Colored fermionic vertex models and symmetric functions, 2021, arXiv:2101.01605. [3] R. J. Baxter. One-dimensional anisotropic Heisenberg chain. Ann. Physics, 70:323–337, 1972. [4] R. J. Baxter. Exactly solved models in statistical mechanics. Academic Press Inc. [Harcourt Brace Jovanovich Publishers], London, 1982. [5] V. V. Bazhanov and A. G. Shadrikov. Trigonometric solutions of the triangle equations, and simple Lie superalgebras. Teoret. Mat. Fiz., 73(3):402–419, 1987. [6] A. A. Belavin, A. M. Polyakov, and A. B. Zamolodchikov. Infinite conformal symmetry of critical fluctuations in two dimensions. J. Statist. Phys., 34(5-6):763–774, 1984. [7] A. Borodin and L. Petrov. Integrable probability: stochastic vertex models and symmetric functions. In Stochastic processes and random matrices, pages 26–131. Oxford Univ. Press, Oxford, 2017. [8] A. Borodin and M. Wheeler. Coloured stochastic vertex models and their spectral theory, 2018, arXiv:1808.01866. [9] B. Brubaker, V. Buciumas, and D. Bump. A Yang-Baxter equation for metaplectic ice. Appendix ([11]) joint with Nathan Gray. Commun. Number Theory Phys., 13(1):101–148, 2019. [10] B. Brubaker, V. Buciumas, D. Bump, and S. Friedberg. Hecke modules from metaplectic ice. Selecta Math. (N.S.), 24(3):2523–2570, 2018. [11] B. Brubaker, V. Buciumas, D. Bump, and N. Gray. Duality for metaplectic ice (appendix to [9]), 2017, arXiv:1709.06500. [12] B. Brubaker, V. Buciumas, D. Bump, and H. P. A. Gustafsson. Colored vertex models and Iwahori Whittaker functions, 2019, arXiv:1906.04140. [13] B. Brubaker, V. Buciumas, D. Bump, and H. P. A. Gustafsson. Vertex operators, solvable lattice models and metaplectic Whittaker functions. Comm. Math. Phys., 380(2):535–579, 2020. [14] B. Brubaker, V. Buciumas, D. Bump, and H. P. A. Gustafsson. Colored five-vertex models and Demazure atoms. J. Combin. Theory Ser. A, 178, 2021. [15] B. Brubaker, D. Bump, G. Chinta, S. Friedberg, and P. E. Gunnells. Metaplectic ice. In Multiple Dirichlet series, L-functions and automorphic forms, volume 300 of Progr. Math., pages 65–92. Birkhäuser/Springer, New York, 2012. [16] B. Brubaker, D. Bump, and S. Friedberg. Schur polynomials and the Yang-Baxter equation. Comm. Math. Phys., 308(2):281–301, 2011. [17] B. Brubaker, D. Bump, and S. Friedberg. Weyl group multiple Dirichlet series: type A combinatorial theory, volume 175 of Annals of Mathematics Studies. Princeton University Press, Princeton, NJ, 2011. [18] B. Brubaker, D. Bump, and A. Licata. Whittaker functions and Demazure operators. J. Number Theory, 146:41–68, 2015. [19] J.-L. Brylinski and P. Deligne. Central extensions of reductive groups by K2. Publ. Math. Inst. Hautes Etudes´ Sci., (94):5–85, 2001. [20] V. Buciumas, T. Scrimshaw, and K. Weber. Colored five-vertex models and Lascoux polynomials and atoms. J. Lond. Math. Soc., 102(3):1047–1066, 2020. 43 [21] D. Bump, P. J. McNamara, and M. Nakasuji. Factorial Schur functions and the Yang-Baxter equation. Comment. Math. Univ. St. Pauli, 63(1-2):23–45, 2014, arXiv:1108.3087. [22] W. Casselman. The unramified principal series of p-adic groups. I. The spherical function. Compositio Math., 40(3):387–406, 1980. [23] W. Casselman and J. Shalika. The unramified principal series of p-adic groups. II. The Whittaker function. Compositio Math., 41(2):207–231, 1980. [24] C. Chevalley. Sur certains groupes simples. Tohoku Math. J. (2), 7:14–66, 1955. [25] G. Chinta and P. E. Gunnells. Constructing Weyl group multiple Dirichlet series. J. Amer. Math. Soc., 23(1):189–215, 2010. [26] G. Chinta, P. E. Gunnells, and A. Puskás.Metaplectic Demazure operators and Whittaker functions. Indiana Univ. Math. J., 66(3):1045–1064, 2017. [27] G. Chinta and O. Offen. A metaplectic Casselman-Shalika formula for GLr. Amer. J. Math., 135(2):403– 441, 2013. [28] I. Corwin and L. Petrov. Stochastic higher spin vertex models on the line. Comm. Math. Phys., 343(2):651– 700, 2016. [29] K. Costello. Integrable lattice models from four-dimensional field theories. In String-Math 2013, volume 88 of Proc. Sympos. Pure Math., pages 3–23. Amer. Math. Soc., Providence, RI, 2014. [30] P. Di Francesco, P. Mathieu, and D. Sénéchal. Conformal field theory. Graduate Texts in Contemporary Physics. Springer-Verlag, New York, 1997. [31] V. G. Drinfeld. Quasi-Hopf algebras. Algebra i Analiz, 1(6):114–148, 1989. [32] A. M. Hamel and R. C. King. Bijective proofs of shifted tableau and alternating sign matrix identities. J. Algebraic Combin., 25(4):417–458, 2007. [33] N. Iwahori and H. Matsumoto. On some Bruhat decomposition and the structure of the Hecke rings of p-adic Chevalley groups. Inst. Hautes Etudes´ Sci. Publ. Math., (25):5–48, 1965. [34] D. A. Kazhdan and S. J. Patterson. Metaplectic forms. Inst. Hautes Etudes´ Sci. Publ. Math., (59):35–142, 1984. [35] T. Kojima. Diagonalization of transfer matrix of supersymmetry Uq(slb(M + 1|N + 1)) chain with a boundary. J. Math. Phys., 54(4):043507, 40, 2013. [36] P. P. Kulish, N. Y. Reshetikhin, and E. K. Sklyanin. Yang-Baxter equations and representation theory. I. Lett. Math. Phys., 5(5):393–403, 1981. [37] H. Matsumoto. Sur les sous-groupes arithmétiquesdes groupes semi-simples déployés. Ann. Sci. Ecole´ Norm. Sup. (4), 2:1–62, 1969. [38] P. J. McNamara. Principal series representations of metaplectic groups over local fields. In Multiple Dirichlet series, L-functions and automorphic forms, volume 300 of Progr. Math., pages 299–327. Birkhäuser/Springer,New York, 2012. [39] P. J. McNamara. The metaplectic Casselman-Shalika formula. Trans. Amer. Math. Soc., 368(4):2913–2937, 2016. [40] C. Moeglin and J.-L. Waldspurger. Spectral decomposition and Eisenstein series, volume 113 of Cambridge Tracts in Mathematics. Cambridge University Press, Cambridge, 1995. Une paraphrase de l’Ecriture´ [A paraphrase of Scripture]. [41] M. Patnaik and A. Puskás. On Iwahori-Whittaker functions for metaplectic groups. Adv. Math., 313:875– 914, 2017. [42] M. M. Patnaik and A. Puskás. Metaplectic covers of Kac-Moody groups and Whittaker functions. Duke Math. J., 168(4):553–653, 2019. [43] J. H. H. Perk and C. L. Schultz. New families of commuting transfer matrices in q-state vertex models. Phys. Lett. A, 84(8):407–410, 1981. [44] A. Puskás.Whittaker functions on metaplectic covers of GL(r), 2016, arXiv:1605.05400. [45] N. Reshetikhin. Multiparameter quantum groups and twisted quasitriangular Hopf algebras. Lett. Math. Phys., 20(4):331–335, 1990. [46] S. Sahi, J. V. Stokman, and V. Venkateswaran. Metaplectic representations of Hecke algebras, Weyl group actions, and associated polynomials, 2018, arXiv:1808.01069. [47] G. Savin. Local Shimura correspondence. Math. Ann., 280(2):185–190, 1988.

44 [48] T. A. Springer. Reductive groups. In Automorphic forms, representations and L-functions (Proc. Sympos. Pure Math., Oregon State Univ., Corvallis, Ore., 1977), Part 1, Proc. Sympos. Pure Math., XXXIII, pages 3–27. Amer. Math. Soc., Providence, R.I., 1979. [49] T. Tokuyama. A generating function of strict Gelfand patterns and some formulas on characters of general linear groups. J. Math. Soc. Japan, 40(4):671–685, 1988. [50] M. H. Weissman. Split metaplectic groups and their L-groups. J. Reine Angew. Math., 696:89–141, 2014. [51] M. Wheeler and P. Zinn-Justin. Littlewood-Richardson coeﬃcients for Grothendieck polynomials from integrability. J. Reine Angew. Math., 757:159–195, 2019. [52] H. Zhang. Fundamental representations of quantum aﬃne superalgebras and R-matrices. Transform. Groups, 22(2):559–590, 2017. [53] R. B. Zhang. Structure and representations of the quantum general linear supergroup. Comm. Math. Phys., 195(3):525–547, 1998.

School of Mathematics, University of Minnesota, Minneapolis, MN 55455 Email address: [email protected]

Department of Mathematical and Statistical Sciences, University of Alberta, Edmonton, AB T6G 2G1, Canada Email address: [email protected]

Department of Mathematics, Stanford University, Stanford, CA 94305-2125 Email address: [email protected]

School of Mathematics, Institute for Advanced Study, Princeton, NJ 08540 Department of Mathematics, Rutgers University, Piscataway, NJ 08854 Department of Mathematical Sciences, University of Gothenburg and Chalmers University of Technology, SE-412 96 Gothenburg, Sweden Email address: [email protected]