Abstract
We introduce a new infinite family of higher-order difference operators that commute with the elliptic Ruijsenaars difference operators of type A. These operators are related to Ruijsenaars’ operators through a formula of Wronski type.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
In Ruijsenaars’s study of relativistic quantum integrable systems [11, 12], he introduced a commuting family of linear difference operators in n variables, denoted by \(D_1,\ldots , D_n\), involving sigma functions as coefficients. In this paper, we construct an explicit infinite family of difference operators \(H_0, H_1, H_2,\ldots \) in the commutative algebra \(\mathbb {C}[D_1,\ldots ,D_n]\) which are related to \(D_r\) (\(r=1,\ldots ,n\)) through a formula of Wronski type. This construction is applicable also to difference operators with trigonometric and rational coefficients. In order to deal with the elliptic, trigonometric and rational cases simultaneously, as in [5, 6] we formulate our results in terms of an entire function [z] satisfying the three-term relation (1.1).
We fix a nonzero entire function [z] in one complex variable \(z\in \mathbb {C}\) satisfying the three-term relation of Hirota type
for any \(\alpha ,\beta ,\gamma \in \mathbb {C}\), where \([\alpha \pm \beta ]=[\alpha +\beta ][\alpha -\beta ]\). We remark that a generic solution of this functional equation is given by
where \(\sigma (z;\Omega )\) denotes the Weierstrass sigma function associated with a period lattice \(\Omega =\mathbb {Z}\omega _1\oplus \mathbb {Z}\omega _2\). It is also satisfied by the functions
which are trigonometric and rational degenerations of the generic solution above. It is known that any solution of (1.1) belongs to one of these three categories. Throughout this paper, we denote by \(D_r=D_r^{x}\) (\(r=1,\ldots ,n\)) the Ruijsenaars operators in n variables \(x=(x_1,\ldots ,x_n)\) with parameters \((\delta ,\kappa )\), associated with [z]. They are defined by
where I runs over all subsets of indices of cardinality r and \(T_{x_i}^{\delta }\) stands for the \(\delta \)-shift operator \(x_i\rightarrow x_i+\delta \) in \(x_i\) for each \(i=1,\ldots ,n\). It is proved by Ruijsenaars [11] that these operators \(D_{r}\) commute with each other, namely
on the basis of a certain functional identity for the sigma function. We denote by \(\mathcal {R}=\mathbb {C}[D_1,\ldots ,D_n]\) the commutative algebra generated by \(D_{r}\) (\(r=1,\ldots ,n\)) and refer to it as the commutative algebra of Ruijsenaars operators (of type \(A_{n-1}\)). We also define \(D_{0}=1\), and \(D_{r}=0\) for \(r>n\). Note also that these operators \(D_{r}\) (\(r=1,\ldots ,n\)) in the trigonometric case are the Macdonald q-difference operators expressed additively.
We define an infinite family of difference operators \(H_l=H_l^{x}\) \((l=0,1,2\ldots )\) by
These operators are expressed briefly as
in terms of the difference product \(\Delta (x)=\prod \nolimits _{1\le i<j\le n}{}[x_i-x_j]\), and the \(\delta \)-shifted factorials
for [z]. Also, for each multi-index \(\mu =(\mu _1,\ldots ,\mu _n)\in \mathbb {N}^n\), \(\mathbb {N}=\{0,1,2,\ldots \}\), we define \(|\mu |=\mu _1+\cdots +\mu _n\) and \(T_{x}^{\mu \delta }=T_{x_1}^{\mu _1\delta }\cdots T_{x_n}^{\mu _n\delta }\).
Theorem 1.1
The linear difference operators \(H_l=H_{l}^{x}\) \((l=0,1,2,\ldots )\) defined as above belong to the commutative algebra \(\mathcal {R}=\mathbb {C}[D_1,\ldots ,D_n]\) of Ruijsenaars operators. In particular, one has
We first remark that the family of difference operators \(H_l\) (\(l=0,1,2,\ldots \)) originates from the kernel identities for Ruijsenaars operators ([6]) and the duality transformation formulas for multiple elliptic hypergeometric series ([5, 10]). Let G(z) be a nonzero meromorphic function on \(\mathbb {C}\) such that \(G(z+\delta )=[z]\,G(z)\), and define the kernel function \(\Phi (x;y)\) of Cauchy type by
for \(x=(x_1,\ldots ,x_n)\) and \(y=(y_1,\ldots ,y_n)\). Then, from [5], Theorem 1.3, we have the kernel identity
for the Ruijsenaars operators \(D_r\) (\(r=0,1,\ldots ,r\)). On the other hand, the duality transformation formula for multiple elliptic hypergeometric series implies
as the special case where \(m=n\) and \(a_i=b_i=\kappa \) (\(i=1,\ldots ,n\)) in the notation of [5], Theorem 2.2. This means that
Namely, the kernel function for the Ruijsenaars operators \(D_r\) (\(r=0,1,\ldots ,n\)) simultaneously intertwines the operators \(H_r\) (\(r=0,1,2\ldots \)). In view of this fact, it would be reasonable to expect that the operators \(H_r\) should already belong to the commutative algebra \(\mathbb {C}[D_1,\ldots ,D_n]\) of Ruijsenaars operators. Theorem 1.1 ensures that it is actually the case.
In this paper, we prove Theorem 1.1 as a consequence of the following recurrence formula of Wronski type for \(H_l\) (\(l=0,1,2,\ldots \)).
Theorem 1.2
The difference operators \(H_l\) \((l=0,1,2,\ldots )\) satisfy the following recurrence formula in relation to \(D_r\) \((r=1,\ldots ,n)\) :
Recall that the elementary symmetric functions
in \(\xi =(\xi _1,\ldots ,\xi _n)\), and the complete homogeneous symmetric functions
are related to each other through the Wronski formulaFootnote 1
See [7]. Theorem 1.2 can be thought of as an operator version of this Wronski formula for symmetric functions. A proof of Theorem 1.2 will be given in Sect. 2 by using a functional identity for [z] (Lemma 2.3).
From the recurrence formulas
we see inductively that \(H_l\) belongs to the commutative algebra \(\mathcal {R}=\mathbb {C}[D_1,\ldots ,D_n]\) of Ruijsenaars operators for all \(l=0,1,2,\ldots \). In fact, we have
The relationship between the two families of difference operators \(D_r\) (\(r=01,2,\ldots \)) and \(H_l\) (\(l=0,1,2,\ldots \)) is described in the following.
Proposition 1.3
For each \(l=0,1,2,\ldots \), the difference operator \(H_l\) is expressed in terms of \(D_r\) \((r=0,1,\ldots )\) by the determinant formula
Conversely,
Proposition 1.4
For each \(l=1,2,\ldots \), the difference operator \(H_l\) is expressed explicitly as
in terms of \(D_r\) \((r=0,1,\ldots )\).
We also summarize the kernel identities relevant to the difference operators \(D_r\) and \(H_r\) for the sake of reference.
Theorem 1.5
(1) For two sets of variables \(x=(x_1,\ldots ,x_n)\) and \(y=(y_1,\ldots ,y_n)\), the kernel function \(\Phi (x,y)\) of Cauchy type in (1.10) satisfies the following two types of kernel identities :
(2) For two sets of variables \(x=(x_1,\ldots ,x_m)\) and \(y=(y_1,\ldots ,y_n)\), let
be the kernel function of dual Cauchy type. Under the balancing condition \(m\kappa +n\delta =0\), \(\Psi (x;y)\) satisfies the kernel identity
where \(\widehat{D}^y_r\) denotes the difference operator obtained from \(D^y_r\) by exchanging the parameters \(\delta \) and \(\kappa \).
Propositions 1.3 and 1.4 are consequences of the recurrence formula of Wronski type (see Sect. 3). After a complement on kernel identities for the Ruijsenaars operators (Sect. 4), we finally give some remarks on the trigonometric case in Sect. 5.
Notes: This paper is based on a collaboration of the authors which was completed as master’s thesis [13] of the second author in Japanese. Also, an earlier version of the present paper, written around 2012, has been circulated among some researchers. For these reasons, some of the results in this paper are already cited in several studies [2, 3, 8, 9] with reference to a private communication or to a paper in preparation.
2 Recurrence formula of Wronski type
In this section, we give a proof of Theorem 1.2. Our goal is to establish the recurrence formula of Wronski type between the two sequences of difference operators \(D_r\) (\(r=0,1,2,\ldots \)) and \(H_s\) (\(s=0,1,2,\ldots \)).
Theorem 2.1
The difference operators \(H_l\) \((l=0,1,2\ldots )\) defined by (1.6) satisfy the recurrence formulas
Since \(D_0=1\), by this theorem we see inductively that \(H_l\) belong to \(\mathbb {C}[D_1,\ldots ,D_n]\) for all \(l=0,1,2,\ldots \).
Theorem 2.2
The difference operators \(H_l\) \((l=0,1,2\ldots )\) belong to the commutative algebra \(\mathbb {C}[D_1,\ldots ,D_n]\) of Ruijsenaars operators. In particular, one has
Proof of Theorem 2.1
We express the difference operators \(D_r\) as
where we define \(\epsilon _I=\sum \nolimits _{i\in I}{}{\epsilon _i}\) by using the unit vectors \(\epsilon _1,\ldots ,\epsilon _n\) of \(\mathbb {N}^n\). Similarly, we express \(H_l\) as
When \(r+s=l\), we compute
where \(\text{ supp }(\lambda )=\{i\in \{1,\ldots ,n\}\ |\ \lambda _i>0\}\). Hence, the recurrence formula (2.1) is equivalent to
saying that
for any \(\lambda \in \mathbb {N}^n\) with \(|\lambda |>0\). We now make the expression \(A_I(x)H_{\lambda -\epsilon _I}(x+\epsilon _I\delta )\) explicit. Setting \(L=\text{ supp }(\lambda )\), we have
Noting that
we can compute \(A_I(x)H_{\lambda -\epsilon _I}(x+\epsilon _I\delta )\) as follows:
Hence, (2.7) is equivalent to
for any \(L\ne \phi \) and \(\lambda \in \mathbb {N}^n\) with \(\text{ supp }(\lambda )=L\). Setting \(\lambda =\epsilon _L+\nu \), we rewrite this in the form
for any \(\nu \in \mathbb {N}^n\) with \(\text{ supp }(\nu )\subseteq L\). Since this formula contains only those variables \(x_i\) with \(i\in L\), we have only to consider the case where \(L=\{1,\ldots ,n\}\) (\(n\ge 1\)):
for any \(\nu \in \mathbb {N}^n\). This identity follows from the following functional identity by the change of variables
\(\square \)
Lemma 2.3
Given two sets of variables \(z=(z_1,\ldots ,z_n)\), \(w=(w_1,\ldots ,w_n)\) and a parameter a, the following identity holds as a meromorphic function in \((z_1,\ldots ,z_n,w_1,\ldots ,w_n)\) for \(n\ge 1\) \(:\)
where \(|z|=\sum _{i=1}^nz_i\), \(|w|=\sum _{k=1}^nw_k\).
Proof
We give a proof of the functional identity (2.15) for the case where \([z]=\sigma (z;\Omega )\) is the Weierstrass sigma function associated with a period lattice \(\Omega \). By the classification of [z], it is not difficult to derive (2.15) for any [z] in this class from the case of \(\sigma (z;\Omega )\), by the limiting procedures from \(\sigma (z;\Omega )\) to \(\sin (\pi z/\omega )\) and z, and by the invariance of (2.15) under the transformation \([z] \rightarrow e^{cz^2}[z]\).
Identity (2.15) for \([z]=\sigma (z;\Omega )\) can be proved by the induction on n. Since it holds trivially for \(n=1\), we assume \(n\ge 2\). We regard the left-hand side of (2.15) as a meromorphic function of \(w_n\) and denote it by \(F(w_n)\) assuming that the other variables are generic. Note first that \(F(w_n)\) is an elliptic function possibly with simple poles at \(w_n\equiv z_1-a,\ldots , z_n-a\) and \(w_n\equiv |z|-|w'|\) modulo the period lattice \(\Omega \), where \(w'=(w_1,\ldots ,w_{n-1})\). We first compute the residue of \(F(w_n)\) at \(w_n=z_n-a\). Nontrivial residues possibly arise from the terms corresponding to I containing n; we parametrize such I’s as \(I=J\cup \{n\}\) with \(J\subseteq \{1,\ldots ,n-1\}\). Then, we have
by the induction hypothesis (\(z'=(z_1,\ldots ,z_{n-1})\)). Since \(F(w_n)\) is symmetric with respect to \((z_1,\ldots ,z_n)\), we see that \(w_n\equiv z_1-a,\ldots ,z_n-a\) are all removable poles of \(F(w_n)\). Hence, \(F(w_n)\) has at most one simple pole in each fundamental parallelogram, which is impossible unless \(F(w_n)\) is a constant function since it is an elliptic function. We next look at the value of \(F(w_n)\) at \(w_n=z_n\). It is computed as
again by the induction hypothesis. This implies that \(F(w_n)\) is identically zero as a meromorphic function of \(w_n\). \(\square \)
Remark 2.4
Lemma 2.3 can be proved in a different way if we make use of the argument of [5]. Recall that in (1.14) of Sect. 1, [5], the following identity is derived from the determinant formula of Frobenius:
where D(z; w) is the Frobenius determinant
and \(E(T_z;u)={\textstyle \prod _{i=1}^{n}}(1+uT_{z_i}^{\delta })\). This means that
By setting \(u=-1\), we obtain
Note that in the limit \(\lambda \rightarrow 0\), each entry of the matrix of the left-hand side has a finite limit
Hence, the left-hand side converges to zero as \(\lambda \rightarrow 0\). This implies that
By replacing each \(z_i\) with \(-z_i\), and \(\delta \) with a, we obtain Lemma 2.3.
3 Explicit relations between the two commuting families
Setting
we rewrite the recurrence formula of Theorem 2.1 as
In the matrix form, this means that
Hence, by Cramer’s formula we obtain
Namely, we have
By symmetry, we also have
The recurrence formula (2.1) can also be written as
Applying this formula repeatedly, we obtain
4 Kernel identities
We recall from [5] the duality transformation formula for multiple elliptic hypergeometric series (of type A): Under the balancing condition \(a_1+\cdots +a_m=b_1+\cdots +b_n\),
where \(x=(x_1,\ldots ,x_m)\) and \(y=(y_1,\ldots ,y_n)\). As we already remarked, when \(m=n\) and \(a_i=b_i=\kappa \) (\(i=1,\ldots ,n\)), this implies
Let \(a_1=\cdots =a_m=\kappa \), \(b_1=\cdots =b_n=-\delta \). Then, this transformation formula implies that under the condition \(m\kappa +n\delta =0\),
By replacing \(y_k\) by \(-y_k\) for \(k=1,\ldots ,n\), we obtain
This means that the dual Cauchy kernel
satisfies
under the condition \(m\kappa +n\delta =0\).
5 The trigonometric cases
In this section, we consider the trigonometric cases where \([x]=e(x/2)-e(-x/2)\) in the notation \(e(u)=\exp (2\pi \sqrt{-1}u)\) of the exponential function. Instead of the parameter \(\delta \), \(\kappa \in \mathbb {C}\), we use the multiplicative parameters \(q=e(\delta )\) and \(t=e(\kappa )\) assuming that \(\mathrm{Im}(\delta )>0\) so that \(|q|<1\). Note that when \(z=e(x)\), we have \([x]=z^{\frac{1}{2}}-z^{-\frac{1}{2}}=-z^{-\frac{1}{2}}(1-z)\), and hence,
in the standard notation \((z;q)_k=(1-z)(1-qz)\cdots (1-q^{k-1}z)\) of q-shifted factorials.
We denote by \(z=(z_1,\ldots ,z_n)\) the multiplicative variables defined by \(z_i=e(x_i)\) (\(i=1,\ldots ,n\)) corresponding to \(x=(x_1,\ldots ,x_n)\). For these z variables, we denote by \(T_{q,z_i}\) the q-shift operator with respect to \(z_i\) (\(i=1,\ldots ,n\)) and set \(T_{q,z}^{\mu }=T_{q,z_1}^{\mu _1}\cdots T_{q,z_n}^{\mu _n}\) for each multi-index \(\mu =(\mu _1,\ldots ,\mu _n)\in \mathbb {N}^n\). In this multiplicative notation, it is convenient to introduce the q-difference operators \(\mathcal {D}_r^{z}\) and \(\mathcal {H}_l^z\) normalized so that
These q-difference operators are given explicitly by
and
The recurrence relation (2.1) of Wronski type is then rewritten as follows:
It is known by [7] that the commuting family of q-difference operators \(\mathcal {D}_r^{z}\) (\(r=0,1,\ldots ,n\)) act on the ring \(\mathbb {C}[z]^{\mathfrak {S}_n}=\mathbb {C}[z_1,\ldots ,z_n]^{\mathfrak {S}_n}\) of symmetric polynomials in \(z=(z_1,\ldots ,z_n)\) and that they are simultaneously diagonalized by the (monic) Macdonald polynomials \(P_\lambda (z)=P_\lambda (z|q,t)\) indexed by partitions \(\lambda =(\lambda _1,\ldots ,\lambda _n)\) with \(l(\lambda )\le n\):
where \(e_r(\xi )\) stands for the elementary symmetric polynomial of degree r for each \(r=0,1,\ldots ,n\), and \(t^\delta q^{\lambda }=(t^{n-1}q^{\lambda _1},t^{n-2}q^{\lambda _2},\ldots ,q^{\lambda _n})\). In terms of the generating function
formula (5.6) is equivalent to
Since \(\mathcal {H}_l^z\in \mathbb {C}[\mathcal {D}_1^z,\ldots ,\mathcal {D}_n^z]\), the q-difference operators \(\mathcal {H}_l^z\) \((l=0,1,2\ldots )\) satisfy
for some symmetric polynomials \(g_l(\xi )\in \mathbb {C}[\xi ]^{\mathfrak {S}_n}\). By the Wronski-type formula (5.5), these polynomials are determined by the recurrence relation
In view of
let us introduce the generating function \(G(\xi ;u)=\sum \limits _{l=0}^{\infty }\,u^l\,g_l(\xi )\). Then, the recurrence formula above is equivalent to the functional equation
namely
Hence, we have
where \((u;q)_\infty =\prod _{i=0}^{\infty }(1-q^i u)\). This means that
We introduce the generation function
for our q-difference operators \(\mathcal {H}_l^z\) (\(l=0,1,2,\ldots \)). Then, the argument above implies that
for any partition \(\lambda =(\lambda _1,\ldots ,\lambda _n)\) with \(l(\lambda )\le n\). Note also that the recurrence formula of Wronski type is equivalent to
Finally, we give comments on the kernel identities for the trigonometric case. Consider two sets of variables \(z=(z_1,\ldots ,z_m)\) and \(w=(w_1,\ldots .w_n)\), assuming that \(m\ge n\). The Cauchy-type kernel for this case is given by
Then, we have the kernel identities
By the kernel function of dual Cauchy type
the two families of q-difference operators are exchanged as follows:
where \(\widehat{\mathcal {D}}^w(u)=\mathcal {D}^{(w|t,q)}(u)\) denotes the q-difference operator obtained from \(\mathcal {D}^w(u)=\mathcal {D}^{(w|q,t)}(u)\) by exchanging q and t.
The three kernel identities (DD), (HH) and (HD) are equivalent to certain special cases of Kajihara’s Euler transformation formula [4]: For two sets of variables \(z=(z_1,\ldots ,z_m)\), \(w=(w_1,\ldots ,w_n)\) and parameters \(a=(a_1,\ldots ,a_m)\), \(b=(b_1,\ldots ,b_n)\),
where \(\alpha =a_1\cdots a_m\) and \(\beta =b_1\cdots b_n\). In fact, one can verify directly that these three kernel identities are equivalent to the following special cases of (5.24), respectively:
References
Aitken, A.C.: Determinants and Matrices. Oliver and Boyd, Edinburgh (1939)
Borodin, A., Corwin, I., Remenik, D.: A classical limit of Noumi’s \(q\)-integral operator. SIGMA Symmetry Integr. Geom. Methods Appl. 11, Paper 098, 7 (2015)
Feigin, B., Hashizume, K., Hoshino, A., Shiraishi, J., Yanagida, K.: A commutative algebra on degenerate \(CP ^1\) and Macdonald polynomials. J. Math. Phys. 50, 095215, 42 pages (2009)
Kajihara, Y.: Euler transformation formula for multiple basic hypergeometric series of type \(A\) and some applications. Adv. Math. 187, 53–97 (2004)
Kajihara, Y., Noumi, M.: Multiple elliptic hypergeometric series. An approach from the Cauchy determinant. Indag. Math. N.S. 14, 395–421 (2003)
Komori, Y., Noumi, M., Shiraishi, J.: Kernel functions of difference operators of Ruijsenaars type and their applications. SIGMA 5, 054, 40 pages (2009)
Macdonald, I.G.: Symmetric Functions and Hall Polynomials, 2nd edn. Oxford University Press, Oxford (1995)
Masuda, Y.: Kernel identities for van Diejen’s \(q\)-difference operators and transformation formulas for multiple basic hypergeometric series. Ramanujan J. 32, 281–314 (2013)
Noumi, M., Shiraishi, J.: A direct approach to the bispectral problem for the Ruijsenaars- Macdonald \(q\)-difference operators. arXiv:1206.5364, 44 pages
Rosengren, H.: New transformations for elliptic hypergeometric series on the root system \(A_n\). Ramanujan J. 12, 155–166 (2006)
Ruijsenaars, S.N.M.: Complete integrability of relativistic Calogero-Moser systems and elliptic function identities. Commun. Math. Phys. 110, 191–213 (1987)
Ruijsenaars, S.N.M.: Hilbert-Schmidt operators vs. Integrable systems of Elliptic Calogero-Moser type I. The eigenfunction identities, Commun. Math. Phys. 286, 629–657 (2009); II. The \(A_{N-1}\) case: First steps, Commun. Math. Phys. 286, 659–680 (2009)
Sano, A.: A commuting family of linear difference operators involving elliptic functions as coefficients (in Japanese), Master’s Thesis, Department of Mathematics, Kobe University, 2008
Turnbull, H.W.: Theory of Equations. Oliver and Boyd, Edinburgh (1939)
Umeda, T.: Application of Koszul complex to Wronski relations for \(U(\mathfrak{gl}_n)\). Comment. Math. Helv. 78, 663–680 (2003)
Acknowledgements
M.N. is grateful to the Knut and Alice Wallenberg Foundation (Grant Number KAW 2019.0525) for funding his guest professorship at KTH.
Funding
Open access funding provided by Royal Institute of Technology.
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Masatoshi Noumi: On leave from: Department of Mathematics, Kobe University, Rokko, Kobe 657-8501, Japan.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Noumi, M., Sano, A. An infinite family of higher-order difference operators that commute with Ruijsenaars operators of type A. Lett Math Phys 111, 91 (2021). https://doi.org/10.1007/s11005-021-01435-9
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1007/s11005-021-01435-9