Completeness of Bethe Ansatz for Gaudin Models with 𝔤𝔩(1|1) Symmetry and Diagonal Twists

Abstract

We studied the Gaudin models with $\mathfrak{gl}\left(1\right|1)$ symmetry that are twisted by a diagonal matrix and defined on tensor products of polynomial evaluation $\mathfrak{gl}\left(1\right|1\left)\right[t]$-modules. Namely, we gave an explicit description of the algebra of Hamiltonians (Gaudin Hamiltonians) acting on tensor products of polynomial evaluation $\mathfrak{gl}\left(1\right|1\left)\right[t]$-modules and showed that a bijection exists between common eigenvectors (up to proportionality) of the algebra of Hamiltonians and monic divisors of an explicit polynomial written in terms of the highest weights and evaluation parameters. In particular, our result implies that each common eigenspace of the algebra of Hamiltonians has dimension one. We also gave dimensions of the generalized eigenspaces.

1. Introduction

In the last half of a century, Gaudin models for simple Lie algebras have been intensively studied by many mathematicians and physicists using various methods, producing numerous spectacular results. For example, the simplicity of the spectrum of Gaudin algebra (Bethe algebra) was used to solve two long-standing conjectures: the transversality conjecture of the intersection of Schubert varieties and the Shapiro–Shapiro conjecture in real algebraic geometry; see [1]. Another example is that the monodromy of the joint eigenvectors of Gaudin algebra was proved to be given by the internal cactus group action on $\mathfrak{g}$-crystals, where $\mathfrak{g}$ is the corresponding finite-dimensional simple Lie algebra; see [2].

In recent years, the Gaudin models for Lie superalgebras have steadily gained attention within the mathematical community. For instance, the algebraic Bethe ansatz for Gaudin models of $\mathfrak{osp}\left(1\right|2)$ symmetry was carried out in [3]. Higher Gaudin Hamiltonians for Gaudin models of $\mathfrak{gl}\left(m\right|n)$ symmetry were constructed in [4] via studying the MacMahon Master Theorem related to Manin matrices. The completeness of Bethe ansatz for Gaudin models of $\mathfrak{gl}\left(m\right|n)$ symmetry that are defined on tensor products of vector representations was proved for the case of generic evaluation parameters in [5]. The relation between ${\mathrm{SPL}}_2$-superopers and the Bethe ansatz equations of $\mathfrak{osp}\left(1\right|2)$ Gaudin model was discussed in [6]. The reproduction procedure for Bethe ansatz equations of $\mathfrak{gl}\left(m\right|n)$ Gaudin models was introduced in [7]. Moreover, it was shown in [7] that the reproduction procedure gives rise to a variety that is isomorphic to the superflag variety. The duality between the quasi-periodic Gaudin model associated with Lie superalgebra $\mathfrak{gl}\left(m\right|n)$ and the quasi-periodic Gaudin model associated with Lie algebra $\mathfrak{gl}\left(k\right)$ was established in [8]. The reproduction procedure for Bethe ansatz equations of Gaudin models associated with orthosymplectic Lie superalgebras was introduced in [9]. In particular, this research developed the missing part of the reproduction procedure when the corresponding Lie algebras are of type D. Finally, in a previous work [10], we gave a complete answer of the periodic (twistless) Gaudin models of $\mathfrak{gl}\left(1\right|1)$ symmetry when the underlying Hilbert space is an arbitrary irreducible tensor product of evaluation polynomial modules. In this paper, we obtained the analogues for quasi-periodic $\mathfrak{gl}\left(1\right|1)$ Gaudin models; namely, we proved the completeness of Bethe ansatz for $\mathfrak{gl}\left(1\right|1)$ Gaudin models with diagonal twists.

The results of this paper are quite similar to those of [10,11], with suitable modifications, following the strategy of [1,12]. Surprisingly, to the best of our knowledge, most of the previous work on Gaudin models for Lie superalgebras was carried out in the periodic case, except, e.g., [8]. Therefore, we also need to establish the results on the algebraic Bethe ansatz for $\mathfrak{gl}\left(1\right|1)$ Gaudin models in the quasi-periodic case; see Section 2.4. In particular, we showed that the Bethe ansatz is complete for generic evaluation parameters; see Theorem 2. Using the completeness of the Bethe ansatz for generic parameters, we were able to describe the image of the algebra of Hamiltonians (Bethe algebra) explicitly and show that the quasi-periodic $\mathfrak{gl}\left(1\right|1)$ Gaudin models are perfectly integrable, cf. [13]. Consequently, we obtained the completeness of the Bethe ansatz for quasi-periodic $\mathfrak{gl}\left(1\right|1)$ Gaudin models with pairwise distinct evaluation parameters.

Note that the perfect integrability for the quasi-periodic $\mathfrak{gl}\left(m\right|n)$ Gaudin models defined on tensor products of symmetric powers of the vector representations was established in [8] [Corollary 5.3] by studying the duality between $\mathfrak{gl}\left(m\right|n)$ and $\mathfrak{gl}\left(k\right)$ Gaudin models and using the known results from [12]. In particular, it gives rise to the perfect integrability for the quasi-periodic $\mathfrak{gl}\left(1\right|1)$ Gaudin models defined on tensor products of polynomial modules. However, an explicit description of the image of Bethe algebra and the complete spectrum of Bethe algebra were not discussed in [8].

The paper is organized as follows. In Section 2, we fix notations and discuss basic facts of the algebraic Bethe ansatz for quasi-periodic $\mathfrak{gl}\left(1\right|1)$ Gaudin models. Then, we recall the space ${\mathcal{V}}^{\mathfrak{S}}$ and Weyl modules and their properties in Section 3. Section 4 contains the main theorems, where we also discuss the higher Gaudin transfer matrices and the relations between higher Gaudin transfer matrices and the first two Gaudin transfer matrices. Section 5 is dedicated to the proofs of main theorems.

2. Preliminaries

2.1. Lie Superalgebra $\mathfrak{gl}\left(1\right|1)$ and Its Representations

A vector superspace $V=V_{\overline{0}}\oplus V_{\overline{1}}$ is a ${\mathbb{Z}}_2$-graded vector space. Elements of $V_{\overline{0}}$ are called even; elements of $V_{\overline{1}}$ are called odd. We write $\left|v\right|\in \{\overline{0},\overline{1}\}$ for the parity of a homogeneous element $v\in V$. Set ${(-1)}^{\overline{0}}=1$ and ${(-1)}^{\overline{1}}=-1$.

Consider the vector superspace ${\mathbb{C}}^{1|1}$, where $dim\left({\mathbb{C}}_{\overline{0}}^{1|1}\right)=1$ and $dim\left({\mathbb{C}}_{\overline{1}}^{1|1}\right)=1$. We chose a homogeneous basis $v_1,v_2$ of ${\mathbb{C}}^{1|1}$ such that $|v_1|=\overline{0}$ and $|v_2|=\overline{1}$. For brevity, we shall write their parities as $|v_i|=|i|$. Denote by $E_{ij}\in \mathrm{End}\left({\mathbb{C}}^{1|1}\right)$ the linear operator of parity $\left|i\right|+\left|j\right|$ such that $E_{ij}{v}_r={\delta }_{jr}{v}_i$ for $i,j,r=1,2$.

The Lie superalgebra $\mathfrak{gl}\left(1\right|1)$ is spanned by elements $e_{ij}$, $i,j=1,2$, with parities $e_{ij}=\left|i\right|+\left|j\right|$, and the supercommutator relations are given by

$$[e_{ij},e_{rs}]={\delta }_{jr}{e}_{is}-{(-1)}^{\left(\right|i|+|j\left|\right)\left(\right|r|+|s\left|\right)}{\delta }_{is}{e}_{rj}.$$

Let $\mathfrak{h}$ be the commutative Lie subalgebra of 𝔤𝔩(1|1) spanned by e₁₁, e₂₂ Denote the universal enveloping algebras of 𝔤𝔩_1|1 and $\mathfrak{h}$ by $\mathrm{U}\left({\mathfrak{gl}}_{1|1}\right)$ and $\mathrm{U}\left(\mathfrak{h}\right)$, respectively.

We call a pair $\mathit{\lambda }=({\lambda }_1,{\lambda }_2)$ of complex numbers a $\mathfrak{gl}\left(1\right|1)$-weight. Set $\left|\mathit{\lambda }\right|={\lambda }_1+{\lambda }_2$. A $\mathfrak{gl}\left(1\right|1)$-weight $\mathit{\lambda }$ is non-degenerate if ${\lambda }_1+{\lambda }_2\ne 0$.

Let M be a $\mathfrak{gl}\left(1\right|1)$-module. A non-zero vector $v\in M$ is called singular if $e_{12}v=0$. Denote the subspace of all singular vectors of M by ${\left(M\right)}^{\mathrm{sing}}$. A non-zero vector $v\in M$ is called of weight $\mathit{\lambda }=({\lambda }_1,{\lambda }_2)$ if $e_{11}v={\lambda }_{1}v$ and $e_{22}v={\lambda }_{2}v$. Denote by ${\left(M\right)}_{\mathit{\lambda }}$ the subspace of M spanned by vectors of weight $\mathit{\lambda }$.

Let $\mathbf{\Lambda }=({\mathit{\lambda }}^{\left(1\right)},\cdots ,{\mathit{\lambda }}^{\left(k\right)})$ be a sequence of $\mathfrak{gl}\left(1\right|1)$-weights. Set $|\mathbf{\Lambda }|={\sum }_{s=1}^k\left|{\mathit{\lambda }}^{\left(s\right)}\right|$.

Denote by $L_{\mathit{\lambda }}$ the irreducible $\mathfrak{gl}\left(1\right|1)$-module generated by an even singular vector $v_{\mathit{\lambda }}$ of weight $\mathit{\lambda }$. Then, $L_{\mathit{\lambda }}$ is two-dimensional if $\mathit{\lambda }$ is non-degenerate and one-dimensional otherwise. Clearly, ${\mathbb{C}}^{1|1}\cong L_{{\omega }_1}$, where ${\omega }_1=(1,0)$, if we identify the action of $e_{ij}$ on ${\mathbb{C}}^{1|1}$ with the operator $E_{ij}$.

A $\mathfrak{gl}\left(1\right|1)$-module M is called a polynomial module if M is a submodule of ${\left({\mathbb{C}}^{1|1}\right)}^{\otimes n}$ for some $n\in {\mathbb{Z}}_{⩾0}$. We say that $\mathit{\lambda }$ is a polynomial weight if $L_{\mathit{\lambda }}$ is a polynomial module. Weight $\mathit{\lambda }=({\lambda }_1,{\lambda }_2)$ is a polynomial weight if and only if ${\lambda }_1,{\lambda }_2\in {\mathbb{Z}}_{⩾0}$ and either ${\lambda }_1>0$ or ${\lambda }_1={\lambda }_2=0$. We also write $L_{({\lambda }_1,{\lambda }_2)}$ for $L_{\mathit{\lambda }}$.

For non-degenerate polynomial weights $\mathit{\lambda }=({\lambda }_1,{\lambda }_2)$ and $\mathit{\mu }=({\mu }_1,{\mu }_2)$, we have

$$L_{({\lambda }_1,{\lambda }_2)}\otimes L_{({\mu }_1,{\mu }_2)}=L_{({\lambda }_1+{\mu }_1,{\lambda }_2+{\mu }_2)}\oplus L_{({\lambda }_1+{\mu }_1-1,{\lambda }_2+{\mu }_2+1)}.$$

2.2. Current Superalgebra $\mathfrak{gl}\left(1\right|1\left)\right[t]$

Denote by $\mathfrak{gl}\left(1\right|1\left)\right[t]$ the Lie superalgebra $\mathfrak{gl}\left(1\right|1)\otimes \mathbb{C}[t]$ of ${\mathfrak{gl}}_{1|1}$-valued polynomials with the point-wise supercommutator. Call $\mathfrak{gl}\left(1\right|1\left)\right[t]$ the current superalgebra of $\mathfrak{gl}\left(1\right|1)$. We identify $\mathfrak{gl}\left(1\right|1)$ with the subalgebra $\mathfrak{gl}\left(1\right|1)\otimes 1$ of constant polynomials in $\mathfrak{gl}\left(1\right|1\left)\right[t]$.

We write $e_{ij}\left[r\right]$ for $e_{ij}\otimes t^r$, $r\in {\mathbb{Z}}_{⩾0}$. A basis of $\mathfrak{gl}\left(1\right|1\left)\right[t]$ is given by $e_{ij}\left[r\right]$, $i,j=1,2$ and $r\in {\mathbb{Z}}_{⩾0}$. They satisfy the supercommutator relations

$$[e_{ij}\left[r\right],e_{kl}\left[s\right]]={\delta }_{jk}{e}_{il}[r+s]-{(-1)}^{\left(\right|i|+|j\left|\right)\left(\right|k|+|l\left|\right)}{\delta }_{il}{e}_{kj}[r+s].$$

In particular, one has

$${\left(e_{12}\left[r\right]\right)}^2={\left(e_{21}\left[r\right]\right)}^2=0,e_{21}\left[r\right]e_{21}\left[s\right]=-e_{21}\left[s\right]e_{21}\left[r\right]$$

(1)

in the universal enveloping superalgebra $\mathrm{U}\left(\mathfrak{gl}\right(1\left|1\right)\left[t\right])$. The universal enveloping superalgebra $\mathrm{U}\left(\mathfrak{gl}\right(1\left|1\right)\left[t\right])$ is a Hopf superalgebra with the coproduct given by

$$\Delta \left(X\right)=X\otimes 1+1\otimes X,\mathrm{for}X\in \mathfrak{gl}\left(1\right|1\left)\right[t].$$

Let $e_{ij}\left(x\right)={\sum }_{r=0}^{\infty }{e}_{ij}\left[r\right]x^{-r-1}$, where x is a formal variable. Then, we have

$$(u-v)[e_{ij}\left(u\right),e_{rs}\left(v\right)]=-[e_{ij},e_{rs}]\left(u\right)+[e_{ij},e_{rs}]\left(v\right).$$

(2)

In particular,

$$[e_{ij}\left(x\right),e_{rs}\left(x\right)]=-{\partial }_x[e_{ij},e_{rs}]\left(x\right).$$

(3)

For each $a\in \mathbb{C}$, there exists an automorphism of $\mathrm{U}\left(\mathfrak{gl}\right(1\left|1\right)\left[t\right])$, ${\rho }_a:e_{ij}\left(x\right)\to e_{ij}(x-a)$. Given a $\mathfrak{gl}\left(1\right|1\left)\right[t]$-module M, denote by $M\left(a\right)$ the pull-back of M through the automorphism ${\rho }_a$.

For each $a\in \mathbb{C}$, we have the evaluation map

$${\mathsf{ev}}_a:\mathrm{U}\left(\mathfrak{gl}\left(1\right|1)\left[t\right]\right)\to \mathrm{U}\left(\mathfrak{gl}\left(1\right|1)\right),e_{ij}\left(x\right)\mapsto e_{ij}/(x-a).$$

For a $\mathfrak{gl}\left(1\right|1)$-module L, denote by $L\left(a\right)$ the $\mathfrak{gl}\left(1\right|1\left)\right[t]$-module obtained by pulling back L through the evaluation map ${\mathsf{ev}}_a$. We call $L\left(a\right)$ an evaluation module at a.

Given any series $\zeta \left(x\right)\in x^{-1}\mathbb{C}\left[x^{-1}\right]$, we have the one-dimensional $\mathfrak{gl}\left(1\right|1\left)\right[t]$-module generated by an even vector v satisfying $e_{ij}\left(x\right)v={\delta }_{ij}{(-1)}^{\left|j\right|}\zeta \left(x\right)v$. We denote this module by ${\mathbb{C}}_{\zeta }$.

If $b_1,\cdots ,b_n$ are pairwise distinct complex numbers and $L_1,\dots ,L_n$ are finite-dimensional irreducible $\mathfrak{gl}\left(1\right|1)$-modules, then the $\mathfrak{gl}\left(1\right|1\left)\right[t]$-module ${⨂}_{s=1}^n{L}_s\left(b_s\right)$ is irreducible.

There is a natural ${\mathbb{Z}}_{⩾0}$-gradation on $\mathrm{U}\left(\mathfrak{gl}\right(1\left|1\right)\left[t\right])$ such that $deg\left(e_{ij}\left[r\right]\right)=r$ which induces the filtration ${\mathcal{F}}_0\mathrm{U}\left(\mathfrak{gl}\left(1\right|1)\left[t\right]\right)\subset {\mathcal{F}}_1\mathrm{U}\left(\mathfrak{gl}\left(1\right|1)\left[t\right]\right)\subset \cdots \subset \mathrm{U}\left(\mathfrak{gl}\left(1\right|1)\left[t\right]\right)$, where ${\mathcal{F}}_s\mathrm{U}\left(\mathfrak{gl}\left(1\right|1)\left[t\right]\right)$ is the subspace of $\mathrm{U}\left(\mathfrak{gl}\right(1\left|1\right)\left[t\right])$ spanned by all elements of degree $⩽s$.

Let M be a ${\mathbb{Z}}_{⩾0}$-graded space with finite-dimensional homogeneous components. Let $M_j\subset M$ be the homogeneous component of degree j. We call the formal power series in variable q,

$$\mathrm{ch}\left(M\right)=\sum _{j=0}^{\infty }dim\left(M_j\right)q^j,$$

(4)

the graded character of M.

2.3. Gaudin Hamiltonians

In this section, we discuss the inhomogeneous Gaudin Hamiltonians. Throughout the paper, we shall fix two complex numbers $\mathit{q}=(q_1,q_2)$. Moreover, we assume that $q_1\ne q_2$; see the end of this section.

Let $\mathit{b}=(b_1,\cdots ,b_k)$ be a sequence of distinct complex numbers and $\mathbf{\Lambda }=({\mathit{\lambda }}^{\left(1\right)},\dots ,{\mathit{\lambda }}^{\left(k\right)})$ a sequence of polynomial $\mathfrak{gl}\left(1\right|1)$-weights, where ${\mathit{\lambda }}^{\left(s\right)}=({\alpha }_s,{\beta }_s)$.

Set $n=|\mathbf{\Lambda }|={\sum }_{s=1}^k({\alpha }_s+{\beta }_s)$ and $L_{\mathbf{\Lambda }}={⨂}_{s=1}^k{L}_{{\mathit{\lambda }}^{\left(s\right)}}$. The quadratic Gaudin Hamiltonians are the linear maps ${\mathcal{H}}_r\in \mathrm{End}\left(L_{\mathbf{\Lambda }}\right)$ given by

$${\mathcal{H}}_r:=q_1{e}_{11}^{\left(r\right)}+q_2{e}_{22}^{\left(r\right)}+\sum _{s=1,s\ne r}^k\frac{e_{11}^{\left(r\right)}{e}_{11}^{\left(s\right)}-e_{12}^{\left(r\right)}{e}_{21}^{\left(s\right)}+e_{21}^{\left(r\right)}{e}_{12}^{\left(s\right)}-e_{22}^{\left(r\right)}{e}_{22}^{\left(s\right)}}{b_r-b_s},1⩽r⩽k.$$

(5)

where $e_{ab}^{\left(r\right)}=1^{\otimes (r-1)}\otimes e_{ab}\otimes 1^{\otimes (k-r)}$.

Lemma 1.

The Gaudin Hamiltonians${\mathcal{H}}_r$

• Are mutually commuting:$[{\mathcal{H}}_r,{\mathcal{H}}_s]=0$for all$r,s$;
• Commute with the action of 𝔥: $[{\mathcal{H}}_r,X]=0$ for all r and $X\in \mathfrak{h}$.

Proof. 

This follows immediately from [5] [Proposition 3.1] for non-twisted (i.e., $q_1=q_2=0$) Gaudin Hamiltonians. □

Instead of working on Gaudin Hamiltonians ${\mathcal{H}}_s$, we work on the generating function of Gaudin Hamiltonians,

$$\mathcal{H}\left(x\right):=\sum _{r=1}^{\infty }{\mathcal{H}}_r{x}^{-r}=q_1{e}_{11}\left(x\right)+q_2{e}_{22}\left(x\right)+\frac{1}{2}\sum _{a,b=1}^2{e}_{ab}\left(x\right)e_{ba}\left(x\right){(-1)}^{\left|b\right|}.$$

(6)

The operator $\mathcal{H}\left(x\right)$ acts on the tensor product of the evaluation $\mathfrak{gl}\left(1\right|1\left)\right[t]$-modules

$$L_{\mathbf{\Lambda }}\left(\mathit{b}\right):=\underset{s=1}{\overset{k}{⨂}}{L}_{{\mathit{\lambda }}^{\left(s\right)}}\left(b_s\right).$$

Note that $L_{\mathbf{\Lambda }}\left(\mathit{b}\right)$ and $L_{\mathbf{\Lambda }}$ are isomorphic as $\mathfrak{gl}\left(1\right|1)$-modules via the identity map; then, we have

$$\mathcal{H}\left(x\right)=\frac{1}{2}\sum _{s=1}^k\frac{{\alpha }_s({\alpha }_s-1)-{\beta }_s({\beta }_s+1)}{{(x-b_s)}^2}\mathrm{Id}+\sum _{s=1}^k\frac{1}{x-b_s}{\mathcal{H}}_s,$$

(7)

as operators in $\mathrm{End}\left(L_{\mathbf{\Lambda }}\right)=\mathrm{End}\left(L_{\mathbf{\Lambda }}\left(\mathit{b}\right)\right)$. We call $\mathcal{H}\left(x\right)$ the Gaudin transfer matrix.

We are interested in finding the eigenvalues and eigenvectors of the Gaudin transfer matrix in $L_{\mathbf{\Lambda }}\left(\mathit{b}\right)$. To be more precise, we call

$$\xi \left(x\right)=\sum _{r=1}^{\infty }{\xi }_r{x}^{-r},{\xi }_r\in \mathbb{C},$$

(8)

an eigenvalue of $\mathcal{H}\left(x\right)$ if there exists a non-zero vector $v\in L_{\mathbf{\Lambda }}\left(\mathit{b}\right)$ such that ${\mathcal{H}}_{r}v={\xi }_{r}v$ for all $r\in {\mathbb{Z}}_{⩾1}$. If $\xi \left(x\right)$ is a rational function, we consider it as a power series in $x^{-1}$ as in (8). The vector v is called an eigenvector of $\mathcal{H}\left(x\right)$ corresponding to eigenvalue $\xi \left(x\right)$. We also define the eigenspace of $\mathcal{H}\left(x\right)$ in $L_{\mathbf{\Lambda }}\left(\mathit{b}\right)$ corresponding to eigenvalue $\xi \left(x\right)$ as ${\bigcap }_{r=1}^{\infty }ker({\mathcal{H}}_r{|}_{L_{\mathbf{\Lambda }}\left(\mathit{b}\right)}-{\xi }_r)$.

It is sufficient to consider $L_{\mathbf{\Lambda }}$ with ${\beta }_s=0$ for all s. Indeed, if $L_{\mathbf{\Lambda }}\left(\mathit{b}\right)$ is an arbitrary tensor product and

$$\xi \left(x\right)=\sum _{s=1}^k\frac{{\beta }_s}{x-b_s},$$

then

$$L_{\mathbf{\Lambda }}\left(\mathit{b}\right)\otimes {\mathbb{C}}_{\xi }\cong L_{\tilde{\mathbf{\Lambda }}}\left(\mathit{b}\right),{\tilde{\mathit{\lambda }}}^{\left(s\right)}=({\alpha }_s+{\beta }_s,0).$$

Identify $L_{\mathbf{\Lambda }}\left(\mathit{b}\right)\otimes {\mathbb{C}}_{\xi }$ with $L_{\mathbf{\Lambda }}\left(\mathit{b}\right)$ as vector spaces. Then, $\mathcal{H}\left(x\right)$ acting on $L_{\mathbf{\Lambda }}\left(\mathit{b}\right)\otimes {\mathbb{C}}_{\xi }$ coincides with $\mathcal{H}\left(x\right)+\zeta \left(x\right)(e_{11}\left(x\right)+e_{22}\left(x\right))+(q_1-q_2)\xi \left(x\right)$ acting on $L_{\mathbf{\Lambda }}\left(\mathit{b}\right)$. Note that the coefficients of $e_{11}\left(x\right)+e_{22}\left(x\right)$ are central in $\mathrm{U}\left(\mathfrak{gl}\right(1\left|1\right)\left[t\right])$ and hence $e_{11}\left(x\right)+e_{22}\left(x\right)$ acts on $L_{\mathbf{\Lambda }}\left(\mathit{b}\right)$ by the scalar series

$$\sum _{s=1}^k\frac{{\alpha }_s+{\beta }_s}{x-b_s};$$

therefore, the problem of the diagonalization of the Gaudin transfer matrix in $L_{\mathbf{\Lambda }}\left(\mathit{b}\right)$ is reduced to diagonalization of the Gaudin transfer matrix in $L_{\tilde{\mathbf{\Lambda }}}\left(\mathit{b}\right)$.

Again, by the fact that the coefficients of $e_{11}\left(x\right)+e_{22}\left(x\right)$ are central, if $q_1=q_2$, then the diagonalization problem of $\mathcal{H}\left(x\right)$ is the same as the one for the homogeneous case $q_1=q_2=0$, which was discussed in [10]. Thus, for the rest of the paper, we shall assume that $q_1\ne q_2$.

Since $L_{\mathit{\lambda }}$ is one-dimensional if $\mathit{\lambda }$ is degenerate, similarly, it suffices to consider the case that all participant $\mathfrak{gl}\left(1\right|1)$-weights are non-degenerate. Hence, we shall always assume throughout the paper that ${\mathit{\lambda }}^{\left(s\right)}$ are non-degenerate for all $1⩽s⩽k$.

2.4. Bethe Ansatz

The main method to find eigenvalues and eigenvectors of the Gaudin transfer matrix in $L_{\mathbf{\Lambda }}$ is the algebraic Bethe ansatz. We give the results for the algebraic Bethe ansatz of quasi-periodic $\mathfrak{gl}\left(1\right|1)$ Gaudin models in this section following e.g., [5] [Section VI].

Fix a non-negative integer l. Let $\mathit{t}=(t_1,\dots ,t_l)$ be a sequence of complex numbers. Define the polynomial $y_{\mathit{t}}={\prod }_{i=1}^l(x-t_i)$. We say that polynomial $y_{\mathit{t}}$ represents $\mathit{t}$.

Set

$${\zeta }_{\mathbf{\Lambda },\mathit{b}}\left(x\right):=q_1-q_2+\sum _{s=1}^k\frac{{\alpha }_s+{\beta }_s}{x-b_s}.$$

(9)

A sequence of complex numbers $\mathit{t}$ is called a solution to the Bethe ansatz equation associated to $\mathbf{\Lambda }$, $\mathit{b}$, l if

$$y_{\mathit{t}}\left(x\right)\mathrm{divides}\mathrm{the}\mathrm{polynomial}{\phi }_{\mathbf{\Lambda },\mathit{b}}\left(x\right):={\zeta }_{\mathbf{\Lambda },\mathit{b}}\left(x\right)\prod _{s=1}^k(x-b_s).$$

(10)

We do not distinguish solutions that differ by a permutation of coordinates (that is represented by the same polynomial).

Let $v_s$ be the highest weight vector of $L_{{\mathit{\lambda }}^{\left(s\right)}}$, and set $|0\rangle =v_1\otimes \cdots \otimes v_k$. We call $|0\rangle$ the vacuum vector.

Define the off-shell Bethe vector ${\mathbb{B}}_l\left(\mathit{t}\right)\in {\left(L_{\mathbf{\Lambda }}\right)}_{(n-l,l)}$ by

$${\mathbb{B}}_l\left(\mathit{t}\right)=e_{21}\left(t_1\right)\cdots e_{21}\left(t_l\right)|0\rangle .$$

(11)

Since $e_{21}\left(x\right)e_{21}\left(u\right)=-e_{21}\left(u\right)e_{21}\left(x\right)$, the order of $t_i$ is not important. Moreover, the off-shell Bethe vector is zero if $t_i=t_j$ for some $1⩽i\ne j⩽l$.

If $\mathit{t}$ is a solution of the Bethe ansatz Equation (10), we call ${\mathbb{B}}_l\left(\mathit{t}\right)$ an on-shell Bethe vector.

Let $\mathit{t}$ be a solution of the Bethe ansatz equation associated to $\mathbf{\Lambda }$, $\mathit{b}$, l.

Theorem 1.

If the on-shell Bethe vector${\mathbb{B}}_l\left(\mathit{t}\right)$is non-zero, then${\mathbb{B}}_l\left(\mathit{t}\right)$is an eigenvector of the Gaudin transfer matrix$\mathcal{H}\left(x\right)$with the corresponding eigenvalue

$${\mathcal{E}}_{y_{\mathit{t}},\mathbf{\Lambda },\mathit{b}}\left(x\right)=\frac{1}{2}{\zeta }_{\mathbf{\Lambda },\mathit{b}}^{\prime }\left(x\right)-{\zeta }_{\mathbf{\Lambda },\mathit{b}}\left(x\right)\frac{y_{\mathit{t}}^{\prime }\left(x\right)}{y_{\mathit{t}}\left(x\right)}+\sum _{r,s=1}^k\frac{{\alpha }_r{\alpha }_s-{\beta }_r{\beta }_s}{2(x-b_r)(x-b_s)}+\sum _{s=1}^k\frac{q_1{\alpha }_s+q_2{\beta }_s}{x-b_s}.$$

(12)

where${\zeta }_{\mathbf{\Lambda },\mathit{b}}\left(x\right)$is given by (9).

Proof. 

By (2) and the fact that coefficients of $e_{11}\left(x\right)+e_{22}\left(x\right)$ are central in $\mathrm{U}\left(\mathfrak{gl}\right(1\left|1\right)\left[t\right])$, we have

$$[\mathcal{H}\left(x\right),e_{21}\left(t\right)]=-\frac{1}{u-t}\left({\zeta }_{\mathbf{\Lambda },\mathit{b}}\left(x\right)e_{21}\left(t\right)-{\zeta }_{\mathbf{\Lambda },\mathit{b}}\left(t\right)e_{21}\left(x\right)\right),$$

as operators on $L_{\mathbf{\Lambda }}\left(\mathit{b}\right)$. Note that if t is a coordinate of a solution of the Bethe ansatz equation, then ${\zeta }_{\mathbf{\Lambda },\mathit{b}}\left(t\right)=0$. Therefore, we have

$$[\mathcal{H}\left(x\right),e_{21}\left(t_i\right)]=-\frac{1}{u-t}{\zeta }_{\mathbf{\Lambda },\mathit{b}}\left(x\right)e_{21}\left(t_i\right)$$

for $1⩽i⩽l$. Hence, we conclude that

$$\begin{array}{c}\hfill \mathcal{H}\left(x\right){\mathbb{B}}_l\left(\mathit{t}\right)=-{\zeta }_{\mathbf{\Lambda },\mathit{b}}\left(x\right)\sum _{j=1}^l\frac{1}{x-t_j}{\mathbb{B}}_l\left(\mathit{t}\right)+e_{21}\left(t_1\right)\cdots e_{21}\left(t_l\right)\mathcal{H}\left(x\right)|0\rangle .\end{array}$$

The theorem now follows from the straightforward computation of the eigenvalue of $\mathcal{H}\left(x\right)$ corresponding to the vector $|0\rangle$. □

Consider another Gaudin transfer matrix

$$\mathcal{T}\left(x\right)=\frac{1}{2}\left({\dot{e}}_{11}\left(x\right)+{\dot{e}}_{22}\left(x\right)\right)+\frac{1}{2}{\left(e_{11}\left(x\right)+e_{22}\left(x\right)\right)}^2+q_1\left(e_{11}\left(x\right)+e_{22}\left(x\right)\right)-\mathcal{H}\left(x\right),$$

(13)

where ${\dot{e}}_{ii}\left(x\right)={\partial }_x\left(e_{ii}\left(x\right)\right)$, $i=1,2$. Then, the eigenvalue of $\mathcal{T}\left(x\right)$ acting on the on-shell Bethe vector ${\mathbb{B}}_l\left(\mathit{t}\right)$ is

$$\begin{array}{cc}\hfill {\mathcal{E}}_{y_{\mathit{t}},\mathbf{\Lambda },\mathit{b}}\left(x\right)=& {\zeta }_{\mathbf{\Lambda },\mathit{b}}\left(x\right)\frac{y_{\mathit{t}}^{\prime }\left(x\right)}{y_{\mathit{t}}\left(x\right)}+\sum _{r,s=1}^k\frac{{\alpha }_r{\beta }_s+{\alpha }_s{\beta }_r+2{\beta }_r{\beta }_s}{2(x-b_r)(x-b_s)}+\sum _{s=1}^k\frac{(q_1-q_2){\beta }_s}{x-b_s}\hfill \\ \hfill =& {\zeta }_{\mathbf{\Lambda },\mathit{b}}\left(x\right)\left(\frac{y_{\mathit{t}}^{\prime }\left(x\right)}{y_{\mathit{t}}\left(x\right)}+\sum _{s=1}^k\frac{{\beta }_s}{x-b_s}\right).\hfill \end{array}$$

(14)

It is important to know if the on-shell Bethe vectors are non-zero.

Proposition 1.

Suppose that the polynomial${\phi }_{\mathbf{\Lambda },\mathit{b}}\left(x\right)$only has simple roots; then, the on-shell Bethe vector${\mathbb{B}}_l\left(\mathit{t}\right)$is nonzero.

Proof. 

Since ${\phi }_{\mathbf{\Lambda },\mathit{b}}\left(x\right)$ only has simple roots, we have $t_i\ne t_j$ for $i\ne j$. Note that $b_s$ are distinct and ${\alpha }_s+{\beta }_s>0$ (since the weights are nondegenerate by our assumption); then, we have $b_s\ne t_i$. Hence, ${\zeta }_{\mathbf{\Lambda },\mathit{b}}\left(t_i\right)=0$. Moreover, we have

$$0\ne {\phi }_{\mathbf{\Lambda },\mathit{b}}^{\prime }\left(t_i\right)={\zeta }_{\mathbf{\Lambda },\mathit{b}}^{\prime }\left(t_i\right)\prod _{s=1}^k(t_i-b_s)+{\zeta }_{\mathbf{\Lambda },\mathit{b}}\left(t_i\right){\left(\prod _{s=1}^k(x-b_s))\right)}^{\prime }{|}_{x=t_i}={\zeta }_{\mathbf{\Lambda },\mathit{b}}^{\prime }\left(t_i\right)\prod _{s=1}^k(t_i-b_s).$$

Therefore, ${\zeta }_{\mathbf{\Lambda },\mathit{b}}^{\prime }\left(t_i\right)\ne 0$.

By (2) and the fact that coefficients of $e_{11}\left(x\right)+e_{22}\left(x\right)$ are central in $\mathrm{U}\left(\mathfrak{gl}\right(1\left|1\right)\left[t\right])$, we have

$$[e_{12}\left(t\right),e_{21}\left(\tilde{t}\right)]=-\frac{1}{t-\tilde{t}}\left({\zeta }_{\mathbf{\Lambda },\mathit{b}}\left(t\right)-{\zeta }_{\mathbf{\Lambda },\mathit{b}}\left(\tilde{t}\right)\right)$$

as operators on $L_{\mathbf{\Lambda }}\left(\mathit{b}\right)$. Therefore, we have $[e_{12}\left(t\right),e_{21}\left(\tilde{t}\right)]=0$ if t and $\tilde{t}$ are distinct coordinates of $\mathit{t}$ while $[e_{12}\left(t\right),e_{21}\left(t\right)]=-{\zeta }_{\mathbf{\Lambda },\mathit{b}}^{\prime }\left(t\right)$. One finds that

$$\begin{array}{c}\hfill e_{12}\left(t_l\right)\cdots e_{12}\left(t_1\right)e_{21}\left(t_1\right)\cdots e_{21}\left(t_l\right)|0\rangle ={(-1)}^l\prod _{i=1}^l{\zeta }_{\mathbf{\Lambda },\mathit{b}}^{\prime }\left(t_i\right)|0\rangle \ne 0,\end{array}$$

completing the proof. □

The conjecture of the completeness of the Bethe ansatz for Gaudin models associated with $\mathfrak{gl}\left(1\right|1)$ was formulated as follows, cf. [7] [Conjecture 8.3].

Conjecture 1.

Suppose all weights${\mathit{\lambda }}^{\left(s\right)}$,$1⩽s⩽k$are polynomial$\mathfrak{gl}\left(1\right|1)$-weights. Then, the Gaudin transfer matrix$\mathcal{H}\left(x\right)$has a simple spectrum in$L_{\mathbf{\Lambda }}\left(\mathit{b}\right)$. There exists a bijective correspondence between the monic divisors y of the polynomial${\phi }_{\mathbf{\Lambda },\mathit{b}}$and the eigenvectors v of the Gaudin transfer matrices (up to multiplication by a non-zero constant). Moreover, this bijection is such that$\mathcal{H}\left(x\right)v={\mathcal{E}}_{y,\mathbf{\Lambda },\mathit{b}}\left(x\right)v$, where${\mathcal{E}}_{y,\mathbf{\Lambda },\mathit{b}}\left(x\right)$is given by (12).

By simple spectrum, we mean that if $v_1$, $v_2$ are eigenvectors of $\mathcal{H}\left(x\right)$ and $v_1\ne c{v}_2$, $c\in {\mathbb{C}}^{\times }$, then the eigenvalues of $\mathcal{H}\left(x\right)$ on $v_1$ and $v_2$ are different.

The conjecture follows from Theorem 4 proved in Section 5.3.

The conjecture is clear for the case when ${\phi }_{\mathbf{\Lambda },\mathit{b}}$ only has simple roots. Note that $dim{L}_{\mathbf{\Lambda }}\left(\mathit{b}\right)=2^k$. If the polynomial ${\phi }_{\mathbf{\Lambda },\mathit{b}}$ has no multiple roots, then ${\phi }_{\mathbf{\Lambda },\mathit{b}}$ has the desired number of distinct monic divisors. Therefore, we have the desired number of on-shell Bethe vectors, which are also nonzero by Proposition 1. By Theorem 1, it implies that we do have an eigenbasis of the Gaudin transfer matrix consisting of on-shell Bethe vectors in $L_{\mathbf{\Lambda }}\left(\mathit{b}\right)$ with different eigenvalues. Thus, the algebraic Bethe ansatz works well for this situation.

Theorem 2.

Suppose that all weights${\mathit{\lambda }}^{\left(s\right)}$,$1⩽s⩽k$are polynomial$\mathfrak{gl}\left(1\right|1)$-weights. If the polynomial${\phi }_{\mathbf{\Lambda },\mathit{b}}$has no multiple roots, then the Gaudin transfer matrix$\mathcal{H}\left(x\right)$is diagonalizable and the Bethe ansatz is complete. In particular, for any givenΛand generic$\mathit{b}$, the Gaudin transfer matrix$\mathcal{H}\left(x\right)$is diagonalizable and the Bethe ansatz is complete.

3. Space ${\mathcal{V}}^{\mathfrak{S}}$ and Weyl Modules

In this section, we discuss the super-analog of ${\mathcal{V}}^S$ in [1] [Section 2.5], cf. [11] [Section 3].

The symmetric group ${\mathfrak{S}}_n$ acts naturally on $\mathbb{C}[z_1,\dots ,z_n]$ by permuting variables. Denote by ${\sigma }_i\left(\mathit{z}\right)$ the i-th elementary symmetric polynomial in $z_1,\dots ,z_n$. The algebra of symmetric polynomials $\mathbb{C}{[z_1,\dots ,z_n]}^{\mathfrak{S}}$ is freely generated by ${\sigma }_1\left(\mathit{z}\right),\dots ,{\sigma }_n\left(\mathit{z}\right)$.

Fix $\ell \in \{0,1,\dots ,n\}$. We have a subgroup ${\mathfrak{S}}_{\ell }\times {\mathfrak{S}}_{n-\ell }\subset {\mathfrak{S}}_n$. Then, ${\mathfrak{S}}_{\ell }$ permutes the first ℓ variables, whereas ${\mathfrak{S}}_{n-\ell }$ permutes the last $n-\ell$ variables. Denote by

$$\mathbb{C}{[z_1,\dots ,z_n]}^{{\mathfrak{S}}_{\ell }\times {\mathfrak{S}}_{n-\ell }}$$

the subalgebra of $\mathbb{C}[z_1,\dots ,z_n]$ consisting of ${\mathfrak{S}}_{\ell }\times {\mathfrak{S}}_{n-\ell }$-invariant polynomials. It is known that $\mathbb{C}{[z_1,\dots ,z_n]}^{{\mathfrak{S}}_{\ell }\times {\mathfrak{S}}_{n-\ell }}$ is a free $\mathbb{C}{[z_1,\dots ,z_n]}^{\mathfrak{S}}$-module of rank $\left(\genfrac{}{}{0pt}{}{n}{\ell }\right)$.

3.1. Definition of ${\mathcal{V}}^{\mathfrak{S}}$

Let $V={\left({\mathbb{C}}^{1|1}\right)}^{\otimes n}$ be the tensor power of the vector representation of $\mathfrak{gl}\left(1\right|1)$. The $\mathfrak{gl}\left(1\right|1)$-module V has weight decomposition

$$V=\underset{\ell =0}{\overset{n}{⨁}}{\left(V\right)}_{(n-\ell ,\ell )}.$$

Let $\mathcal{V}$ be the space of polynomials in variables $\mathit{z}=(z_1,z_2,\dots ,z_n)$ with coefficients in V,

$$\mathcal{V}=V\otimes \mathbb{C}[z_1,z_2,\dots ,z_n].$$

The space V is identified with the subspace $V\otimes 1$ of constant polynomials in $\mathcal{V}$. The space $\mathcal{V}$ has a natural grading induced from the grading on $\mathbb{C}[z_1,\dots ,z_n]$ with $deg\left(z_i\right)=1$. Namely, the degree of an element $v\otimes p$ in $\mathcal{V}$ is given by the degree of the polynomial p, $deg(v\otimes p)=degp$. Clearly, the space $\mathrm{End}\left(\mathcal{V}\right)$ has a gradation structure induced from that on $\mathcal{V}$.

Let $P^{(i,j)}$ be the graded flip operator that acts on the i-th and j-th factors of V. Let $s_1$, $s_2$, …, $s_{n-1}$ be the simple permutations of the symmetric group ${\mathfrak{S}}_n$. Define the ${\mathfrak{S}}_n$-action on $\mathcal{V}$ by the rule:

$$s_i:\mathit{f}(z_1,\dots ,z_n)\mapsto P^{(i,i+1)}\mathit{f}(z_1,\dots ,z_{i+1},z_i,\dots ,z_n),$$

for $\mathit{f}(z_1,\dots ,z_n)\in \mathcal{V}$. Note that the ${\mathfrak{S}}_n$-action respects the gradation on $\mathcal{V}$. Denote the subspace of all vectors in $\mathcal{V}$ invariant with respect to the ${\mathfrak{S}}_n$-action by ${\mathcal{V}}^{\mathfrak{S}}$.

Clearly, the $\mathfrak{gl}\left(1\right|1)$-action on $\mathcal{V}$ commutes with the ${\mathfrak{S}}_n$-action on $\mathcal{V}$ and preserves the grading. Therefore, ${\mathcal{V}}^{\mathfrak{S}}$ is a graded $\mathfrak{gl}\left(1\right|1)$-module. Hence, we have the weight decomposition for both ${\mathcal{V}}^{\mathfrak{S}}$ and ${\left({\mathcal{V}}^{\mathfrak{S}}\right)}^{\mathrm{sing}}$:

$${\mathcal{V}}^{\mathfrak{S}}=\underset{\ell =0}{\overset{n}{⨁}}{\left({\mathcal{V}}^{\mathfrak{S}}\right)}_{(n-\ell ,\ell )},{\left({\mathcal{V}}^{\mathfrak{S}}\right)}^{\mathrm{sing}}=\underset{\ell =0}{\overset{n}{⨁}}{\left({\mathcal{V}}^{\mathfrak{S}}\right)}_{(n-\ell ,\ell )}^{\mathrm{sing}}.$$

Note that ${\left({\mathcal{V}}^{\mathfrak{S}}\right)}_{(n-\ell ,\ell )}$ and ${\left({\mathcal{V}}^{\mathfrak{S}}\right)}_{(n-\ell ,\ell )}^{\mathrm{sing}}$ are also graded $\mathbb{C}{[z_1,\dots ,z_n]}^{\mathfrak{S}}$-modules.

The space $\mathcal{V}$ is a $\mathfrak{gl}\left(1\right|1\left)\right[t]$-module where $e_{ij}\left[r\right]$ acts by

$$\begin{array}{cc}\hfill e_{ij}\left[r\right]& (p(z_1,\dots ,z_n)w_1\otimes \cdots \otimes w_n)\hfill \\ \hfill =& p(z_1,\dots ,z_n)\sum _{s=1}^n{(-1)}^{\left(\right|w_1|+\cdots +|w_{s-1}\left|\right)\left(\right|i|+\left|j\right|)}{z}_s^r{w}_1\otimes \cdots \otimes e_{ij}{w}_s\otimes \cdots \otimes w_n,\hfill \end{array}$$

(15)

for $p(z_1,\dots ,z_n)\in \mathbb{C}[z_1,\dots ,z_n]$ and $w_s\in {\mathbb{C}}^{1|1}$.

Lemma 2.

The$\mathfrak{gl}\left(1\right|1\left)\right[t]$-action on$\mathcal{V}$commutes with the${\mathfrak{S}}_n$-action on$\mathcal{V}$. Both$\mathcal{V}$and${\mathcal{V}}^{\mathfrak{S}}$are graded$\mathfrak{gl}\left(1\right|1\left)\right[t]$-modules.

3.2. Properties of ${\mathcal{V}}^{\mathfrak{S}}$ and ${\left({\mathcal{V}}^{\mathfrak{S}}\right)}^{\mathrm{sing}}$

In this section, we recall properties of ${\mathcal{V}}^{\mathfrak{S}}$ and ${\left({\mathcal{V}}^{\mathfrak{S}}\right)}^{\mathrm{sing}}$ from [11] [Section 3].

Lemma 3.

The space${\left({\mathcal{V}}^{\mathfrak{S}}\right)}_{(n-\ell ,\ell )}$is a free$\mathbb{C}{[z_1,\dots ,z_n]}^{\mathfrak{S}}$-module of rank$\left(\genfrac{}{}{0pt}{}{n}{\ell }\right)$. In particular, the space${\mathcal{V}}^{\mathfrak{S}}$is a free$\mathbb{C}{[z_1,\dots ,z_n]}^{\mathfrak{S}}$-module of rank$2^n$.

Set $v^{+}=v_1^{\otimes n}=v_1\otimes \cdots \otimes v_1$.

Lemma 4.

The$\mathfrak{gl}\left(1\right|1\left)\right[t]$-module${\mathcal{V}}^{\mathfrak{S}}$is a cyclic module generated by$v^{+}$.

Lemma 5.

The set

$$\{e_{21}\left[r_1\right]e_{21}\left[r_2\right]\cdots e_{21}\left[r_{\ell }\right]v^{+}|0⩽r_1<r_2<\cdots <r_{\ell }⩽n-1\}$$

(16)

is a free generating set of${\left({\mathcal{V}}^{\mathfrak{S}}\right)}_{(n-\ell ,\ell )}$over$\mathbb{C}{[z_1,\dots ,z_n]}^{\mathfrak{S}}$.

Lemma 6.

The space${\left({\mathcal{V}}^{\mathfrak{S}}\right)}_{(n-\ell ,\ell )}^{\mathrm{sing}}$is a free$\mathbb{C}{[z_1,\dots ,z_n]}^{\mathfrak{S}}$-module of rank$\left(\genfrac{}{}{0pt}{}{n-1}{\ell }\right)$with a free generating set given by

$$\{e_{12}\left[0\right]e_{21}\left[0\right]e_{21}\left[r_1\right]\cdots e_{21}\left[r_{\ell }\right]v^{+},1⩽r_1<r_2<\cdots <r_{\ell }⩽n-1\}.$$

(17)

In particular, the space${\left({\mathcal{V}}^{\mathfrak{S}}\right)}^{\mathrm{sing}}$is a free$\mathbb{C}{[z_1,\dots ,z_n]}^{\mathfrak{S}}$-module of rank$2^{n-1}$.

Set ${\left(q\right)}_r={\prod }_{i=1}^r(1-q^i)$.

Proposition 2.

We have

$$\mathrm{ch}\left({\left({\mathcal{V}}^{\mathfrak{S}}\right)}_{(n-\ell ,\ell )}\right)=\frac{q^{\ell (\ell -1)/2}}{{\left(q\right)}_{\ell }{\left(q\right)}_{n-\ell }},\mathrm{ch}\left({\left({\mathcal{V}}^{\mathfrak{S}}\right)}_{(n-\ell ,\ell )}^{\mathrm{sing}}\right)=\frac{q^{\ell (\ell +1)/2}}{{\left(q\right)}_{\ell }{\left(q\right)}_{n-1-\ell }(1-q^n)}.$$

Given $\mathit{a}=(a_1,\dots ,a_n)\in {\mathbb{C}}^n$, let $I_{\mathit{a}}$ be the ideal of $\mathbb{C}{[z_1,\dots ,z_n]}^{\mathfrak{S}}$ generated by ${\sigma }_i\left(\mathit{z}\right)-\mathit{a}$, $i=1,\dots ,n$. Then, for any $\mathit{a}$, by Lemmas 2 and 3, the quotient space ${\mathcal{V}}^{\mathfrak{S}}/I_{\mathit{a}}{\mathcal{V}}^{\mathfrak{S}}$ is a $\mathfrak{gl}\left(1\right|1\left)\right[t]$-module of dimension $2^n$ over $\mathbb{C}$. Denote by ${\overline{v}}^{+}$ the image of $v^{+}$ under this quotient.

3.3. Weyl Modules

In this section, we recall a special family of Weyl modules for $\mathfrak{gl}\left(1\right|1\left)\right[t]$ and their properties from [10] [Section 3.3].

Let $\eta \left(x\right)$ be a monic polynomial of degree m with complex coefficients, where $m\in {\mathbb{Z}}_{⩾0}$,

$$\eta \left(x\right)=\sum _{i=0}^m{\gamma }_i{x}^i,{\gamma }_m=1.$$

Denote by $W_{\eta }$ the $\mathfrak{gl}\left(1\right|1\left)\right[t]$-module generated by an even vector w subject to the relations:

$$e_{11}\left(x\right)w={\eta }^{\prime }\left(x\right)/\eta \left(x\right)w,e_{22}\left(x\right)w=e_{12}\left(x\right)w=0,$$

(18)

$$\sum _{i=0}^m{\gamma }_i{e}_{21}\left[i\right]w=0.$$

(19)

It is convenient to write (19) as $(e_{21}\otimes \eta \left(t\right))w=0$.

Clearly, we have $dim{W}_{\eta }⩽2^m$ by the PBW theorem and (1), (19). The module $W_{\eta }$ is the universal $\mathfrak{gl}\left(1\right|1\left)\right[t]$-module satisfying (18), (19), which we call a Weyl module.

If $\eta \left(x\right)={(x-b)}^m$, we write $W_{\eta }$ as $W_m\left(b\right)$.

Lemma 7.

Let$\mathit{a}=(0,\dots ,0)\in {\mathbb{C}}^n$. Then,${\mathcal{V}}^{\mathfrak{S}}/I_{\mathit{a}}{\mathcal{V}}^{\mathfrak{S}}$is isomorphic to$W_n\left(0\right)$as$\mathfrak{gl}\left(1\right|1\left)\right[t]$-modules.

In particular, we have $dim{W}_m\left(b\right)=2^m$.

Lemma 8.

Let$\eta \left(x\right)={\prod }_{s=1}^k{(x-b_s)}^{n_s}$, where$b_s\ne b_r$for$1⩽s\ne r⩽k$. Then,$W_{\eta }$is isomorphic to${⨂}_{s=1}^k{W}_{n_s}\left(b_s\right)$as$\mathfrak{gl}\left(1\right|1\left)\right[t]$-modules.

Given sequences $\mathit{n}=(n_1,\dots ,n_k)$ of non-negative integers and $\mathit{b}=(b_1,\dots ,b_s)$ of distinct complex numbers, by Lemma 8, we call ${⨂}_{s=1}^k{W}_{n_s}\left(b_s\right)$ the Weyl module associated with $\mathit{n}$ and $\mathit{b}$.

Given $\mathit{a}=(a_1,\dots ,a_n)\in {\mathbb{C}}^n$, define $k\in {\mathbb{Z}}_{>0}$, $b_s\in \mathbb{C}$ and $n_s\in {\mathbb{Z}}_{>0}$ for $1⩽s⩽k$ by

$$x^n+\sum _{i=1}^n{(-1)}^i{a}_i{x}^{n-i}=\prod _{s=1}^k{(x-b_s)}^{n_s},$$

(20)

where $b_1,\dots ,b_k$ are distinct. Note that $n={\sum }_{s=1}^k{n}_s$.

Lemma 9.

The$\mathfrak{gl}\left(1\right|1\left)\right[t]$-module${\mathcal{V}}^{\mathfrak{S}}/I_{\mathit{a}}{\mathcal{V}}^{\mathfrak{S}}$is isomorphic to${⨂}_{s=1}^k{W}_{n_s}\left(b_s\right)$.

We also need the following statements.

Lemma 10.

Let $b\in \mathbb{C}$. We have the following properties for $W_m\left(b\right)$.

• As a$\mathfrak{gl}\left(1\right|1)$-module,$W_m\left(b\right)$is isomorphic to${\left({\mathbb{C}}^{1|1}\right)}^{\otimes m}$.
• A$\mathfrak{gl}\left(1\right|1\left)\right[t]$-module M is an irreducible subquotient of$W_m\left(b\right)$if and only if M has the form$L_{\mathit{\lambda }}\left(b\right)$, where$\mathit{\lambda }$is a polynomial weight such that$\left|\mathit{\lambda }\right|=m$.

Corollary 1.

A$\mathfrak{gl}\left(1\right|1\left)\right[t]$-module M is an irreducible subquotient of${⨂}_{s=1}^k{W}_{n_s}\left(b_s\right)$if and only if M has the form${⨂}_{s=1}^k{L}_{{\mathit{\lambda }}^{\left(s\right)}}\left(b_s\right)$, where${\mathit{\lambda }}^{\left(s\right)}$is a polynomial weight such that$|{\mathit{\lambda }}^{\left(s\right)}|=n_s$for each$1⩽s⩽k$.

4. Main Theorems

4.1. The Algebra ${\mathcal{O}}_l$

Let ${\Omega }_l$ be the n-dimensional affine space with coordinates $f_1,\dots ,f_l$, $g_1$, …, $g_{n-l}$. Introduce two polynomials

$$f\left(x\right)=x^l+\sum _{i=1}^l{f}_i{x}^{l-i},g\left(x\right)=x^{n-l}+\sum _{i=1}^{n-l}{g}_i{x}^{n-l-i}.$$

(21)

Denote by ${\mathcal{O}}_l$ the algebra of regular functions on ${\Omega }_l$, namely

$${\mathcal{O}}_l=\mathbb{C}[f_1,\dots ,f_l,g_1,\dots ,g_{n-l-1},g_{n-l}].$$

Define the degree function by

$$deg{f}_i=i,deg{g}_j=j,$$

for all $i=1,\dots ,l$ and $j=1,\dots ,n-l$. The algebra ${\mathcal{O}}_l$ is graded with the graded character given by

$$\mathrm{ch}\left({\mathcal{O}}_l\right)=\frac{1}{{\left(q\right)}_l{\left(q\right)}_{n-l}}.$$

(22)

Let ${\mathcal{F}}_0{\mathcal{O}}_l\subset {\mathcal{F}}_1{\mathcal{O}}_l\subset \cdots \subset {\mathcal{O}}_l$ be the increasing filtration corresponding to this grading, where ${\mathcal{F}}_s{\mathcal{O}}_l$ consists of elements of a degree of at most s.

Let ${\Sigma }_1,\dots ,{\Sigma }_n$ be the elements of ${\mathcal{O}}_l$ such that

$$(q_1-q_2)f\left(x\right)g\left(x\right)=(q_1-q_2)x^n+\sum _{i=1}^n{(-1)}^i\left((q_1-q_2){\Sigma }_i-(n+1-i){\Sigma }_{i-1}\right)x^{n-i},$$

(23)

where ${\Sigma }_0=1$. The homomorphism

$${\pi }_l:\mathbb{C}{[z_1,\dots ,z_n]}^{\mathfrak{S}}\to {\mathcal{O}}_l,{\sigma }_i\left(\mathit{z}\right)\mapsto {\Sigma }_i,i=1,\dots ,n,$$

(24)

is injective and induces a $\mathbb{C}{[z_1,\dots ,z_n]}^{\mathfrak{S}}$-module structure on ${\mathcal{O}}_l$.

Express $f^{\prime }\left(x\right)g\left(x\right)$ as follows:

$$(q_1-q_2)f^{\prime }\left(x\right)g\left(x\right)=(q_1-q_2)l{x}^{n-1}+\sum _{i=1}^{n-1}{G}_i{x}^{n-1-i},$$

(25)

where $G_i\in {\mathcal{O}}_l$.

Lemma 11.

The elements$G_i$and${\Sigma }_j$,$i=1,\dots ,n-1$,$j=1,\dots ,n$generate the algebra${\mathcal{O}}_l$.

Lemma 12.

We have$G_i\in {\mathcal{F}}_i{\mathcal{O}}_l\setminus {\mathcal{F}}_{i-1}{\mathcal{O}}_l$and${\Sigma }_j\in {\mathcal{F}}_j{\mathcal{O}}_l\setminus {\mathcal{F}}_{j-1}{\mathcal{O}}_l$,$i=1,\dots ,n-1$,$j=1,\dots ,n$.

4.2. Bethe Algebra

We call the unital subalgebra of $\mathrm{U}\left(\mathfrak{gl}\right(1\left|1\right)\left[t\right])$ generated by the coefficients of

$$e_{11}\left(x\right)+e_{22}\left(x\right),\mathcal{H}\left(x\right)=q_1{e}_{11}\left(x\right)+q_2{e}_{22}\left(x\right)+\frac{1}{2}\sum _{i,j=1}^2{e}_{ij}\left(x\right)e_{ji}\left(x\right){(-1)}^{\left|j\right|}$$

the Bethe algebra. We denote the Bethe algebra by $\mathcal{B}$. Note that the coefficients of $e_{11}\left(x\right)+e_{22}\left(x\right)$ generate the center of $\mathrm{U}\left(\mathfrak{gl}\right(1\left|1\right)\left[t\right])$.

Lemma 13 ([4]).

The Bethe algebra$\mathcal{B}$is commutative. The Bethe algebra$\mathcal{B}$commutes with the subalgebra$\mathrm{U}(\subset \mathrm{U}(\mathfrak{gl}\left(1\right|1\left)\right[t\left]\right)$.

Being a subalgebra of $\mathrm{U}\left(\mathfrak{gl}\right(1\left|1\right)\left[t\right])$, the Bethe algebra $\mathcal{B}$ acts on any $\mathfrak{gl}\left(1\right|1\left)\right[t]$-module M. Since $\mathcal{B}$ commutes with U(𝔥), the Bethe algebra preserves the subspace ${\left(M\right)}_{\mathit{\lambda }}$ for any weight $\mathit{\lambda }$. If $K\subset M$ is a $\mathcal{B}$-invariant subspace, then we call the image of $\mathcal{B}$ in $\mathrm{End}\left(K\right)$ the Bethe algebra associated withK.

Let $\mathit{a}=(a_1,\dots ,a_n)\in {\mathbb{C}}^n$. Define $k\in {\mathbb{Z}}_{>0}$, a sequence of positive integers $\mathit{n}=(n_1,\dots ,n_k)$ and a sequence of distinct complex numbers $\mathit{b}=(b_1,\dots ,b_k)$ by (20). Let $\mathbf{\Lambda }=({\mathit{\lambda }}^{\left(1\right)},\dots ,{\mathit{\lambda }}^{\left(k\right)})$ be a sequence of polynomial $\mathfrak{gl}\left(1\right|1)$-weights such that $|{\mathit{\lambda }}^{\left(s\right)}|=n_s$.

We study the action of the Bethe algebra $\mathcal{B}$ on the following $\mathcal{B}$-modules:

$${\mathcal{M}}_l={\left({\mathcal{V}}^{\mathfrak{S}}\right)}_{(n-l,l)},{\mathcal{M}}_{l,\mathit{a}}={\left(\underset{s=1}{\overset{k}{⨂}}{W}_{n_s}\left(b_s\right)\right)}_{(n-l,l)},{\mathcal{M}}_{l,\mathbf{\Lambda },\mathit{b}}={\left(\underset{s=1}{\overset{k}{⨂}}{L}_{{\mathit{\lambda }}^{\left(s\right)}}\left(b_s\right)\right)}_{(n-l,l)}.$$

Denote the Bethe algebras associated with ${\mathcal{M}}_l$, ${\mathcal{M}}_{l,\mathit{a}}$, ${\mathcal{M}}_{l,\mathbf{\Lambda },\mathit{b}}$ by ${\mathcal{B}}_l$, ${\mathcal{B}}_{l,\mathit{a}}$, ${\mathcal{B}}_{l,\mathbf{\Lambda },\mathit{b}}$, respectively. For any element $X\in \mathcal{B}$, we denote by $X\left(\mathit{z}\right)$, $X\left(\mathit{a}\right)$, $X(\mathbf{\Lambda },\mathit{b})$ the respective linear operators.

Since, by Lemma 4, the $\mathfrak{gl}\left(1\right|1\left)\right[t]$-module ${\mathcal{V}}^{\mathfrak{S}}$ is generated by $v_1^{\otimes n}=v_1\otimes \cdots \otimes v_1$, the series $e_{11}\left(x\right)+e_{22}\left(x\right)$ acts on ${\mathcal{V}}^{\mathfrak{S}}$ by multiplication by the series

$$\sum _{i=1}^n\frac{1}{x-z_i}=\sum _{i=1}^n\sum _{j=0}^{\infty }{z}_i^j{x}^{-j-1}.$$

Therefore, there exist unique central elements $C_1,\dots ,C_n$ of $\mathrm{U}\left(\mathfrak{gl}\right(1\left|1\right)\left[t\right])$ of minimal degrees such that each $C_i$ acts on ${\mathcal{V}}^{\mathfrak{S}}$ by multiplication by ${\sigma }_i\left(\mathit{z}\right)$.

Define $B_i\in \mathcal{B}$ by

$$\left(x^n+\sum _{i=1}^n{(-1)}^i{C}_i{x}^{n-i}\right)\mathcal{T}\left(x\right)=x^n\sum _{i=1}^{\infty }{B}_i{x}^{-i},$$

(26)

where $\mathcal{T}\left(x\right)$ is defined in (13).

Lemma 14.

We have$B_i\left(\mathit{z}\right)=0$for$i>n$and$B_1\left(\mathit{z}\right)=(q_1-q_2)l$.

Proof. 

Let $V\left(\mathit{c}\right)={⨂}_{i=1}^n{\mathbb{C}}^{1|1}\left(c_i\right)$, where $c_i\in \mathbb{C}$. Note that $B_i\left(\mathit{z}\right)$ is a polynomial in $\mathit{z}$ with values in $\mathrm{End}\left({\left(V\right)}_{(n-l,l)}\right)$. For any sequence of complex numbers $\mathit{c}=(c_1,\dots ,c_n)$, we can evaluate $B_i\left(\mathit{z}\right)$ at $\mathit{z}=\mathit{c}$ to an operator on ${\left(V\left(\mathit{c}\right)\right)}_{(n-l,l)}$. By Theorem 2, the Gaudin transfer matrix $\mathcal{H}\left(x\right)$ is diagonalizable and the Bethe ansatz is complete for ${\left(V\left(\mathit{c}\right)\right)}_{(n-l,l)}$ when $\mathit{c}\in {\mathbb{C}}^n$ is generic. Hence, by (14) and (26), $\left(x^n+{\sum }_{i=1}^n{(-1)}^i{C}_i{x}^{n-i}\right)\mathcal{T}\left(x\right)$ acts on ${\left(V\left(\mathit{c}\right)\right)}_{(n-l,l)}$ as a polynomial in x for generic $\mathit{c}$. In particular, it implies that $B_i$, $i>n$ acts on ${\left(V\left(\mathit{c}\right)\right)}_{(n-l,l)}$ by zero for generic $\mathit{c}$. Therefore, $B_i\left(\mathit{z}\right)$, $i>n$ is identically zero.

By the same reasoning, one shows that $B_1\left(\mathit{z}\right)=(q_1-q_2)l$. Alternatively, it also follows from $B_1=(q_1-q_2)e_{22}$. □

Lemma 15.

The elements$B_i\left(\mathit{z}\right)$and$C_j\left(\mathit{z}\right)$, for$1<i⩽n$and$1⩽j⩽n$, generate the algebra${\mathcal{B}}_l$.

Proof. 

It follows from the definition of $\mathcal{B}$, (26) and Lemma 14. □

One can restrict the filtration on $\mathrm{U}\left(\mathfrak{gl}\right(1\left|1\right)\left[t\right])$ to the Bethe algebra, ${\mathcal{F}}_0\mathcal{B}\subset {\mathcal{F}}_1\mathcal{B}\subset \cdots \subset \mathcal{B}$.

Lemma 16.

We have$B_i\in {\mathcal{F}}_{i-1}\mathcal{B}/{\mathcal{F}}_{i-2}\mathcal{B}$and$C_j\in {\mathcal{F}}_j\mathcal{B}/{\mathcal{F}}_{j-1}\mathcal{B}$for$1<i⩽n$and$1⩽j⩽n$.

4.3. Main Theorems

Recall from Proposition 2 that there exists a unique vector (up to proportionality) of degree $l(l-1)/2$ in ${\mathcal{M}}_l$ explicitly given by

$${\mathfrak{u}}_l:=e_{21}\left[0\right]e_{21}\left[1\right]\cdots e_{21}[l-1]v^{+};$$

see Lemma 5.

Any commutative algebra $\mathcal{A}$ is a module over itself induced by left multiplication. We call it the regular representation of $\mathcal{A}$. The dual space ${\mathcal{A}}^{\ast }$ is naturally an $\mathcal{A}$-module, which is called the coregular representation. A bilinear form $(·|·):\mathcal{A}\otimes \mathcal{A}\to \mathbb{C}$ is called invariant if $\left(ab\right|c)=(a\left|bc\right)$ for all $a,b,c\in \mathcal{A}$. A finite-dimensional commutative algebra $\mathcal{A}$ admitting an invariant non-degenerate symmetric bilinear form $(·|·):\mathcal{A}\otimes \mathcal{A}\to \mathbb{C}$ is called a Frobenius algebra. The regular and coregular representations of a Frobenius algebra are isomorphic.

Let M be an $\mathcal{A}$-module and $\mathcal{E}:\mathcal{A}\to \mathbb{C}$ a character; then, the $\mathcal{A}$-eigenspace associated to $\mathcal{E}$ in M is defined by ${\bigcap }_{a\in \mathcal{A}}ker(a{|}_M-\mathcal{E}\left(a\right))$. The generalized $\mathcal{A}$-eigenspace associated to $\mathcal{E}$ in M is defined by ${\bigcap }_{a\in \mathcal{A}}\left({\bigcup }_{m=1}^{\infty }ker{(a{|}_M-\mathcal{E}\left(a\right))}^m\right)$.

Theorem 3.

The action of the Bethe algebra${\mathcal{B}}_l$on${\mathcal{M}}_l$has the following properties.

• The map${\eta }_l:G_i\mapsto B_{i+1}\left(\mathit{z}\right)$,${\Sigma }_j\mapsto C_j\left(\mathit{z}\right)$,$i=1,\dots ,n-1$,$j=1,\dots ,n$extends uniquely to an isomorphism${\eta }_l:{\mathcal{O}}_l\to {\mathcal{B}}_l$of filtered algebras. Moreover, the isomorphism${\eta }_l$is an isomorphism of$\mathbb{C}{[z_1,\dots ,z_n]}^{\mathfrak{S}}$-modules.
• The map${\rho }_l:{\mathcal{O}}_l\mapsto {\mathcal{M}}_l$,$F\mapsto {\eta }_l\left(F\right){\mathfrak{u}}_l$is an isomorphism of filtered vector spaces identifying the${\mathcal{B}}_l$-module${\mathcal{M}}_l$with the regular representation of${\mathcal{O}}_l$.

Theorem 3 is proved in Section 5.

Let $\mathit{a}=(a_1,\dots ,a_n)\in {\mathbb{C}}^n$. Define $k\in {\mathbb{Z}}_{>0}$, a sequence of positive integers $\mathit{n}=(n_1,\dots ,n_k)$ and a sequence of distinct complex numbers $\mathit{b}=(b_1,\dots ,b_k)$ by (20). Let $\mathbf{\Lambda }=({\mathit{\lambda }}^{\left(1\right)},\dots ,{\mathit{\lambda }}^{\left(k\right)})$ be a sequence of non-degenerate polynomial weights such that $|{\mathit{\lambda }}^{\left(s\right)}|=n_s$ for each $1⩽s⩽k$.

Theorem 4.

The action of the Bethe algebra${\mathcal{B}}_{l,\mathbf{\Lambda },\mathit{b}}$on${\mathcal{M}}_{l,\mathbf{\Lambda },\mathit{b}}$has the following properties.

• The Bethe algebra${\mathcal{B}}_{l,\mathbf{\Lambda },\mathit{b}}$is isomorphic to

  $$\mathbb{C}{[w_1,\dots ,w_k]}^{{\mathfrak{S}}_l\times {\mathfrak{S}}_{k-l}}/{\langle {\sigma }_i\left(\mathit{w}\right)-{\epsilon }_i\rangle }_{i=1,\dots ,k},$$

  where${\epsilon }_i$is given by

  $${\phi }_{\mathbf{\Lambda },\mathit{b}}\left(x\right):=\prod _{s=1}^k(x-b_s)\left(q_1-q_2+\sum _{s=1}^k\frac{n_s}{x-b_s}\right)=(q_1-q_2)\left(x^k+\sum _{i=1}^k{(-1)}^i{\epsilon }_i{x}^{k-i}\right)$$

  and${\sigma }_i\left(\mathit{w}\right)$are elementary symmetric functions in$w_1,\dots ,w_k$.
• The Bethe algebra${\mathcal{B}}_{l,\mathbf{\Lambda },\mathit{b}}$is a Frobenius algebra. Moreover, the${\mathcal{B}}_{l,\mathbf{\Lambda },\mathit{b}}$-module${\mathcal{M}}_{l,\mathbf{\Lambda },\mathit{b}}$is isomorphic to the regular representation of${\mathcal{B}}_{l,\mathbf{\Lambda },\mathit{b}}$.
• The Bethe algebra${\mathcal{B}}_{l,\mathbf{\Lambda },\mathit{b}}$is a maximal commutative subalgebra in${\mathcal{M}}_{l,\mathbf{\Lambda },\mathit{b}}$of dimension$\left(\genfrac{}{}{0pt}{}{k}{l}\right)$.
• Every$\mathcal{B}$-eigenspace in${\mathcal{M}}_{l,\mathbf{\Lambda },\mathit{b}}$has dimension one.
• The$\mathcal{B}$-eigenspaces in${\mathcal{M}}_{l,\mathbf{\Lambda },\mathit{b}}$bijectively correspond to the monic degree l divisors$y\left(x\right)$of the polynomial${\phi }_{\mathbf{\Lambda },\mathit{b}}\left(x\right)$. Moreover, the eigenvalue of$\mathcal{H}\left(x\right)$corresponding to the monic divisor y is described by${\mathcal{E}}_{y,\mathbf{\Lambda },\mathit{b}}\left(x\right)$; see (12).
• Every generalized$\mathcal{B}$-eigenspace in${\mathcal{M}}_{l,\mathbf{\Lambda },\mathit{b}}$is a cyclic$\mathcal{B}$-module.
• The dimension of the generalized$\mathcal{B}$-eigenspace associated to${\mathcal{E}}_{y,\mathbf{\Lambda },\mathit{b}}\left(x\right)$is

  $$\prod _{a\in \mathbb{C}}\left(\genfrac{}{}{0pt}{}{{\mathrm{Mult}}_a\left({\phi }_{\mathbf{\Lambda },\mathit{b}}\right)}{{\mathrm{Mult}}_a\left(y\right)}\right),$$

  where${\mathrm{Mult}}_a\left(p\right)$is the multiplicity of a as a root of the polynomial p.

Theorem 4 is proved in Section 5.

Note that its results are quite parallel to that of XXX spin chains; see [11] [Theorem 4.11].

4.4. Higher Gaudin Transfer Matrices

To define higher Gaudin transfer matrices, we first recall basics about pseudo-differential operators. Let $\mathcal{A}$ be a differential superalgebra with an even derivation $\partial :\mathcal{A}\to \mathcal{A}$. For $r\in {\mathbb{Z}}_{>0}$, denote the r-th derivative of $a\in \mathcal{A}$ by $a_{\left[r\right]}$. Define the superalgebra of pseudo-differential operators $\mathcal{A}\left(\left({\partial }^{-1}\right)\right)$ as follows. Elements of $\mathcal{A}\left(\left({\partial }^{-1}\right)\right)$ are Laurent series in ${\partial }^{-1}$ with coefficients in $\mathcal{A}$, and the product is given by

$$\partial {\partial }^{-1}={\partial }^{-1}\partial =1,{\partial }^{r}a=\sum _{s=0}^{\infty }\left(\genfrac{}{}{0pt}{}{r}{s}\right)a_{\left[s\right]}{\partial }^{r-s},r\in \mathbb{Z},a\in \mathcal{A},$$

where

$$\left(\genfrac{}{}{0pt}{}{r}{s}\right)=\frac{r(r-1)\cdots (r-s+1)}{s!}.$$

Let

$${\mathcal{A}}_x^{m|n}=\mathrm{U}\left(\mathfrak{gl}\left(1\right|1)\left[t\right]\right)\left(\left(x^{-1}\right)\right)=\left\{\sum _{r=-\infty }^s{g}_r{x}^r,r\in \mathbb{Z},g_r\in \mathrm{U}\left(\mathfrak{gl}\left(1\right|1)\left[t\right]\right)\right\}.$$

Consider the operator in $\mathrm{End}\left({\mathbb{C}}^{1|1}\right)\otimes {\mathcal{A}}_x^{m|n}\left(\left({\partial }_x^{-1}\right)\right)$,

$$\mathfrak{Z}(x,{\partial }_x):=\sum _{a,b=1}^2{E}_{ab}\otimes \left({\delta }_{ab}({\partial }_x-q_a)-e_{ab}\left(x\right){(-1)}^{\left|a\right|}\right),$$

which is a Manin matrix; see [4] [Lemma 3.1] and [8] [Lemma 4.2]. Define the Berezinian—see [14]—of $\mathfrak{Z}(x,{\partial }_x)$ by

$$\begin{array}{cc}& \mathrm{Ber}\left(\mathfrak{Z}(x,{\partial }_x)\right)\hfill \\ \hfill =& \left({\partial }_x-q_1-e_{11}\left(x\right)\right){\left({\partial }_x-q_2+e_{22}\left(x\right)+e_{21}\left(x\right){\left({\partial }_x-q_1-e_{11}\left(x\right)\right)}^{-1}{e}_{12}\left(x\right)\right)}^{-1}.\hfill \end{array}$$

(27)

Denote the Berezinian by $\mathfrak{D}(x,{\partial }_x)$ and expand it as an element in ${\mathcal{A}}_x^{m|n}\left(\left({\partial }_x^{-1}\right)\right)$,

$$\mathfrak{D}(x,{\partial }_x)=\sum _{r=0}^{\infty }{(-1)}^r{\mathcal{G}}_r\left(x\right){\partial }_x^{-r}.$$

(28)

We call the series ${\mathcal{G}}_r\left(x\right)\in {\mathcal{A}}_x^{m|n}$, $r\in {\mathbb{Z}}_{⩾0}$ the higher Gaudin transfer matrices. In particular, we call ${\mathcal{G}}_1\left(x\right)$ and ${\mathcal{G}}_2\left(x\right)$ the first and second Gaudin transfer matrices, respectively.

Example 1.

We have ${\mathcal{G}}_0\left(x\right)=1$,

$${\mathcal{G}}_1\left(x\right)=q_1-q_2+e_{11}\left(x\right)+e_{22}\left(x\right),$$

$${\mathcal{G}}_2\left(x\right)=\left(q_1-q_2+e_{11}\left(x\right)+e_{22}\left(x\right)\right)(-q_2+e_{22}\left(x\right))-e_{21}\left(x\right)e_{12}\left(x\right).$$

Remark 1.

In principle, the Bethe algebra should be the unital subalgebra of$\mathrm{U}\left(\mathfrak{gl}\right(1\left|1\right)\left[t\right])$generated by coefficients${\mathcal{G}}_r\left(x\right)$,$r\in {\mathbb{Z}}_{>0}$, cf. [15]. However, it turns out that the first two transfer matrices already give (almost) complete information about the Bethe algebra; see the discussion below.

Now, we describe the eigenvalues of higher Gaudin transfer matrices acting on the on-shell Bethe vector.

Let $\mathbf{\Lambda }=({\mathit{\lambda }}^{\left(1\right)},\dots ,{\mathit{\lambda }}^{\left(k\right)})$ be a sequence of $\mathfrak{gl}\left(1\right|1)$-weights and $\mathit{b}=(b_1,\dots ,b_k)$ a sequence of distinct complex numbers, where ${\mathit{\lambda }}^{\left(s\right)}=({\alpha }_s,{\beta }_s)$. Let $\mathit{t}=(t_1,\dots ,t_l)$, where $0⩽l<k$. Suppose that $y_{\mathit{t}}$ divides the polynomial ${\phi }_{\mathbf{\Lambda },\mathit{b}}$ (namely $\mathit{t}$ satisfies the Bethe ansatz equation); see (10).

Theorem 5 ([16] [Theorem 5.2]).

If$t_i\ne t_j$for$1⩽i<j⩽l$, then

$$\mathfrak{D}(x,{\partial }_x){\mathbb{B}}_l\left(\mathit{t}\right)={\mathbb{B}}_l\left(\mathit{t}\right)\left({\partial }_x-q_1-\sum _{s=1}^k\frac{{\alpha }_s}{x-b_s}+\frac{y_{\mathit{t}}^{\prime }}{y_{\mathit{t}}}\right){\left({\partial }_x-q_2+\sum _{s=1}^k\frac{{\beta }_s}{x-b_s}+\frac{y_{\mathit{t}}^{\prime }}{y_{\mathit{t}}}\right)}^{-1}.$$

(29)

The theorem is a differential analog of [11] [Theorem 6.4]. Note that the pseudo-differential operator in the right-hand side of (29), denoted by ${\mathfrak{D}}_{y,\mathbf{\Lambda },\mathit{b}}$, was introduced [7] [Section 5.3]. This theorem is generalized to the $\mathfrak{gl}\left(m\right|n)$ case in [16] [Theorem 5.2] where, on the right-hand side, the pseudo-differential operator describing the eigenvalues of higher Gaudin transfer matrices should be replaced by the pseudo-differential operator in [7] [Equation (6.5)]. This generalization is a classical limit of [17] [Conjecture 5.15] and [16] [Corollary 3.6] that connects the rational difference operator introduced in [18] [Equation (5.6)] with the eigenvalues of higher transfer matrices on the on-shell Bethe vector for XXX spin chains associated with $\mathfrak{gl}\left(m\right|n)$. The method used in the proof of [16] [Theorem 5.2] is motivated by [19,20] via the nested algebraic Bethe ansatz introduced in [21].

Remark 2.

As shown in [7] [Lemma 5.7], the odd reflection of${\mathfrak{D}}_{y,\mathbf{\Lambda },\mathit{b}}$, cf. [7] [Equation (3.1)], which comes from the study of the fermionic reproduction procedure of the Bethe ansatz equation, is compatible with the odd reflection of Lie superalgebras. The difference analog of this fact was used in [22] to investigate the relations between the odd reflections of the super Yangian of type A and the fermionic reproduction procedure of the Bethe ansatz equation for XXX spin chains.

We conclude this section by discussing the connections between ${\mathcal{G}}_i\left(x\right)$, $i⩾3$ and ${\mathcal{G}}_1\left(x\right)$, ${\mathcal{G}}_2\left(x\right)$.

Let

$$\mu \left(x\right)=q_1+\sum _{s=1}^k\frac{{\alpha }_s}{x-b_s}-\frac{y_{\mathit{t}}^{\prime }}{y_{\mathit{t}}},\nu \left(x\right)=-q_2+\sum _{s=1}^k\frac{{\beta }_s}{x-b_s}+\frac{y_{\mathit{t}}^{\prime }}{y_{\mathit{t}}}.$$

For simplicity, we do not write the dependence of $\mu \left(x\right)$ and $\nu \left(x\right)$ on $\mathbf{\Lambda },\mathit{b},\mathit{t}$ explicitly. Then, the eigenvalue of $\mathfrak{D}(x,{\partial }_x)$ acting on ${\mathbb{B}}_l\left(\mathit{t}\right)$ is given by

$$({\partial }_x-\mu \left(x\right)){({\partial }_x+\nu \left(x\right))}^{-1}=1-(\mu \left(x\right)+\nu \left(x\right)){({\partial }_x+\nu \left(x\right))}^{-1}.$$

(30)

Hence, the eigenvalues of ${\mathcal{G}}_i\left(x\right)$ are essentially only determined by $\mu \left(x\right)+\nu \left(x\right)$ and $\nu \left(x\right)$. Comparing (28) and the expansion of (30), we have

$${\mathcal{G}}_1\left(x\right){\mathbb{B}}_l\left(\mathit{t}\right)=(\mu \left(x\right)+\nu \left(x\right)){\mathbb{B}}_l\left(\mathit{t}\right),{\mathcal{G}}_2\left(x\right){\mathbb{B}}_l\left(\mathit{t}\right)=(\mu \left(x\right)+\nu \left(x\right))\nu \left(x\right){\mathbb{B}}_l\left(\mathit{t}\right);$$

(31)

see also (14). Therefore, the spectrum of all higher transfer matrices are simply determined by that of the first two transfer matrices, which justifies our definition of Bethe algebra.

Lemma 17.

Let the complex parameters$c_1,\dots ,c_m$and the positive integer m vary. Then, the kernels of the representations${⨂}_{i=1}^m{\mathbb{C}}^{1|1}\left(c_i\right)$of$\mathrm{U}\left(\mathfrak{gl}\right(1\left|1\right)\left[t\right])$have a zero intersection.

Proof. 

The proof is contained in the proof of [23] [Proposition 1.7]. □

Corollary 2.

We have

$$\mathfrak{D}(x,{\partial }_x)=\left({\partial }_x-{\mathcal{G}}_1\left(x\right)+\frac{{\mathcal{G}}_2\left(x\right)}{{\mathcal{G}}_1\left(x\right)}\right){\left({\partial }_x+\frac{{\mathcal{G}}_2\left(x\right)}{{\mathcal{G}}_1\left(x\right)}\right)}^{-1}.$$

(32)

Proof. 

By Lemma 17, it suffices to check that the left-hand side and the right-hand side of (32) act identically on a basis of ${⨂}_{i=1}^m{\mathbb{C}}^{1|1}\left(c_i\right)$ for all $m\in {\mathbb{Z}}_{>0}$ and generic $\mathit{c}=(c_1,\dots ,c_m)$.

By Theorem 2, there is a basis of ${⨂}_{i=1}^m{\mathbb{C}}^{1|1}\left(c_i\right)$ consisting of on-shell Bethe vectors for generic $\mathit{c}$. Therefore, the statement follows from Theorem 5, (30) and (31). □

5. Proof of Main Theorems

In this section, we prove the main theorems. For completeness, we provide all details, even if they are parallel to those in [10] [Section 5].

5.1. The First Isomorphism

Proof of Theorem 3.

We first show that the homomorphism defined by ${\eta }_l$ is well-defined.

Consider the tensor product $V\left(\mathit{c}\right)={⨂}_{i=1}^n{\mathbb{C}}^{1|1}\left(c_i\right)$, where $c_i\in \mathbb{C}$, and the corresponding Bethe ansatz equation associated to weight $(n-l,l)$. Let $\mathit{t}$ be a solution with distinct coordinates and ${\mathbb{B}}_l\left(\mathit{t}\right)$ be the corresponding on-shell Bethe vector. Denote ${\mathcal{E}}_{i,\mathit{t}}$ the eigenvalues of $B_i$ acting on ${\mathbb{B}}_l\left(\mathit{t}\right)$; see Theorem 1 and Equation (14).

Define a character $\pi :{\mathcal{O}}_l\to \mathbb{C}$ by sending

$$f\left(x\right)\mapsto y_{\mathit{t}}\left(x\right),{\Sigma }_n\mapsto \prod _{i=1}^n{c}_i,$$

$$g\left(x\right)\mapsto \frac{1}{(q_1-q_2)y_{\mathit{t}}\left(x\right)}\prod _{i=1}^n(x-c_i)\left(q_1-q_2+\sum _{i=1}^n\frac{1}{x-c_i}\right).$$

Then,

$$\pi \left({\Sigma }_i\right)={\sigma }_i\left(\mathit{c}\right),\pi \left(G_i\right)={\mathcal{E}}_{i,\mathit{t}},$$

(33)

by (23) and by (12), (14), (25), respectively.

Now, let $P(G_i,{\Sigma }_j)$ be a polynomial in $G_i,{\Sigma }_j$ such that $P(G_i,{\Sigma }_j)$ is equal to zero in ${\mathcal{O}}_l$. It suffices to show that $P(B_i\left(\mathit{z}\right),C_j\left(\mathit{z}\right))$ is equal to zero in ${\mathcal{B}}_l$.

Note that $P(B_i\left(\mathit{z}\right),C_j\left(\mathit{z}\right))$ is a polynomial in $z_1,\dots ,z_n$ with values in $\mathrm{End}\left({\left(V\right)}_{(n-l,l)}\right)$. For any sequence $\mathit{c}$ of complex numbers, we can evaluate $P(B_i\left(\mathit{z}\right),C_j\left(\mathit{z}\right))$ at $\mathit{z}=\mathit{c}$ to an operator on ${\left(V\left(\mathit{c}\right)\right)}_{(n-l,l)}$. By Theorem 2, the transfer matrix $\mathcal{T}\left(x\right)$ is diagonalizable and the Bethe ansatz is complete for ${\left(V\left(\mathit{c}\right)\right)}_{(n-l,l)}$ when $\mathit{c}\in {\mathbb{C}}^n$ is generic. Hence, by (33), the value of $P(B_i\left(\mathit{z}\right),C_j\left(\mathit{z}\right))$ at $\mathit{z}=\mathit{c}$ is also equal to zero for generic $\mathit{c}$. Therefore, $P(B_i\left(\mathit{z}\right),C_j\left(\mathit{z}\right))$ is identically zero and the map ${\eta }_l$ is well-defined.

Let us now show that the map ${\eta }_l$ is injective. Let $P(G_i,{\Sigma }_j)$ be a polynomial in $G_i,{\Sigma }_j$ such that $P(G_i,{\Sigma }_j)$ is non-zero in ${\mathcal{O}}_l$. Then, the value at a generic point of ${\Omega }_l$ (e.g., the non-vanishing points of $P(G_i,{\Sigma }_j)$ such that f and g are relatively prime and have only simple zeros) is not equal to zero. Moreover, at those points, the transfer matrix $\mathcal{T}\left(x\right)$ is diagonalizable and the Bethe ansatz is complete again by Theorem 2. Therefore, again by (33), the polynomial $P(B_i\left(\mathit{z}\right),C_j\left(\mathit{z}\right))$ is a non-zero element in ${\mathcal{B}}_l$. Thus, the map ${\eta }_l$ is injective.

The surjectivity of ${\eta }_l$ follows from Lemma 15. Hence, ${\eta }_l$ is an isomorphism of algebras.

The fact that ${\eta }_l$ is an isomorphism of graded algebra respecting the gradation follows from Lemmas 12 and 16. This completes the proof of part (i).

The kernel of ${\rho }_l$ is an ideal of ${\mathcal{O}}_l$. If we identify ${\sigma }_i\left(\mathit{z}\right)$ with ${\Sigma }_i$, then the algebra ${\mathcal{O}}_l$ contains the algebra $\mathbb{C}{[z_1,\dots ,z_n]}^{\mathfrak{S}}$; see (24). The kernel of ${\rho }_l$ intersects $\mathbb{C}{[z_1,\dots ,z_n]}^{\mathfrak{S}}$ trivially. Therefore, the kernel of ${\rho }_l$ is trivial as well. Hence, ${\rho }_l$ is an injective map. Comparing Equation (22) and Proposition 2, we have $\mathrm{ch}\left({\mathcal{M}}_l\right)=q^{l(l-1)/2}\mathrm{ch}\left({\mathcal{O}}_l\right)$. Thus, ${\rho }_l$ is an isomorphism of graded vector spaces, which shifts the degree by $l(l-1)/2$, completing the proof of part (ii). □

5.2. The Second Isomorphism

Let $\mathit{a}=(a_1,\dots ,a_n)$ be a sequence of complex numbers. Define $k\in {\mathbb{Z}}_{>0}$, a sequence of positive integers $\mathit{n}=(n_1,\dots ,n_k)$, and a sequence of distinct complex numbers $\mathit{b}=(b_1,\dots ,b_k)$ by (20). Let $I_{l,\mathit{a}}^{\mathcal{O}}$ be the ideal of ${\mathcal{O}}_l$ generated by the elements ${\Sigma }_i-a_i$, $i=1,\dots ,n$, where ${\Sigma }_1,\dots ,{\Sigma }_{n-1}$ are defined in (23). Let ${\mathcal{O}}_{l,\mathit{a}}$ be the quotient algebra

$${\mathcal{O}}_{l,\mathit{a}}={\mathcal{O}}_l/I_{l,\mathit{a}}^{\mathcal{O}}.$$

Let $I_{l,\mathit{a}}^{\mathcal{B}}$ be the ideal of ${\mathcal{B}}_l$ generated by $C_i\left(\mathit{z}\right)-a_i$, $i=1,\dots ,n$. Consider the subspace

$$I_{l,\mathit{a}}^{\mathcal{M}}=I_{l,\mathit{a}}^{\mathcal{B}}{\mathcal{M}}_l={\left(I_{\mathit{a}}{\mathcal{V}}^{\mathfrak{S}}\right)}_{(n-l,l)},$$

where $I_{\mathit{a}}$, as before, is the ideal of $\mathbb{C}{[z_1,\dots ,z_n]}^{\mathfrak{S}}$ generated by ${\sigma }_i\left(\mathit{z}\right)-a_i$, $i=1,\dots ,n$.

Lemma 18.

We have

$${\eta }_l\left(I_{l,\mathit{a}}^{\mathcal{O}}\right)=I_{l,\mathit{a}}^{\mathcal{B}},{\rho }_l\left(I_{l,\mathit{a}}^{\mathcal{O}}\right)=I_{l,\mathit{a}}^{\mathcal{M}},{\mathcal{B}}_{l,\mathit{a}}={\mathcal{B}}_l/I_{l,\mathit{a}}^{\mathcal{B}},{\mathcal{M}}_{l,\mathit{a}}={\left({\mathcal{V}}^{\mathfrak{S}}\right)}_{(n-l,l)}^{\mathrm{sing}}/I_{l,\mathit{a}}^{\mathcal{M}}.$$

Proof. 

The lemma follows from Theorem 3 and Lemma 9. □

By Lemma 18, the maps ${\eta }_l$ and ${\rho }_l$ induce the maps

$${\eta }_{l,\mathit{a}}:{\mathcal{O}}_{l,\mathit{a}}\to {\mathcal{B}}_{l,\mathit{a}},{\rho }_{l,\mathit{a}}:{\mathcal{O}}_{l,\mathit{a}}\to {\mathcal{M}}_{l,\mathit{a}}.$$

The map ${\eta }_{l,\mathit{a}}$ is an isomorphism of algebras. Since ${\mathcal{B}}_{l,\mathit{a}}$ is finite-dimensional, by, e.g., [1] [Lemma 3.9], ${\mathcal{O}}_{l,\mathit{a}}$ is a Frobenius algebra, and so is ${\mathcal{B}}_{l,\mathit{a}}$. The map ${\rho }_{l,\mathit{a}}$ is an isomorphism of vector spaces. Moreover, it follows from Theorem 3 and Lemma 18 that ${\rho }_{l,\mathit{a}}$ identifies the regular representation of ${\mathcal{O}}_{l,\mathit{a}}$ with the ${\mathcal{B}}_{l,\mathit{a}}$-module ${\mathcal{M}}_{l,\mathit{a}}$.

The statement of this section implies, by, e.g., [13] [Lemma 1.3], the following. Set

$${\zeta }_{\mathit{n},\mathit{b}}\left(x\right)=q_1-q_2+\sum _{s=1}^k\frac{n_s}{x-b_s},{\psi }_{\mathit{n},\mathit{b}}\left(x\right):={\zeta }_{\mathit{n},\mathit{b}}\left(x\right)\prod _{r=1}^k{(x-b_r)}^{n_s}.$$

Theorem 6.

Suppose that$\mathit{b}=(b_1,\dots ,b_k)$is a sequence of distinct complex numbers. Then, the Gaudin transfer matrix$\mathcal{H}\left(x\right)$has a simple spectrum in${\left({⨂}_{s=1}^k{W}_{n_s}\left(b_s\right)\right)}^{\mathrm{sing}}$. There exists a bijective correspondence between the monic divisors y of the polynomial${\psi }_{\mathit{n},\mathit{b}}$and the eigenvectors$v_y$of the Gaudin transfer matrix$\mathcal{H}\left(x\right)$(up to multiplication by a non-zero constant). Moreover, this bijection is such that

$$\mathcal{H}\left(x\right)v_y=\left(\frac{1}{2}{\zeta }_{\mathit{n},\mathit{b}}^{\prime }\left(x\right)-{\zeta }_{\mathit{n},\mathit{b}}\left(x\right)\frac{y^{\prime }\left(x\right)}{y\left(x\right)}+\frac{1}{2}\left(\sum _{s=1}^k\frac{n_s}{x-b_s}\right)\left(\sum _{s=1}^k\frac{n_s}{x-b_s}+2q_1\right)\right)v_y.$$

Remark 3.

Fix$l\in {\mathbb{Z}}_{⩾0}$and set$\mathit{t}=(t_1,\dots ,t_l)$. Let${\mathbf{y}}_{\mathit{t}}$represent$\mathit{t}$. Then, the Bethe ansatz equation for$\left({⨂}_{s=1}^k{W}_{n_s}\left(b_s\right)\right)$is

$$y_{\mathit{t}}\left(x\right)\mathit{divides}\mathit{the}\mathit{polynomial}{\psi }_{\mathit{n},\mathit{b}}\left(x\right).$$

Note that, in this case,$y_{\mathit{t}}$may have multiple roots. If there are multiple roots in$y_{\mathit{t}}$, then the corresponding on-shell Bethe vector is zero. Therefore, an actual eigenvector should be obtained via an appropriate derivative as pointed out in [7] [Section 8.2].

5.3. The Third Isomorphism

Recall from Section 2.3 that, without a loss of generality, we can assume that ${\beta }_s=0$, $1⩽s⩽k$. In this case, ${\alpha }_s=n_s$, $1⩽s⩽k$.

Lemma 19.

There exists a surjective$\mathfrak{gl}\left(1\right|1\left)\right[t]$-module homomorphism from${⨂}_{s=1}^k{W}_{n_s}\left(b_k\right)$to${⨂}_{s=1}^k{L}_{{\mathit{\lambda }}^{\left(s\right)}}\left(b_k\right)$that maps a vacuum vector to a vacuum vector.

Proof. 

It follows from Lemma 10 and our assumption that ${\beta }_s=0$ for all $1⩽s⩽k$. □

By Lemma 9, the surjective $\mathfrak{gl}\left(1\right|1\left)\right[t]$-module homomorphism

$$\underset{s=1}{\overset{k}{⨂}}{W}_{n_s}\left(b_k\right)\twoheadrightarrow \underset{s=1}{\overset{k}{⨂}}{L}_{{\mathit{\lambda }}^{\left(s\right)}}\left(b_k\right)$$

induces a surjective $\mathfrak{gl}\left(1\right|1\left)\right[t]$-module homomorphism

$${\mathcal{V}}^{\mathfrak{S}}\twoheadrightarrow \underset{s=1}{\overset{k}{⨂}}{L}_{{\mathit{\lambda }}^{\left(s\right)}}\left(b_k\right).$$

The second map then induces a projection of the Bethe algebras ${\mathcal{B}}_l\twoheadrightarrow {\mathcal{B}}_{l,\mathbf{\Lambda },\mathit{b}}$. We describe the kernel of this projection. We consider the corresponding ideal in the algebra ${\mathcal{O}}_l$.

Suppose that $l⩽k$. Define the polynomial $h\left(x\right)$ by

$$h\left(x\right)=\prod _{s=1}^k{(x-b_s)}^{n_s-1}.$$

Divide the polynomial $g\left(x\right)$ in (21) by $h\left(x\right)$ and let

$$p\left(x\right)=x^{k-l}+p_1{x}^{k-l-1}+\cdots +p_{k-l-1}x+p_{k-l},$$

(34)

$$r\left(x\right)=r_1{x}^{n-k-1}+r_2{x}^{n-k-2}+\cdots +r_{n-k-1}x+r_{n-k}$$

(35)

be the quotient and the remainder, respectively. Clearly, $p_i,r_j\in {\mathcal{O}}_l$.

Denote by $I_{l,\mathbf{\Lambda },\mathit{b}}^{\mathcal{O}}$ the ideal of ${\mathcal{O}}_l$ generated by $r_1,\dots ,r_{n-k}$ and the coefficients of polynomial

$${\phi }_{\mathbf{\Lambda },\mathit{b}}\left(x\right)-(q_1-q_2)p\left(x\right)f\left(x\right)=\prod _{s=1}^k(x-b_s)\left(q_1-q_2+\sum _{s=1}^k\frac{n_s}{x-b_s}\right)-(q_1-q_2)p\left(x\right)f\left(x\right).$$

Let ${\mathcal{O}}_{l,\mathbf{\Lambda },\mathit{b}}$ be the quotient algebra

$${\mathcal{O}}_{l,\mathbf{\Lambda },\mathit{b}}={\mathcal{O}}_l/I_{l,\mathbf{\Lambda },\mathit{b}}^{\mathcal{O}}.$$

Clearly, if ${\mathcal{O}}_{l,\mathbf{\Lambda },\mathit{b}}$ is finite-dimensional, then it is a Frobenius algebra.

Let $I_{l,\mathbf{\Lambda },\mathit{b}}^{\mathcal{B}}$ be the image of $I_{l,\mathbf{\Lambda },\mathit{b}}^{\mathcal{O}}$ under the isomorphism ${\eta }_l$.

Lemma 20.

The ideal$I_{l,\mathbf{\Lambda },\mathit{b}}^{\mathcal{B}}$is contained in the kernel of the projection${\mathcal{B}}_l\twoheadrightarrow {\mathcal{B}}_{l,\mathbf{\Lambda },\mathit{b}}$.

Proof. 

We treat $\mathit{b}=(b_1,\dots ,b_k)$ as variables. Note that the elements of $I_{l,\mathbf{\Lambda },\mathit{b}}^{\mathcal{B}}$ act on ${\mathcal{M}}_{l,\mathbf{\Lambda },\mathit{b}}$ as polynomials in $\mathit{b}$ with values in $\mathrm{End}\left({\left(L_{\mathbf{\Lambda }}\right)}_{(n-l,l)}\right)$. Therefore it suffices to show it for generic $\mathit{b}$. Let $\mathfrak{f}\left(x\right)$ be the image of $f\left(x\right)$ under ${\eta }_l$. The condition that $I_{l,\mathbf{\Lambda },\mathit{b}}^{\mathcal{B}}$ vanishes is equivalent to the condition that ${\phi }_{\mathbf{\Lambda },\mathit{b}}\left(x\right)$ is divisible by $\mathfrak{f}\left(x\right)$.

By Theorem 2, there exists an eigenbasis of the operator $\mathcal{T}\left(x\right)$ in ${\mathcal{M}}_{l,\mathbf{\Lambda },\mathit{b}}$ for generic $\mathit{b}$. Clearly, a solution of the Bethe ansatz equation associated to $\mathbf{\Lambda },\mathit{b},l$ is also a solution to the Bethe ansatz equation for ${\mathcal{M}}_{l,\mathit{a}}$; see Theorem 6 and Remark 3. Moreover, the expressions of corresponding on-shell Bethe vectors coincide (with different vacuum vectors). By Lemma 19 and Theorems 1 and 6, ${\phi }_{\mathbf{\Lambda },\mathit{b}}\left(x\right)$ is divisible by $\mathfrak{f}\left(x\right)$ for generic $\mathit{b}$ since the eigenvalue of $\mathfrak{f}\left(x\right)$ corresponds to $y_{\mathit{t}}\left(x\right)$ in (12). Therefore, $I_{l,\mathbf{\Lambda },\mathit{b}}^{\mathcal{B}}$ vanishes for generic $\mathit{b}$, thus completing the proof. □

Therefore, we have the epimorphism

$${\mathcal{O}}_{l,\mathbf{\Lambda },\mathit{b}}\cong {\mathcal{B}}_l/I_{l,\mathbf{\Lambda },\mathit{b}}^{\mathcal{B}}\twoheadrightarrow {\mathcal{B}}_{l,\mathbf{\Lambda },\mathit{b}}.$$

(36)

We claim that the surjection in (36) is an isomorphism by checking $dim{\mathcal{O}}_{l,\mathbf{\Lambda },\mathit{b}}=dim{\mathcal{B}}_{l,\mathbf{\Lambda },\mathit{b}}$.

Lemma 21.

We have${\displaystyle dim{\mathcal{O}}_{l,\mathbf{\Lambda },\mathit{b}}=\left(\genfrac{}{}{0pt}{}{k}{l}\right)}$.

Proof. 

Note that $\mathbb{C}[p_1,\dots ,p_{k-l},r_1,\dots ,r_{n-k}]\cong \mathbb{C}[g_1,\dots ,g_{n-l}]$, where $p_i$ and $r_j$ are defined in (34) and (35). It is not hard to check that

$${\mathcal{O}}_{l,\mathbf{\Lambda },\mathit{b}}\cong \mathbb{C}[f_1,\dots ,f_l,p_1,\dots ,p_{k-l}]/{\tilde{I}}_{l,\mathbf{\Lambda },\mathit{b}}^{\mathcal{O}},$$

(37)

where ${\tilde{I}}_{l,\mathbf{\Lambda },\mathit{b}}^{\mathcal{O}}$ is the ideal of $\mathbb{C}[f_1,\dots ,f_l,p_1,\dots ,p_{k-l}]$ generated by the coefficients of the polynomial ${\phi }_{\mathbf{\Lambda },\mathit{b}}\left(x\right)-(q_1-q_2)p\left(x\right)f\left(x\right)$.

Introduce new variables $\mathit{w}=(w_1,\dots ,w_k)$ such that

$$f\left(x\right)=\prod _{i=1}^l(x-w_i),p\left(x\right)=\prod _{i=1}^{k-l}(x-w_{l+i}).$$

Let $\mathbf{\epsilon }=({\epsilon }_1,\dots ,{\epsilon }_k)$ be complex numbers such that

$${\phi }_{\mathbf{\Lambda },\mathit{b}}\left(x\right)=\prod _{s=1}^k(x-b_s)\left(q_1-q_2+\sum _{s=1}^k\frac{n_s}{x-b_s}\right)=(q_1-q_2)\left(x^k+\sum _{i=1}^k{(-1)}^i{\epsilon }_i{x}^{k-i}\right).$$

Then,

$$\mathbb{C}[f_1,\dots ,f_l,p_1,\dots ,p_{k-l}]/{\tilde{I}}_{l,\mathbf{\Lambda },\mathit{b}}^{\mathcal{O}}\cong \mathbb{C}{[w_1,\dots ,w_{k-1}]}^{{\mathfrak{S}}_l\times {\mathfrak{S}}_{k-l}}/{\langle {\sigma }_i\left(\mathit{w}\right)-{\epsilon }_i\rangle }_{i=1,\dots ,k}.$$

(38)

The lemma now follows from the fact that $\mathbb{C}{[w_1,\dots ,w_{k-1}]}^{{\mathfrak{S}}_l\times {\mathfrak{S}}_{k-l}}$ is a free $\mathbb{C}{[w_1,\dots ,w_{k-1}]}^{\mathfrak{S}}$-module of rank $\left(\genfrac{}{}{0pt}{}{k}{l}\right)$. □

Note that we have the projection ${\left({\mathcal{V}}^{\mathfrak{S}}\right)}_{(n-l,l)}\twoheadrightarrow {\mathcal{M}}_{l,\mathbf{\Lambda },\mathit{b}}$. Since, by Theorem 3, the Bethe algebra ${\mathcal{B}}_l$ acts on ${\left({\mathcal{V}}^{\mathfrak{S}}\right)}_{(n-l,l)}$ cyclically, the Bethe algebra ${\mathcal{B}}_{l,\mathbf{\Lambda },\mathit{b}}$ acts on ${\mathcal{M}}_{l,\mathbf{\Lambda },\mathit{b}}$ cyclically as well. Therefore, we have

$$dim{\mathcal{B}}_{l,\mathbf{\Lambda },\mathit{b}}=dim{\mathcal{M}}_{l,\mathbf{\Lambda },\mathit{b}}=\left(\genfrac{}{}{0pt}{}{k}{l}\right).$$

(39)

Proof of Theorem 4.

Part (i) follows from Lemma 21 and (36)–(39). Clearly, we have that ${\mathcal{B}}_{l,\mathbf{\Lambda },\mathit{b}}\cong {\mathcal{O}}_{l,\mathbf{\Lambda },\mathit{b}}$ is a Frobenius algebra. Moreover, the map ${\rho }_l$ from Theorem 3 induces a map

$${\rho }_{l,\mathbf{\Lambda },\mathit{b}}:{\mathcal{O}}_{l,\mathbf{\Lambda },\mathit{b}}\to {\mathcal{M}}_{l,\mathbf{\Lambda },\mathit{b}}$$

that identifies the regular representation of ${\mathcal{O}}_{l,\mathbf{\Lambda },\mathit{b}}$ with the ${\mathcal{B}}_{l,\mathbf{\Lambda },\mathit{b}}$-module ${\mathcal{M}}_{l,\mathbf{\Lambda },\mathit{b}}$. Therefore, part (ii) is proved.

Since ${\mathcal{B}}_{l,\mathbf{\Lambda },\mathit{b}}$ is a Frobenius algebra, the regular and coregular representations of the algebra ${\mathcal{B}}_{l,\mathbf{\Lambda },\mathit{b}}$ are isomorphic to each other. Parts (iii)–(vi) follow from the general facts about the coregular representations; see, e.g., [1] [Section 3.3] or [13] [Lemma 1.3].

Due to part (iv), it suffices to consider the algebraic multiplicity of every eigenvalue. It is well known that roots of a polynomial continuously depend on its coefficients. Hence, the eigenvalues of $\mathcal{T}\left(x\right)$ continuously depend on $\mathit{b}$. Part (vii) follows from the deformation argument and Theorem 2. □

6. Conclusions

In this paper, we investigated the $\mathfrak{gl}\left(1\right|1)$ Gaudin models that are twisted by a diagonal matrix G and defined on tensor products of polynomial evaluation $\mathfrak{gl}\left(1\right|1\left)\right[t]$-modules. Our results generalize all of the results of [10] to the twisted case. Meanwhile, we gave an explicit description of the algebra of Hamiltonians acting on tensor products of polynomial evaluation $\mathfrak{gl}\left(1\right|1\left)\right[t]$-modules by generators and relations. Moreover, we showed that there exists a bijection between common eigenvectors (up to proportionality) of the algebra of Hamiltonians and monic divisors of an explicit polynomial written in terms of the highest weights and evaluation parameters. In particular, our result implies that each common eigenspace of the algebra of Hamiltonians has dimension one. We also gave dimensions of the generalized eigenspaces. Our results give a confirmed answer to the completeness of the Bethe ansatz in the case of $\mathfrak{gl}\left(1\right|1)$ Gaudin models. We expect our results tp be an essential step towards understanding the more general $\mathfrak{gl}\left(m\right|n)$ Gaudin models.