## Exact solutions by integrals of the non-stationary elliptic Calogero-Sutherland equation

Farrokh Atai 1,* and Edwin Langmann 2, †

1 Department of Mathematics, Kobe University, Rokko, Kobe 657-8501, Japan 2 Department of Physics, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden

August 2, 2019

## Abstract

We use generalized kernel functions to construct explicit solutions by integrals of the non-stationary Schr¨dinger equation for the Hamiltonian of the elliptic Calogeroo Sutherland model (also known as elliptic Knizhnik-Zamolodchikov-Bernard equation). Our solutions provide integral represenations of elliptic generalizations of the Jack polyomials.

MSC class: 33E30, 32A26, 81Q05, 16R60

Keywords: Elliptic integrable systems; quantum Calogero-Moser-Sutherland systems; elliptic Knizhnik-Zamolodchikov-Bernard equation; kernel functions; affine analogue of Jack polynomials; integral representations of special function.

## 1 Introduction

In this paper we construct explicit solutions of the partial differential equation

$$\begin{array} { c } \cdot & \cdot & \cdot & \cdot & \cdot \\ & \left ( \kappa \frac { i } { \pi } \frac { \partial } { \partial \tau } - \sum _ { j = 1 } ^ { n } \frac { \partial ^ { 2 } } { \partial x _ { j } ^ { 2 } } + \sum _ { j < k } ^ { n } 2 g ( g - 1 ) \wp ( x _ { j } - x _ { k } | \pi, \pi \tau ) \right ) \psi ( x _ { 1 }, \dots, x _ { n } ; \tau ) = 0 \\ \end{array} \quad \begin{array} { c } \cdot & \cdot & \cdot & \cdot & \cdot & \cdot \\ & \left ( \kappa \frac { i } { \pi } \frac { \partial } { \partial \tau } - \sum _ { j = 1 } ^ { n } \frac { \partial ^ { 2 } } { \partial x _ { j } ^ { 2 } } + \sum _ { j < k } ^ { n } 2 g ( g - 1 ) \wp ( x _ { j } - x _ { k } | \pi, \pi \tau ) \right ) \psi ( x _ { 1 }, \dots, x _ { n } ; \tau ) = 0 \\ \end{array} \quad \begin{array} { c } \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot \\ & \left ( \kappa \frac { i } { \pi } \frac { \partial } { \partial \tau } - \sum _ { j = 1 } ^ { n } \frac { \partial } { \partial x _ { j } ^ { 2 } } + \sum _ { j < k } ^ { n } - \sum _ { j = 1 } ^ { n } \frac { \partial } { \partial x _ { j } ^ { 2 } } + \sum _ { j < k } ^ { n } - \sum _ { j = 1 } ^ { n } \frac { \partial } { \partial x _ { j } ^ { 2 } } \end{array}$$

where ℘ x π, πτ ( | ) is the usual Weierstrass function with periods (2 π, 2 πτ ), x j ∈ -[ π, π ] are variables on a circle ( j = 1 , . . . , n ), g &gt; 1 / 2, and κ = kg with integer k ≥ 1 (the definition of the function ℘ can be found in Appendix A; our notation otherwise is explained at the end of this section). We refer to this as non-stationary elliptic Calogero-Sutherland (eCS) equation since, for -i τ &gt; 0 and i κ &gt; 0, it can be interpreted as non-stationary Schr¨dinger equation o for the Hamiltonian defining the eCS model [1, 2, 3] 1 with the parameter i κ corresponding to a time scale. 2

∗ Electronic address: farrokh@math.kobe-u.ac.jp

† Electronic address: langmann@kth.se

1 Calogero-Sutherland model is short for A n -1 quantum Calogero-Moser-Sutherland model .

2 In physics units /planckover2pi1 = 2 m = 1, the time variable in this Schr¨dinger equation is o t = -πτ/κ &gt; 0.

Our solutions have the form

/negationslash

$$\psi ( x _ { 1 }, \dots, x _ { n } ; \tau ) & = C ( \tau ) \left ( \prod _ { j \neq k } ^ { n } \theta \left ( z _ { j } / z _ { k } ; p \right ) \right ) ^ { g / 2 } \mathcal { P } ( z _ { 1 }, \dots, z _ { n } ; p ) \quad ( z _ { j } = e ^ { i x _ { j } }, \ p = e ^ { 2 \pi i \tau } ), \\ \text{with}$$

with the usual multiplicative theta function, and P ( z , . . . , z 1 n ; p ) are symmetric functions in the variables ( z , . . . , z 1 n ) that reduce to symmetric Laurent polynomials in the limit p → 0. 3 Our main result are explicit integrals representing the latter functions.

$$\theta ( z ; p ) \equiv ( 1 - z ) \prod _ { \ell = 1 } ^ { \infty } ( 1 - p ^ { \ell } z ) ( 1 - p ^ { \ell } / z ) \quad ( z \in \mathbb { C }, \ | p | < 1 ) \\ \text{ usual multiplicative theta function, and } \mathcal { P } ( z _ { 1 }, \dots, z _ { n } ; p ) \text{ are symmetric functions in the}$$

In particular, for κ = g , we construct such solutions ψ λ ( x , . . . , x 1 n ; τ ) of (1) labeled by integer vectors λ = ( λ , . . . , λ 1 n ) ∈ Z n , and the corresponding integrals we thus obtain are

/negationslash

/negationslash where the integration contours are circles with radii ρ a &gt; 1 such that | p ρ | a -1 &lt; ρ a &lt; ρ a -1 for all a , with arbitrary ρ n -1 in the range 1 &lt; ρ n -1 &lt; 1 / p | | . We show that these integrals define analytic functions in the region | p ρ | n -1 &lt; z | j | &lt; ρ n -1 and | p | &lt; 1. Moreover, when λ 1 ≥ λ 2 ≥ · · · ≥ λ n , then the limiting case p → 0 of (3) reduces (essentially 4 ) to an integral representation of the Jack polynomials [4] obtained by Awata, Matsuo, Odake &amp; Shiraishi in [5] in generalization of a result by Mimachi &amp; Yamada [6], up to explicitly known nonzero constants (see Section 3). For integer vectors λ ∈ Z n that do not satisfy the condition λ 1 ≥ λ 2 · · · ≥ λ n , the functions P λ ( z , . . . , z 1 n ; p ) vanish in the limit p → 0, but they are non-trivial for non-zero p .

$$\text{integer vectors} \, \lambda = ( \lambda _ { 1 }, \dots, \lambda _ { n } ) \in \mathbb { Z } ^ { \prime \prime }, \, \text{and the corresponding integrals we thus obtain are} \\ \mathcal { P } _ { \lambda } ( z _ { 1 }, \dots, z _ { n } ; p ) & = \left ( \prod _ { j = 1 } ^ { n } z _ { j } ^ { \lambda _ { n } } \right ) \left ( \prod _ { a = 1 } ^ { n - 1 } \prod _ { j = 1 } ^ { a } \oint _ { | \xi _ { a j } | = \rho _ { a } } \frac { d \xi _ { a j } } { 2 \pi i \xi _ { a j } ^ { \lambda _ { a - 1 } } } \xi _ { a j } ^ { \lambda _ { a - 1 } - 1 } \right ) \\ & \quad \times \left ( \frac { \prod _ { j \neq k } ^ { n - 1 } \theta ( \xi _ { n - 1, j } / \xi _ { n - 1, k } ; p ) ^ { g } } { \prod _ { j = 1 } ^ { n - 1 } \prod _ { k = 1 } ^ { n - 1 } \theta ( z _ { j } / \xi _ { n - 1, k } ; p ) ^ { g } } \right ) \left ( \prod _ { a = 1 } ^ { n - 2 } \prod _ { j \neq k } ^ { a } \theta ( \xi _ { a j } / \xi _ { a k } ; p ) ^ { g } } \right ) \quad ( 3 ) \\ \text{where the integration contours are circles with radii} \, \rho _ { a } > 1 \, \text{such that} \, | p | \rho _ { a - 1 } < \rho _ { a } < \rho _ { a - 1 } \\ \text{for all} \, a, \, \text{with arbitrary} \, \rho _ { n - 1 } \, \text{in the range} \, 1 < \rho _ { n - 1 } < 1 / | p |. \, \text{We show that these integrals}$$

Our solution for κ = g is complete in the sense that we obtain elliptic generalizations of all Jack polynomials. Our general results are solutions for κ = kg for integer k ≥ 1, but for k ≥ 2 we only get elliptic generalizations for some Jack polynomials; a precise statement of our results for all integers k ≥ 1 is given in Theorem 3.1.

Remark 1.1 . The integral representations of the Jack polynomials obtained from our result in the limit p → 0 differ from the results in [5] for non-integer g in the definition of the integration contours. More specifically, our integration contours are always non-intersecting concentric circles (see Fig. 2(a) in [5]), whereas, for non-integer g , the ones in [5] are concentric circles that are pinched so as to intersect at one point (see Fig. 2(b) in [5]). To convince critical readers that our simpler integration contours work even for non-integer g , we give in Appendix B a concise proof of the result in [5] explaining this technical point. 5

3 For simplicity, we have chosen not to specify the normalization function C τ ( ) in the introduction. However, this function is important to understand the limit p → 0; see Section 3 for details.

5 We are grateful to J. Shiraishi for valuable comments prompting us to write Appendix B.

4 There is a difference in the definition of the integration contours; see Remark 1.1 for details. Moreover, if λ n &lt; 0, then we get a Jack polynomial up to a factor ( z 1 · · · z n ) λ n ; see Remark 2.1.

We are motivated by Shiraishi's non-stationary eCS functions [7, Section 7] defined by explicit formal power series and which, as Shiraishi conjectured, provide exact eigenfunctions of the non-stationary eCS equation, for arbitrary κ ∈ C \{ } 0 and g ∈ C . As discussed in Section 5, we hope that our results can provide a useful test of Shiraishi's conjecture.

The special case n = 2 of our results is equivalent to an integral representation of elliptic generalizations of the Gegenbauer polynomials obtained in [8]; see [9, Section 3.2] in the special case where ˜ g ν = g for all ν = 0 1 2 3 and [10] for similar such results. , , , The key to our generalization of this result are generalized kernel function identities for the eCS model obtained in [11].

The non-stationary eCS equation is also known as elliptic Knizhnik-ZamolodchikovBernard (KZB) equation [12, 13] in the literature, 6 and integral representations of its solutions were obtained by Etingof &amp; Kirillov Jr. [15] and Felder &amp; Varchenko [16] using a representation theoretic approach; see also [17, 18]. These previous results are for g ∈ Z ≥ 1 and κ ∈ C \ { 0 } , whereas our results hold true for arbitrary values g &gt; 1 / 2 and κ ∈ g Z ≥ 1 . Moreover, our method is different in that it is based on generalized kernel functions rather than representation theory, and our results are more explicit in the sense that they are expressed in terms of more elementary functions and integrations.

Plan: In the rest of this paper we review the integral representations of the Jack polynomials (Section 2), give a mathematically precise formulation of our results (Section 3), prove our results (Section 4), and conclude with some remarks to put our results in perspective (Section 5). For the convenience of the reader we include two appendices with definitions and properties of standard special functions that we use (Appendix A) and a short proof of the integral representations of the Jack polynomials (Appendix B). A third appendix contains some technical details in our proof in Section 4 (Appendix C).

Notation: Throughout this paper, n ∈ Z ≥ 1 , τ ∈ C with imaginary part /Ifractur ( τ ) &gt; 0, and i ≡ + √ -1. We write ∑ n j&lt;k for the sum over all j, k = 1 2 , , . . . , n such that j &lt; k , etc. For ξ a complex variable, the symbol ∮ | ξ | = ρ dξ ( · · · ) is used for the integral over the counterclockwise oriented circle | ξ | = ρ of radius ρ &gt; 0. We often follow the notation and conventions in Macdonald's book [19]. In particular, we write x short for ( x , . . . , x 1 n ) in the following, and similarly for z , etc. If z = ( z , . . . , z 1 n ) ∈ C n then 1 /z is short for (1 /z , . . . , 1 1 /z n ). For L ∈ Z ≥ 1 , r = ( r , . . . , r 1 L ) ∈ Z L and s = ( s , . . . , s 1 L ) ∈ Z L ≥ 1 , we sometimes write

$$( r _ { 1 } ^ { s _ { 1 } }, \dots, r _ { L } ^ { s _ { L } } ) \text{ short for } \underbrace { ( r _ { 1 }, \dots, r _ { 1 }, \dots, r _ { L }, \dots, r _ { L } ) } _ { s _ { 1 } \text{ times} } ), \\ \text{mple, (3^{\,2}, ( - 1 ) ^ { 3 }, 2 ^ { 1 } ) \text{ is short for (3,3, -1, -1, -1, 2)}.$$

for example, (3 2 , ( -1) 3 , 2 ) is short for (3 1 , 3 , -1 , -1 , -1 2). ,

Acknowledgements: We thank M. Halln¨ as, M. Noumi, H. Rosengren, B. Shapiro, and J. Shiraishi for helpful discussions. We are grateful to M. Noumi and J. Shiraishi for their interest, encouragement, and advise. We gratefully acknowledge partial financial support

6 To be more precise: it is a special case of the KZB equations corresponding to g = sl n and genus 1; see e.g. [14].

by the Stiftelse Olle Engkvist Byggm¨stare (contract 184-0573). E.L. acknowledges support a by VR Grant No. 2016-05167. The work of F.A. was carried out as a JSPS International Research Fellow.

## 2 Jack polynomials (review)

For the convenience of non-expert readers, we shortly review a definition and some properties of the Jack polynomials (Section 2.1). We also present the integral representations of the Jack polynomials obtained in [5], elaborating on the integration contours issue described in Remark 1.1 (Section 2.2 and Appendix B).

The rest of this section can be skipped without loss of continuity.

## 2.1 Definition and properties

For n ∈ Z ≥ 1 , we denote by P n the set of all partitions of length ≤ n , i.e., λ ∈ P n if λ = ( λ , . . . , λ 1 n ) ∈ Z n satisfies λ 1 ≥ · · · ≥ λ n ≥ 0. We recall that, for λ ∈ P n and z = ( z , . . . , z 1 n ) ∈ C n , the monomial symmetric polynomials are defined as

$$m _ { \lambda } ( z ) \equiv \sum _ { a } z _ { 1 } ^ { a _ { 1 } } \cdots z _ { n } ^ { a _ { n } } \\ \intertext { n o r m u t a t i o n } \, a - l a. \quad a \ \.$$

with the sum over all distinct permutations a = ( a , . . . , a 1 n ) of λ = ( λ , . . . , λ 1 n ).

Let Λ n ≡ C [ z , . . . , z 1 n ] S n be the ring of symmetric polynomials u z , ( 1 · · · , z n ) in the variables z = ( z , . . . , z 1 n ) ∈ C n with complex coefficients. 7 The monomial symmetric polynomials m z λ ( ) obviously constitute a basis in Λ n labeled by partitions λ ∈ P n . The Jack polynomials P (1 /g ) λ,n ( z ) constitute another such basis depending on a parameter 1 /g &gt; 0. They can be defined by a Gram-Schmidt orthogonalization for the following scalar product,

/negationslash

$$\langle & f, g \rangle _ { n } ^ { \prime } \equiv \frac { 1 } { n! } \oint _ { | z _ { 1 } | = 1 } \frac { d z _ { 1 } } { 2 \pi i z _ { 1 } } \cdots \oint _ { | z _ { n } | = 1 } \frac { d z _ { n } } { 2 \pi i z _ { n } } \left ( \prod _ { j \neq k } ^ { n } ( 1 - z _ { j } / z _ { k } ) \right ) ^ { g } f ( z ) \overline { g ( z ) } \quad ( f, g \in \Lambda _ { n } ) \quad ( 5 ) \\ & \text{the how random} \, \text{comlaw communication} \, \text{mixing the normal} \, \text{matrix normal} \, \text{and}$$

(the bar denotes complex conjugation) using the monomial symmetric polynomials and dominance partial ordering ,

$$\sum _ { - \. } \mu \leq \lambda \Leftrightarrow \mu _ { 1 } + \dots + \mu _ { j } \leq \lambda _ { 1 } + \dots + \lambda _ { j } \quad \forall j = 1, 2, \dots, n \quad ( \mu, \lambda \in \mathbf P _ { n } ). \\ - \. \quad \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \ dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots$$

To be precise, it was proved by Macdonald in [20] that the following two conditions define the Jack polynomials P (1 /g ) λ ∈ Λ n uniquely,

$$( \cdot )$$

7 We work with complex numbers since we are motivated by quantum mechanics.

$$\L y m o m a s { \tau } _ { \lambda } & \quad \subset \lambda _ { n } \text{ unique}, \\ ( a ) & \quad P _ { \lambda } ^ { ( 1 / g ) } = m _ { \lambda } + \sum _ { \mu < \lambda } v _ { \lambda \mu } m _ { \mu } \text{ for some } v _ { \lambda \mu } \in \mathbb { R }, \\ ( b ) & \quad \left \langle P _ { \lambda } ^ { ( 1 / g ) }, m _ { \mu } \right \rangle _ { n } ^ { \prime } = 0 \text{ for all } \mu < \lambda. \\ \frac { \text{ with complex numbers since we are motivated by quantum mechanics.}$$

(In the following we write P (1 /g ) λ,n rather than P (1 /g ) λ to avoid possible confusion.)

We quote three well-known properties of the Jack polynomials (they are all stated and proved in Macdonald's book [19]). First, the following orthogonality relation with respect to the scalar product in (5), with the quadratic norms

$$\left \langle P _ { \lambda, n } ^ { ( 1 / g ) }, P _ { \mu, n } ^ { ( 1 / g ) } \right \rangle _ { n } ^ { \prime } = \delta _ { \lambda \mu } \mathcal { N } _ { \lambda, n } ( g ) \\ \text{norms}$$

$$\mathcal { N } _ { \lambda, n } ( g ) \equiv \prod _ { j < k } ^ { n } \frac { \Gamma ( \lambda _ { j } - \lambda _ { k } + g ( k - j + 1 ) ) \Gamma ( \lambda _ { j } - \lambda _ { k } + g ( k - j + 1 ) + 1 ) } { \Gamma ( \lambda _ { j } - \lambda _ { k } + g ( k - j ) ) \Gamma ( \lambda _ { j } - \lambda _ { k } + g ( k - j ) + 1 ) } \quad \quad ( 7 ) \\ \hat { \lambda } \quad. \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \end{array}$$

( δ λµ is the Kronecker symbol equal to 1 for λ = µ and 0 otherwise, and Γ is the Euler gamma function). Second, a generating function due to Stanley [21], 8

$$& \prod _ { j = 1 } ^ { n } \prod _ { k = 1 } ^ { m } \frac { 1 } { ( 1 - z _ { j } / \xi _ { k } ) ^ { g } } = \sum _ { \lambda \in P _ { m } } \frac { 1 } { b _ { \lambda } ( g ) } P _ { \lambda, n } ^ { ( 1 / g ) } ( z ) P _ { \lambda, m } ^ { ( 1 / g ) } ( 1 / \xi ) \quad ( m \leq n ) \\. & \quad \quad \, _ { - \, \ e m } \,. \,.$$

for n, m ∈ Z ≥ 1 , with

$$b _ { \lambda } ( g ) \equiv \prod _ { ( j, k ) \in \lambda } \frac { \lambda _ { j } - k + g ( \lambda _ { k } ^ { \prime } - j + 1 ) } { \lambda _ { j } - k + 1 + g ( \lambda _ { k } ^ { \prime } - j ) } \quad ( \lambda \in \mathbf P _ { m } ) \\.$$

where ∏ ( j,k ) ∈ λ here is short for ∏ m j =1 ∏ λ j k =1 ( λ ′ is the partition conjugate to λ ). Third, the following Pieri relation of the Jack polynomials,

$$( z _ { 1 } \cdots z _ { n } ) ^ { r } P _ { \lambda, n } ^ { ( 1 / g ) } ( z ) = P _ { \lambda + ( r ^ { n } ), n } ^ { ( 1 / g ) } ( z ) \quad ( \lambda \in \mathbf P _ { n }, r \in \mathbb { Z } _ { \geq 1 } )$$

with ( r n ) short for ( r, . . . , r ) ( n times).

Remark 2.1 . Using (10), one can naturally extend the definition of the Jack polynomials P (1 /g ) λ,n ( z ) to all ordered integer vectors λ , i.e., to all λ = ( λ , . . . , λ 1 n ) ∈ Z n satisfying the condition

$$\lambda _ { 1 } \geq \dots \geq \lambda _ { n } \\. \,. \,. \,. \,.$$

(that is, one can drop the restriction that λ n is non-negative). Thus, for λ n &lt; 0, we define 9

$$P _ { \lambda, n } ^ { ( 1 / g ) } ( z ) \equiv ( z _ { 1 } \cdots z _ { n } ) ^ { \lambda _ { n } } P _ { \lambda - ( \lambda _ { n } ^ { n } ), n } ^ { ( 1 / g ) } ( z ).$$

In the following we use this generalized definition of the Jack polynomials even though, for λ n &lt; 0, they are Laurent polynomials. Note that the orthogonality relations in (6)-(7) hold true also for these generalized Jack polynomials, and N λ,n ( g ) = N λ +( r n ) ,n ( g ) for all r ∈ Z .

8 Note that, below, we make no distinction between λ ∈ P m and λ = ( λ , . . . , λ 1 m , 0 n -m ) ∈ P n , i.e., we identify partitions that differ only by a string of zeros.

9 Note that λ -( λ n n ) always is a partition under the stated conditions.

## 2.2 Integral representations

$$\text{For } L \in \mathbb { Z } _ { \geq 1 }, \, r & = ( r _ { 1 }, \dots, r _ { L } ) \in \mathbb { Z } ^ { L } \text{ and } s = ( s _ { 1 }, \dots, s _ { L } ) \in \mathbb { Z } ^ { L } _ { \geq 1 }, \, \text{let} \\ \lambda ^ { ( a ) } & \equiv ( r _ { 1 } ^ { s _ { 1 } }, \dots, r _ { a } ^ { s _ { a } } ), \quad N _ { a } \equiv s _ { 1 } + \dots + s _ { a } \quad ( a = 1, \dots, L ), \quad n \equiv N _ { L }, \quad \quad ( 1 1 ) \$$

and

/negationslash

/negationslash with integration contour radii ρ a &gt; 0 restricted as follows,

$$\mathcal { P } _ { r, s, L } ( z _ { 1 }, \dots, z _ { n } ) \equiv & \left ( \prod _ { j = 1 } ^ { n } z _ { j } ^ { r _ { L } } \right ) \left ( \prod _ { a = 1 } ^ { L - 1 } \prod _ { j = 1 } ^ { N _ { a } } \oint _ { | \xi _ { a j } | = \rho _ { a } } \frac { d \xi _ { a j } } { 2 \pi \mathbf 1 \xi _ { a j } } \xi _ { a j } ^ { r _ { a } - r _ { a + 1 } } \right ) \\ \times & \left ( \frac { \prod _ { j \neq k } ^ { N _ { L - 1 } } ( 1 - \xi _ { N _ { L - 1 }, j } / \xi _ { N _ { L - 1 }, k } ) ^ { g } } { \prod _ { j = 1 } ^ { n } \prod _ { k = 1 } ^ { N _ { L - 1 } } ( 1 - z _ { j } / \xi _ { N _ { L - 1 }, k } ) ^ { g } } \right ) \left ( \prod _ { a = 1 } ^ { L - 2 } \frac { \prod _ { j \neq k } ^ { N _ { a } } ( 1 - \xi _ { a j } / \xi _ { a k } ) ^ { g } } { \prod _ { j = 1 } ^ { N _ { a + 1 } } \prod _ { k = 1 } ^ { N _ { a } } ( 1 - \xi _ { a + 1, j } / \xi _ { a k } ) ^ { g } } \right ) \ ( 1 2 a ) \\ \text{with integration contour radii} \, \rho _ { a } > 0 \text{ restricted as follows,} \\ & \quad \cdot \quad \, - \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \end{$$

$$\rho _ { a + 1 } < \rho _ { a } \ \ ( a = 1, \dots, L - 2 ), \ \ 0 < | z _ { j } | < \rho _ { L - 1 } \ \ ( j = 1, \dots, n ), \quad \ \ ( 1 2 \text{b} )$$

and ρ L -1 &gt; 0 arbitrary. Note that, if and only if r = ( r , . . . , r 1 L ) satisfies the condition r 1 ≥ · · · ≥ r L , then λ ( a ) -( r N a a +1 ) is a partitions of length ≤ N a for all a = 1 , . . . , L -1, and therefore is well-defined ( N λ,n ( g ) and b λ ( g ) are given in (7) and (9), respectively). Moreover, then the Jack polynomial P (1 /g ) λ,n ( z ) with λ = λ ( L ) exists and is non-zero.

$$C _ { L } ( r ; s ) & \equiv \prod _ { a = 1 } ^ { L - 1 } \frac { N _ { a }! \mathcal { N } _ { \lambda ^ { ( a ) }, N _ { a } } ( g ) } { b _ { \lambda ^ { ( a ) } - ( r _ { a + 1 } ^ { N _ { a } } ) } ( g ) } } \\ _ { x } ( a ) \text{ and } b _ { x } ( a ) \text{ are given in } ( 7 ) \text{ and } ( 9 ). \text{respectively} ). \text{ Moreover, then}$$

Proposition 2.2 (Awata et al. [5]) . Let L ∈ Z ≥ 1 , r ∈ Z L , s ∈ Z L ≥ 1 , λ ( a ) as in (11) , and ρ L -1 &gt; 0 arbitrary. Then the following hold true.

- (a) The integrals P r,s,L ( z , . . . , z 1 n ) in (12a) -(12b) are well-defined analytic functions of ( z , . . . , z 1 n ) in the region 0 &lt; z | j | &lt; ρ L -1 ( j = 1 , . . . , n ) .
- (b) If r = ( r , . . . , r 1 L ) ∈ Z L satisfies the condition r 1 ≥ · · · ≥ r L , then

$$\mathcal { P } _ { r, s, L } ( z _ { 1 }, \dots, z _ { n } ) = C _ { L } ( r ; s ) P _ { \lambda, n } ^ { ( 1 / g ) } ( z _ { 1 }, \dots, z _ { n } ) \ \ ( \lambda = \lambda ^ { ( L ) } ).$$

Otherwise, P r,s,L ( z , . . . , z 1 n ) = 0 .

Proof. See Appendix B.

Every partition λ = ( λ , . . . , λ 1 n ) ∈ P n can be represented as

$$\lambda = ( r _ { 1 } ^ { s _ { 1 } }, \dots, r _ { L } ^ { s _ { L } } )$$

with L = n , r = λ , and s = (1 ), and thus n P (1 /g ) λ,n ( z ) = P λ, (1 n ) ,n ( z , . . . , z 1 n ) /C n ( λ ; (1 n )) provides an integral representation for every Jack polynomials. If λ is such that all λ j are distinct this is the only integral representation of the corresponding Jack polynomial

provided by Proposition 2.2. However, if some of the λ j are identical, there are several different ways to write λ as in (14), and thus one gets from Proposition 2.2 several different integral representations P (1 /g ) λ,n ( z ) = P r,s,L ( z , . . . , z 1 n ) /C L ( r ; s ) of one and the same Jack polynomial. The number of integrations in such a representation is

$$\sum _ { a = 1 } ^ { L - 1 } N _ { a } & = \sum _ { j = 1 } ^ { L - 1 } ( L - j ) s _ { j }, \\. \cdots \cdots, \cdots \sigma. \cdots,$$

and thus different such representations have different complexity.

Note that the integral in (3) is a natural elliptic generalization of the one in (12a) in the special case L = n , r = λ , and s = (1 ). n

## 3 Result

We give a mathematically precise formulation of the main result in this paper.

We introduce the Euler operator

$$D _ { n } ( x ) & \equiv - \mathfrak { i } \sum _ { j = 1 } ^ { n } \frac { \partial } { \partial x _ { j } }. \\. &. &. &. &. &..$$

To simplify some formulas later on we introduce the notation and

$$V ( x ; \tau ) \equiv \wp ( x | \pi, \pi \tau ) + \frac { \eta _ { 1 } ( \tau ) } { \pi } = \sum _ { m \in \mathbb { Z } } \frac { 1 } { 4 \sin ^ { 2 } \left ( \frac { 1 } { 2 } x + \pi m \tau \right ) } \quad ( x \in \mathbb { C } ) \\ \text{and} \\ / \cdots \quad \cdots \quad \infty \quad \cdot$$

$$\frac { \eta _ { 1 } ( \tau ) } { \pi } \equiv \frac { 1 } { 1 2 } + \sum _ { m = 1 } ^ { \infty } \frac { 1 } { 2 \sin ^ { 2 } ( m \pi \tau ) } \\ \text{Moreover, we find it convenient to write the non-stationary eCS equation}$$

$$& \text{as} \quad \left ( \kappa \frac { \text{i} } { \pi } \frac { \partial } { \partial \tau } - \sum _ { j = 1 } ^ { n } \frac { \partial ^ { 2 } } { \partial x _ { j } ^ { 2 } } + \sum _ { j < k } ^ { n } 2 g ( g - 1 ) V ( x _ { j } - x _ { k } ; \tau ) \right ) \psi ( x ; \tau ) = E ( \tau ) \psi ( x ; \tau ) \quad \quad ( 1 8 ) \\ & \text{where $E(\tau)$ is an eigenvalue that we allow to be non-zero.}$$

(see Appendix A). Moreover, we find it convenient to write the non-stationary eCS equation as where E τ ( ) is an eigenvalue that we allow to be non-zero.

To explain the significance of this eigenvalue, we note that (18) is invariant under transformations

$$\psi ( x ; \tau ) \rightarrow C ( \tau ) \psi ( x ; \tau ), \ \ E ( \tau ) \rightarrow E ( \tau ) + \kappa \frac { \mathrm i } { \pi } \frac { 1 } { C ( \tau ) } \frac { \partial C ( \tau ) } { \partial \tau }$$

for arbitrary differentiable non-zero functions C τ ( ). Thus, for non-zero κ , one can transform E τ ( ) to zero, or any other value, by such a transformation. In particular, the differential equation in (18) can be transformed to (1) by a change of the normalization of ψ x τ ( ; ).

However, we find it convenient to fix the normalization of ψ x τ ( ; ) so that the limits κ → 0 and /Ifractur ( τ ) → + ∞ are well-defined and non-trivial. We thus have to work with (18) rather than with (1). As will be seen, our normalization is such that E τ ( ) is equal to the eigenvalues of the trigonometric Calogero-Sutherland model [2] in the limit /Ifractur ( τ ) → + ∞ .

As mentioned in Section 1, we obtain explicit solutions of the non-stationary eCS equation in (18) for κ = kg and k ∈ Z ≥ 1 that provide elliptic generalizations of the Jack polynomials P (1 /g ) λ,n ( z ). When k = 1 then this solution can be constructed for all variable numbers n ∈ Z ≥ 1 and λ ∈ Z n and is therefore complete in this sense, whereas, for k ≥ 2, n and λ are restricted as follows, or

$$n = k L, \quad \lambda = ( \underbrace { r _ { 1 }, \dots, r _ { 1 } } _ { k \text{ times} }, \dots, \underbrace { r _ { L }, \dots, r _ { L } } _ { k \text{ times} } )$$

$$n = 1 + k ( L - 1 ), \quad \lambda = ( r _ { 1 }, \underbrace { r _ { 2 }, \dots, r _ { 2 } } _ { k \text{ times} }, \dots, \underbrace { r _ { L }, \dots, r _ { L } } _ { k \text{ times} } ) \\ n = k L, \quad \lambda = ( r _ { 1 }, \dots, r _ { 1 }, \dots, r _ { L }, \dots, r _ { L } )$$

︸ ︷︷ ︸ ︸ ︷︷ ︸ for arbitrary L ∈ Z ≥ 1 and integer vectors r = ( r , . . . , r 1 L ) ∈ Z L . Using the notation in (4) we can distinguish these two cases by a parameter s 1 equal to 1 in the first case and k in the second case. Then both cases can be summarized as follows,

$$n = s _ { 1 } + k ( L - 1 ), \quad \lambda = ( \lambda _ { 1 }, \dots, \lambda _ { n } ) = ( r _ { 1 } ^ { s _ { 1 } }, r _ { 2 } ^ { k }, \dots, r _ { L } ^ { k } ).$$

This suggests to label our solutions by the triple ( r, s, L ) with s = ( s , k 1 L -1 ).

The pertinent integrals are given by

/negationslash

/negationslash with N a = s 1 +( a -1) k ( a = 1 , . . . , L -1) and integration radii ρ a constrained as follows,

$$\Pi _ { r, s, L } ( z _ { 1 }, \dots, z _ { n } ; p ) & \equiv \left ( \prod _ { j = 1 } ^ { n } z _ { j } ^ { r _ { L } } \right ) \left ( \prod _ { a = 1 } ^ { L - 1 } \prod _ { j = 1 } ^ { N _ { a } } \oint _ { | \epsilon _ { a j } | = \rho _ { a } } \frac { d \xi _ { a j } } { 2 \pi i \xi _ { a j } } \xi _ { a j } ^ { r _ { a - r _ { a + 1 } } } \right ) \\ & \times \left ( \frac { \prod _ { j \neq k } ^ { N _ { L - 1 } } \theta ( \xi _ { N _ { L - 1 } 1 } j ) } { \prod _ { j = 1 } ^ { N _ { L - 1 } } \theta ( z _ { j } ) } \left ( \prod _ { a = 1 } ^ { L - 2 } \prod _ { j \neq k } \theta ( \xi _ { a j } / \xi _ { a k } ; p ) ^ { g } } { \prod _ { a = 1 } ^ { N _ { a + 1 } } \prod _ { j = 1 } ^ { N _ { a } } \theta ( \xi _ { a + 1, j } / \xi _ { a k } ; p ) ^ { g } } \right ) \ ( 2 1 a ) \\ \text{with } N _ { a } = s _ { 1 } + ( a - 1 ) k \ ( a = 1, \dots, L - 1 ) \text{ and integration radii } \rho _ { a } \text{ constrained as follows,}$$

$$| p | \rho _ { a } < \rho _ { a + 1 } < \rho _ { a } \ \ ( a = 1, \dots, L - 2 ), \ \ | p | \rho _ { L - 1 } < | z _ { j } | < \rho _ { L - 1 } \ \ \ ( 2 1 b )$$

for some arbitrary ρ L -1 in the range 1 &lt; ρ L -1 &lt; 1 / p | | . Recall that the Jack polynomials are well-defined for all ordered integer vectors λ (see Remark 2.1).

Theorem 3.1. Let g &gt; 1 / 2 , k ∈ Z ≥ 1 , L ∈ Z ≥ 1 , r = ( r , . . . , r 1 L ) ∈ Z L , s 1 = 1 for k = 1 and s 1 ∈ { 1 , k } for k ≥ 2 , s = ( s , k 1 L -1 ) , and ρ L -1 arbitrary in the range 1 &lt; ρ L -1 &lt; 1 / p | | . Then the following hold true.

(a) The integrals P r,s,L ( z , . . . , z 1 n ; p ) in (21a) -(21b) are analytic functions of ( z, p ) = ( z , . . . , z 1 n , p ) in the following region ⊂ C n +1 ,

$$| p | \rho _ { L - 1 } < | z _ { j } | < \rho _ { L - 1 } \ \ ( j = 1, \dots, n ), \ \ | p | < 1.$$

## (b) The functions

/negationslash

$$\psi _ { r, s, L } ( x _ { 1 }, \dots, x _ { n } ; \tau ) \equiv \left ( \prod _ { j \neq k } ^ { n } \theta \left ( z _ { j } / z _ { k } ; p \right ) \right ) ^ { g / 2 } \mathcal { P } _ { r, s, L } ( z _ { 1 }, \dots, z _ { n } ; p ) \quad ( z _ { j } = e ^ { i x _ { j } }, \, p = e ^ { 2 \pi i r } ) \\ \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \ dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots$$

are solutions of the non-stationary eCS equation in (18) for κ = kg with the following eigenvalue determined by the integer vector λ in (20) ,

$$E _ { r, s, L } ( \tau ) = \sum _ { j = 1 } ^ { n } \left ( \lambda _ { j } + \frac { 1 } { 2 } g ( n + 1 - 2 j ) \right ) ^ { 2 } + g ^ { 2 } n ( n - 1 ) \left ( \frac { \eta _ { 1 } ( \tau ) } { \pi } - \frac { 1 } { 1 2 } \right ). \quad \quad ( 2 2 ) \\ M o r e o n e r$$

Moreover,

$$\ e r, D _ { n } ( x ) \psi _ { r, s, L } ( x ; \tau ) = d _ { r, s, L } \psi _ { \lambda } ( x ; \tau ), \quad d _ { r, s, L } = \sum _ { j = 1 } ^ { n } \lambda _ { j }.$$

(c) If r 1 ≥ ··· ≥ r L , then the integer vector λ in (20) is ordered, and the integral P r,s,L ( z , . . . , z 1 n ; p ) in (21a) -(21b) defines an elliptic generalization of the Jack polynomial P (1 /g ) λ,n ( z ) in the following sense,

$$\lim _ { p \to 0 } \mathcal { P } _ { r, s, L } ( z _ { 1 }, \dots, z _ { n } ; p ) = C _ { L } ( r ; s ) P _ { \lambda, n } ^ { ( 1 / g ) } ( z _ { 1 }, \dots, z _ { n } )$$

with the constant C L ( r ; s ) in (13) . Otherwise, lim p → 0 P r,s,L ( z , . . . , z 1 n ; p ) = 0 .

Thus, by varying ρ L -1 in the allowed range, the integral in (21a)-(21b) defines a function in the region | p | &lt; z j &lt; 1 / p | | ( j = 1 , . . . , n ), | p | &lt; 1. It is interesting to note that the following integration radii satisfy the conditions in (21b),

$$\rho _ { a } = \rho _ { L - 1 } \epsilon ^ { 1 - L + a } \ \ ( a = 1, \dots, L - 1 ),$$

for arbitrary /epsilon1 and ρ L -1 in the range | p | &lt; /epsilon1 &lt; 1 and 1 &lt; ρ L -1 &lt; 1 / p | | .

Part (c) of Theorem 3.1 is a corollary of the integral representation of the Jack polynomials [5] reviewed in Section 2.2. Indeed, lim p → 0 θ z ( ; p ) = (1 -z ) implies that the integrals in (21a)-(21b) reduce to the ones in (12a)-(12b) in the limit p → 0, and thus Part (c) is a simple consequence of Part (b) of Proposition 2.2.

We therefore only need to prove Parts (a) and (b) of this theorem. This is done in Section 4.

## 4 Proof

In this section we prove Parts (a) and (b) of Theorem 3.1 (Part (c) is a corollary of Proposition 2.2). After some preliminaries in Section 4.1, we give an outline in Section 4.2 where we explain our proof strategy, suppressing technical details. The detailed proof is given in Section 4.3 with some technical parts deferred to Appendix C.

## 4.1 Definitions and preliminary remarks

Recall the definitions of θ z ( ; p ), V ( x τ ; ), and η 1 ( τ ) /π in (2), (16), and (17), respectively.

We recall the definition of the Euler operator D n ( x ) in (15). We also introduce the short-hand notations

/negationslash

$$\Psi _ { n } ( x ; \tau ) \equiv \left ( \prod _ { j \neq k } ^ { n } \theta \left ( z _ { j } / z _ { k } ; p \right ) \right ) ^ { g / 2 } \quad ( z _ { j } = e ^ { i x _ { j } }, \ p = \mathrm e ^ { 2 \pi i \tau } )$$

and

With that we can write the non-stationary eCS equation in (18) short as follows,

$$H _ { n } ( x ; \tau ) \equiv - \sum _ { j = 1 } ^ { n } \frac { \partial ^ { 2 } } { \partial x _ { j } ^ { 2 } } + \sum _ { j < k } ^ { n } 2 g ( g - 1 ) V ( x _ { j } - x _ { k } ; \tau ). \intertext {hat, we can write the non-stationary eCS equation in (1.8) short, as follows}$$

$$\left ( \kappa \frac { \mathrm i } { \pi } \frac { \partial } { \partial \tau } + H _ { n } ( x ; \tau ) \, \right ) \psi ( x ; \tau ) = E ( \tau ) \psi ( x ; \tau ). \quad \right ) \quad \left ( 2 6 \right ) \\ \ast \, \ \nu \quad \dots \quad \neq \dots \quad \dots \quad \dots \quad \dots \quad \dots \,.$$

Moreover, Ψ ( n x τ ; ) is a common factor in all our solutions, i.e., the function of interest P ( z ; p ) is such that ψ x τ ( ; ) = Ψ ( n x τ ; ) P ( z ; p ).

To motivate our notation we mention in passing that the differential operator H n ( x τ ; ) in (25) has a natural physical interpretation as Hamiltonian of a quantum many-body system describing n particles moving on a circle and interacting with two-body interactions given by the Weierstrass ℘ -function, ψ x τ ( ; ) is a quantum mechanical wave function, and D n ( x ) in (15) can be interpreted as momentum operator.

## 4.2 Outline of proof of Theorem 3.1

We explain the strategy we use to construct the solutions in Theorem 3.1.

From now on we simplify notation by suppressing the τ -dependence in our equations if there is no danger of confusion, e.g., we write H n ( x ) short for H n ( x τ ; ) etc.

A simple but important observation is the following: since the operators κ i π ∂τ ∂ + H n ( x ) and D n ( x ) obviously commute, we can construct our solutions ψ x ( ) of (26) such that they are also eigenfunctions of D n ( x ).

The key to our result is a generalized kernel function K NM ( x, y ) for the elliptic CalogeroSutherland model obtained in [11]. This function is explicitly known (see (38)), and it satisfies the functional identities [11]

$$\left ( D _ { N } ( x ) + D _ { M } ( y ) \right ) K _ { N M } ( x, y ) = 0,$$

$$\left ( ( N - M ) g \frac { \mathrm i } { \pi } \frac { \partial } { \partial \tau } + H _ { N } ( x ) - H _ { M } ( y ) - c _ { N M } \right ) K _ { N M } ( x, y ) = 0 \quad \text{(28)} \quad \text{if the known constant $c_{x}$ see Lemma $4$ 1 for the nrecise statement}$$

with a known constant c NM ; see Lemma 4.1 for the precise statement.

We fix N -M = k ∈ Z ≥ 1 and use this generalized kernel function to construct an integral transform that maps a given solution ψ M ( y ) of

$$D _ { M } ( y ) \psi _ { M } ( y ) & = d _ { M } \psi _ { M } ( y ), \quad \left ( k g \frac { \mathrm i } { \pi } \frac { \partial } { \partial \tau } + H _ { M } ( y ) \ \right ) \psi _ { M } ( y ) = E _ { M } \psi _ { M } ( y ) \\ \dots \, \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \dots \quad \cdot \end{cases}$$

with known eigenvalues E M and d M , to a solution of

$$D _ { N } ( x ) \psi _ { N } ( x ) = d _ { N } \psi _ { N } ( x ), \quad \left ( k g \frac { \mathfrak { i } } { \pi } \frac { \partial } { \partial \tau } + H _ { N } ( x ) \right ) \psi _ { N } ( x ) = E _ { N } \psi _ { N } ( x ) \quad \quad ( 3 0 ) \\ \dots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots$$

with N = M + k and eigenvalues d N and E N that can be computed. Since this new solution ψ N ( x ) has a larger number of variables, we can use such integral transforms to inductively construct solutions with an increasing sequence of variable numbers N a = N 1 +( a -1) k , a = 2 , . . . , L , starting out with a variable number N 1 which is such that explicit solutions can be constructed by other methods.

To obtain a variety of solutions it is important have free parameters in these integral transformations. We thus note that, if ψ N ( x ), d N , and E N are solutions of the equations in (30), then

$$\psi _ { N } ^ { Q } ( x ) = e ^ { i Q \sum _ { j = 1 } ^ { N } x _ { j } } \psi _ { N } ( x ), \ \ d _ { N } ^ { Q } = d _ { N } + N Q, \ \ E _ { N } ^ { Q } = E _ { N } + 2 Q d _ { N } + N Q ^ { 2 } \quad \ ( 3 1 )$$

are solutions of the same equations for arbitrary real Q (this can be checked by straightforward computations). This motivates the following ansatz for the integral transform:

$$( \mathcal { K } _ { N M } ^ { Q Q ^ { \prime } } \psi _ { M } ) & ( x ) = C _ { N M } e ^ { i Q \sum _ { j = 1 } ^ { N } x _ { j } } \int \frac { d ^ { M } y } { ( 2 \pi ) ^ { M } } K _ { N M } ( x, y ) e ^ { - i Q ^ { \prime } \sum _ { k = 1 } ^ { M } y _ { j } } \psi _ { M } ( y ) \\ \dots \quad. \quad \sim \sim \dots \sim \dots \dots \dots \dots \dots \dots \dots$$

with real parameters Q,Q ′ and C NM a non-zero constant to be determined. We stress that (32) is formal since we neither specified the integration domain nor the parameters Q,Q ′ , and these details are important to get well-defined solutions. However, if we ignore such details for now, we can proceed with formal computations using (27)-(29) and dropping boundary terms obtained by partial integrations used to transfer the action of the differential operator (details on this computation are given in the proof of Lemma 4.2 below). We thus find that ψ N ( x ) = ( K QQ ′ NM ψ M )( x ) satisfies (30) with

$$d _ { N } = d _ { M } + N Q - M Q ^ { \prime }, \ \ E _ { N } = E _ { M } + 2 ( Q - Q ^ { \prime } ) d _ { M } + M ( Q ^ { \prime } ) ^ { 2 } - 2 M Q Q ^ { \prime } + N Q ^ { 2 } + c _ { N M } ; \ ( 3 3 )$$

see Lemma 4.2 below for the mathematically precise statement.

Turning now to our inductive scheme, we use two possible base cases N = N 1 , namely N 1 = 1 and N 1 = k (for k = 1 these two cases coincide). As explained later on, solutions of (30) in these cases are given by

$$N _ { 1 } = 1 \colon \ \psi _ { N _ { 1 } } ( x _ { 1 } ) = z _ { 1 } ^ { r _ { 1 } }, \ \ d _ { N _ { 1 } } = r _ { 1 }, \ \ E _ { N _ { 1 } } = r _ { 1 } ^ { 2 }$$

and

$$N _ { 1 } = k \, \colon \quad \psi _ { N _ { 1 } } ( x ) = \Psi _ { k } ( x ) ( z _ { 1 } \cdots z _ { k } ) ^ { r _ { 1 } }, \quad d _ { N _ { 1 } } = k r _ { 1 }, \quad E _ { N$$

for arbitrary integer r 1 , with Ψ ( k x ) the function in (24) for n = k and c k 0 the constant in (28) for ( N,M ) = ( k, 0) (recall that z j = e i x j ). The induction step from N a to N a +1 = N a + k is given by the integral transform described above,

$$\psi _ { N _ { a } } ( x ) = \left ( \mathcal { K } _ { N _ { a } N _ { a - 1 } } ^ { Q _ { a } Q ^ { \prime } _ { a } } \psi _ { N _ { a - 1 } } \right ) ( x ) \\ \text{We find that, at each such step$$

for a = 2 3 , , . . . , L . We find that, at each such step, the possible parameters Q a and Q ′ a are determined by integer parameters r a and r ′ a , and one can set r ′ a = r a , without loss of generality. Our solutions ψ r,s,L ( x ) are obtained as ψ N L ( x ).

With the outline of our methods explained, we now turn to making these arguments precise.

## 4.3 Proof of Theorem 3.1

We prove Theorem 3.1 following the outline in Section 4.2.

## 4.3.1 Generalized kernel function identity

The result in this section was obtained in Ref. [11].

Recall the definitions of D n ( x ) and H n ( x ) in (15) and (25), respectively. We define

$$\vartheta ( x ) \equiv 2 \sin \left ( \frac { 1 } { 2 } x \right ) \prod _ { m = 1 } ^ { \infty } \left ( 1 - 2 p ^ { m } \cos ( x ) + p ^ { 2 m } \right ) \quad ( x \in \mathbb { C }, \ p = e ^$$

proportional to the Jacobi theta function ϑ 1 ( 1 2 x τ | ) (see Appendix A). This allows us to write the function in (24) as

$$\Psi _ { N } ( x ) & = \left ( \prod _ { j < k } ^ { N } \vartheta \left ( x _ { j } - x _ { k } \right ) ^ { 2 } \right ) ^ { g / 2 } \quad ( x \in \mathbb { C } ^ { N } ) \quad \quad \quad$$

(this is true since ϑ x ( j -x k ) 2 = θ z /z ( j k ; p θ z /z ) ( k j ; p ) for | z /z j k | = 1; note that we always assume that this latter condition is fulfilled). With that we can define the generalized kernel function mentioned in the previous section as

$$K _ { N M } ( x, y ) \equiv \frac { \Psi _ { N } ( x ) \Psi _ { M } ( y ) } { \left ( \prod _ { j = 1 } ^ { N } \prod _ { k = 1 } ^ { M } \vartheta \left ( x _ { j } - y _ { k } \right ) \right ) ^ { g } } \\ \dots, x _ { N } \right ) \in \mathbb { C } ^ { N } \text{ and } y \in ( y _ { 1 }, \dots, y _ { M } ) \in \mathbb { C } ^ { M }.$$

with x = ( x , . . . , x 1 N ) ∈ C N and y ∈ ( y , . . . , y 1 M ) ∈ C M .

Lemma 4.1. For all N,M ∈ Z ≥ 0 and g ∈ C , the function K NM ( x, y ) in (38) satisfies the functional identities in (27) -(28) with the constant

$$c _ { N M } = g ^ { 2 } [ N ( N - 1 ) - M ( M - 1 ) ] \frac { \eta _ { 1 } } { \pi } + \frac { 1 } { 1 2 } g ^ { 2 } ( N - M ) [ N ( N - 1 ) + M ( M - 1 ) - 2 N M ].$$

Proof. This result is stated and proved in Ref. [11]; see first part of the Theorem. For the convenience of the reader we give a table translating the symbols used here to the ones used in Ref. [11]:

| Symbols used here    | ϑ ( x )   | g   | τ       | p   | η 1 π   | H N ( x )   | D N ( x )   | K NM ( x, y )   |
|----------------------|-----------|-----|---------|-----|---------|-------------|-------------|-----------------|
| Symbols used in [11] | 2 θ ( x ) | λ   | i β 2 π | q 2 | c 0     | H λ,N ( x ) | - P N ( x ) | F N,M ( x ; y ) |

(we ignore an irrelevant multiplicative factor between K NM ( x, y ) and F N,M ( x y ; ) coming from the factor 2 between ϑ x ( ) and θ ( x )).

A careful reader will notice that the proof in [11] is for the function Ψ N ( x ) given by

$$\left ( \prod _ { j < k } ^ { N } \vartheta \left ( x _ { j } - x _ { k } \right ) \right ) ^ { g } \\ \cdots \vdots \vdots \int _ { 0 } \in \vdots \int _ { 1 } \in \vdots \int _ { 2 } \in \vdots \vdots \vdots$$

differing from the function Ψ N ( x ) in (37) by a phase factor. However, this phase factor is locally constant (it only changes if one crosses branch cuts), and it therefore does not affect the proof of Lemma 4.1 in [11].

## 4.3.2 Integral operator

To make the formal integral in (32) precise, we replace

$$\dots \quad \int \frac { d ^ { M } y } { ( 2 \pi ) ^ { M } } \to \int _ { [ - \pi - i \epsilon, \pi - i \epsilon ] ^ { M } } \frac { d ^ { M } y } { ( 2 \pi ) ^ { M } } \equiv \int _ { - \pi - i \epsilon } ^ { \pi - i \epsilon } \frac { d y _ { 1 } } { 2 \pi } \cdots \int _ { - \pi - i \epsilon } ^ { \pi - i \epsilon } \frac { d y _ { M } } { 2 \pi } \quad \quad \quad ( 4 0 ) \\ \dots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \end{cases}$$

with /epsilon1 a parameter further constrained below. This parameter /epsilon1 allows us to deform the integration paths in the complex plane so as to avoid singularities of the integrand, and by this we can make the computations described after (32) mathematically precise. By straightforward computations we obtain explicit formulas for the transformed function ψ N ( x ) = ( K QQ ′ NM ψ M )( x ), which is well-defined provided ψ M ( y ) satisfies certain conditions. We also find that it is convenient to choose the constant in (32) as

$$C _ { N M } = e ^ { \pi i g N M / 2 }.$$

We summarize our result as follows.

Lemma 4.2. Let g &gt; 1 / 2 , M,N ∈ Z ≥ 1 with N &gt; M ,

$$\psi _ { M } ( y ; \tau ) = \Psi _ { M } ( y ; \tau ) \mathcal { P } _ { M } ( \xi ; p ) \quad ( \xi = ( e ^ { i y _ { 1 } }, \dots, e ^ { i y _ { M } } ), \quad p = e ^ { 2 \pi i \tau } ) \quad \quad \quad$$

with P M ( ξ ; p ) an analytic function of ( ξ, p ) = ( ξ , . . . , ξ 1 M , p ) in the following region ⊂ C M +1 ,

$$\rho _ { 0 } | p | < | \xi _ { j } | < \rho _ { 0 } \ \ ( 1 \leq j \leq M ), \ \ | p | < 1$$

for some ρ 0 &gt; 0 , and

$$Q = r - M g / 2, \quad Q ^ { \prime } = r ^ { \prime } - N g / 2$$

with r, r ′ ∈ Z . Then the integral transform ( K QQ ′ NM ψ M )( x τ ; ) of ψ M ( y ; τ ) in (32) becomes well-defined by the substitution in (40) provided that ρ ≡ e /epsilon1 is in the range | p ρ | 0 &lt; ρ &lt; ρ 0 . Furthermore, with z = (e i x 1 , . . . , e i x N ) where

$$\bar { \nu }, \quad & ( \mathcal { K } _ { N M } ^ { Q Q ^ { \prime } } \psi _ { M } ) ( x ; \tau ) = \Psi _ { N } ( x ; \tau ) ( \hat { \mathcal { K } } _ { N M } ^ { r r ^ { \prime } } \mathcal { P } _ { M } ) ( z ; p$$

/negationslash is an analytic function of ( z, p ) = ( z , . . . , z 1 N , p ) in the following region ⊂ C N +1 ,

$$( \hat { K } _ { N M } ^ { r r ^ { \prime } } \mathcal { P } _ { M } ) ( z ; p ) = \left ( \prod _ { j = 1 } ^ { N } z _ { j } ^ { r } \right ) \left ( \prod _ { k = 1 } ^ { M } \oint _ { | \xi _ { k } | = \rho } \frac { d$$

$$| p | \rho < | z _ { j } | < \rho \quad ( 1$$

Finally, if ψ M ( y ; τ ) satisfies the equations in (29) with known eigenvalues d M and E M , then ψ N ( x τ ; ) = ( K QQ ′ NM ψ M )( x τ ; ) satisfies the equations in (30) with eigenvalues d N and E N given in (33) , (39) and (44) .

Proof. Evaluating (32), using (38), (40) and (42), we obtain

$$( \mathcal { K } _ { N M } ^ { Q Q ^ { \prime } }$$

We change variables to z j = e i x j and ξ k = e i y k and recall the definition of Ψ N ( x ) in (24). Inserting ϑ x ( j -y k ) = e π / i 2 e -i( x j -y k ) / 2 θ z /ξ ( j k ), which follow from the definitions in (2) and (36) for | p | &lt; z /ξ | j k | &lt; 1, we obtain, by straightforward computations,

/negationslash

$$( \mathcal { K } _ { N M } ^ { Q Q ^ { \prime } }$$

with ρ = e . /epsilon1 We insert C NM in (41) and Q,Q ′ in (44) to obtain (45)-(46).

$$\ll \kappa _ { N M } \gamma _ { N } / \ll \k$$

The proof that the integral transform in (46) defines a analytic function is deferred to Appendix C.1.

We are left to prove that, if ψ M ( y ) satisfies (29), then ψ N ( x ) in (32) satisfies (30). To simplify notation, we write in the following ∫ short for ∫ [ - -π i /epsilon1,π -i /epsilon1 ] M , and we use the short hand notations | x | = ∑ j x j and similarly for | y | . We start with

$$\begin{array} {$$

using (27) in the last equality, where the interchange of integration and differentiations in the second equality is justified by the Leibniz integral rule 10 using the properties of the integrand established above (we use {· · · } to indicate that derivatives act only within this bracket, here and in the following). We compute this by partial integrations using that the boundary term vanishes, i.e.,

$$\begin{array} { c } \check { \check { \cdot } } \, \check { \cdot } \, \check { \cdot } \, \check { \cdot } \, \check { \cdot } \, \check { \cdot } \, \check { \cdot } \, \check { \cdot } \, \check { \cdot } \, \check { \cdot } \, \check { \cdot } \, \check { \cdot } \, \check { \cdot } \, \check { \cdot } \, \check { \cdot } \, \check { \cdot } \, \check { \cdot } \, \check { \cdot } \, \check { \cdot } \, \check { \cdot } \, \check { \cdot } \, \check { \cdot } \, \check { \cdot } \, \check { \cdot } \\ \int \frac { d ^ { M } y } { ( 2 \pi ) ^ { M } } \left ( \{ D _ { M } ( y ) K _ { N M } ( x, y ) \} \text{e} ^ { - i Q ^ { \prime } | y | } \psi _ { M } ( y ) + K _ { N M } ( x, y ) D _ { M } ( y ) \text{e} ^ { - i Q ^ { \prime } | y | } \psi _ { M } ( y ) \right ) = 0 \ \ ( 4 8 ) \\ \text{(this is owned in Annendix C.2.)}. \ \text{We thus obtain} \end{array}$$

(this is proved in Appendix C.2). We thus obtain

$$D _ { i }$$

$$= ( d _ { M } + N Q - M Q ^ { \prime } ) \psi _ { N } ( x ) \\ \dots \dots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots$$

$$\left ( \begin{array} { c } \left ( \begin{array} { c } \left ( \end{array} { c } \right ) \end{array} \right ) \, \psi \, \nu \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \sum \, \dots \\ D _ { N } ( x ) \psi _ { N } ( x ) = C _ { N M } \right ) \right ) ^ { i Q | x | } \int \frac { d ^ { M } y } { ( 2 \pi ) ^ { M } } K _ { N M } ( x, y ) \left ( D _ { M } ( y ) + N Q \right ) \mathrm e ^ { - i Q ^ { \prime } | y | } \left ( D _ { M } ( y ) + N Q - M Q ^ { \prime } \right ) \psi _ { M } ( y ) \\ = ( d _ { M } + N Q - M Q ^ { \prime } ) \psi _ { N } ( x ) \end{array}$$

using D N ( y ψ ) M ( y ) = d M ψ M ( y ) in the last step. This proves the first identity in (30) with d N in (33).

In a similar manner,

$$& \text{In a similar manner,} \\ & \left ( k g \frac { i } { \pi } \frac { \partial } { \partial \tau } + H _ { N } ( x ) \right ) \psi _ { N } ( x ) = C _ { M M } e ^ { i Q | x | } \int \frac { d ^ { M } y } { ( 2 \pi ) ^ { M } } \left ( k g \frac { i } { \pi } \frac { \partial } { \partial \tau } + H _ { N } ( x ) \\ & + 2 Q D _ { N } ( x ) + N Q ^ { 2 } \right ) K _ { N M } ( x, y ) e ^ { - i Q ^ { \prime } | y | } \psi _ { M } ( y ) = C _ { N M } e ^ { i Q | x | } \int \frac { d ^ { M } y } { ( 2 \pi ) ^ { M } } \left [ \left \{ \left ( H _ { M } ( y ) + c _ { N M } \\ & \quad - 2 Q D _ { M } ( y ) + N Q ^ { 2 } \right ) K _ { N M } ( x, y ) \right \} e ^ { - i Q ^ { \prime } | y | } \psi _ { M } ( y ) + K _ { N M } ( x, y ) k g \frac { i } { \pi } \frac { \partial } { \partial \tau } e ^ { - i Q ^ { \prime } | y | } \psi _ { M } ( y ) \right ] \\ & \text{using the Leibniz integral rule and (27)-(28) \text{ to get the second identity (recall that } k = } \text{$\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\ tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau $} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$} \text{ $\tau$}$$

using the Leibniz integral rule and (27)-(28) to get the second identity (recall that k = N -M ). Again, we proceed by partial integrations using that that the boundary terms which arise vanish, i.e., (48) and

$$& \int \frac { d ^ { M } y } { ( 2 \pi ) ^ { M } } \left ( \{ H _ { M } ( y ) K _ { N M } ( x, y ) \} \, \text{e} ^ { - i Q ^ { \prime } | y | } \psi _ { M } ( y ) - K _ { N M } ( x, y ) H _ { M } ( y ) e ^ { - i Q ^ { \prime } | y | } \psi _ { M } ( y ) \right ) = 0 \ \ ( 4 9 ) \\ & \text{(this is nroved in Amendiv C. 2.)} \, \text{This.}$$

(this is proved in Appendix C.2). Thus,

$$& ( \text{this is proved in Appendix C.2}). \text{ Thus,} \\ & \quad \left ( k g \frac { i } { \pi } \frac { \partial } { \partial \tau } + H _ { N } ( x ) \right ) \psi _ { N } ( x ) = C _ { M } e ^ { i Q | x | } \int \frac { d ^ { M } y } { ( 2 \pi ) ^ { M } } K _ { N M } ( x, y ) \left ( k g \frac { i } { \pi } \frac { \partial } { \partial \tau } + H _ { M } ( y ) + c _ { N M } \\ & \quad + 2 Q D _ { M } ( y ) + N Q ^ { 2 } \right ) e ^ { - i Q ^ { \prime } | y | } \psi _ { M } ( y ) = C _ { N M } e ^ { i Q | x | } \int \frac { d ^ { M } y } { ( 2 \pi ) ^ { M } } K _ { N M } e ^ { - i Q ^ { \prime } | y | } \left ( k g \frac { i } { \pi } \frac { \partial } { \partial \tau } + H _ { M } ( y ) \\ & \quad - 2 Q ^ { \prime } D _ { M } ( y ) + M ( Q ^ { \prime } ) ^ { 2 } + c _ { N M } + 2 Q D _ { M } ( y ) - 2 M Q Q ^ { \prime } + N Q ^ { 2 } \right ) \psi _ { M } ( y ) \\ & \quad = ( E _ { M } + 2 ( Q - Q ^ { \prime } ) d _ { M } + M ( Q ^ { \prime } ) ^ { 2 } - 2 M Q Q ^ { \prime } + N Q ^ { 2 } + c _ { N M } ) \psi _ { N } ( x ) \\ & \quad \text{using (29) in the last step. This proves the second identity in (30) with E _ { N } \text{ in (33)}. \quad \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \tau \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \ \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \quad - \text{$\Delta$} \ }$$

using (29) in the last step. This proves the second identity in (30) with E N in (33).

10 See e.g. [22, Section 4.2] for a precise mathematical formulation of the Leibniz integral rule. Note that, in all cases where we use this rule in this paper, the pertinent function in the integrand satisfies analyticity properties that are much stronger than what is needed to prove the validity of this rule.

## 4.3.3 Induction

We are now ready to prove Theorem 3.1 by induction over L .

For L = 1 we need to prove that

$$\text{For } L = 1 \text{ we need to prove that} \\ & \psi _ { N _ { 1 } } ( x ) = \Psi _ { N _ { 1 } } ( x ) \left ( \prod _ { j = 1 } ^ { N _ { 1 } } z _ { j } ^ { r _ { 1 } } \right ), \quad d _ { N _ { 1 } } = N _ { 1 } r _ { 1 }, \\ & E _ { N _ { 1 } } = \sum _ { j = 1 } ^ { N _ { 1 } } ( r _ { 1 } + \frac { 1 } { 2 } g ( N _ { 1 } + 1 - 2 j ) \right ) ^ { 2 } + g ^ { 2 } N _ { 1 } ( N _ { 1 } - 1 ) \left ( \frac { \eta _ { 1 } } { \pi } - \frac { 1 } { 1 2 } \right ) \\ \text{solve } ( 3 0 ) \text{ for } N = N _ { 1 }. \text{ For } r _ { 1 } = 0 \text{ this follows from Lemma } 4. 1 \text{ in the special case } N = N _ { 1 }, \\ \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{ $\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_} \text{$\_}$$

solve (30) for N = N 1 . For r 1 = 0 this follows from Lemma 4.1 in the special case N = N 1 , M = 0 (this is non-trivial only in case N 1 = k &gt; 1). The general case of non-zero r 1 follows from this by a simple computation explained in Section 4.2; see (31) (note that ∏ j z r 1 j = exp(i r 1 ∑ j x j ), and that E N 1 in (34b) and (50) agree).

For L &gt; 1 one can construct the function ψ r,s,L ( x , . . . , x 1 n ; τ ) in Theorem 3.1 inductively using the integral operator in (45)-(46) as follows: ψ λ ( x ) = ψ N L ( x ) with

$$\psi _ { N _ { a } } ( x ) = ( \mathcal { K } _ { N _ { a } N _ { a - 1 } } ^ { Q _ { a } Q _ { a } ^ { \prime } } \psi _ { N _ { a - 1 } } ) ( x ), \ \ Q _ { a } = r _ { a } - N _ { a - 1 } \frac { g } { 2 }, \ \ Q _ { a } ^ { \prime } = r _ { a } - N _ { a } \frac { g } { 2 }$$

for a = 2 3 , , . . . , L , with ψ N 1 ( x ) in (50). Then the repeated application of Lemma 4.2 guarantees that the function ψ N a ( x ) has the claimed properties.

We are left to show that d λ = d N L and E λ = E N L where

$$d _ { N _ { a } } = d _ { N _ { a - 1 } } + N _ { a } Q _ { a } - N _ { a - 1 } Q _ { a } ^ { \prime }, \quad & & ( 5 2 a ) \\. \, \_ \quad \dots \, \_ \dots \, \_ \dots \, \alpha \quad \_ \dots \, \quad \dots$$

$$E _ { N _ { a } } & = E _ { N _ { a - 1 } } + 2 ( Q _ { a } - Q _ { a } ^ { \prime } ) d _ { N _ { a - 1 } } + N _ { a - 1 } ( Q _ { a } ^ { \prime } ) ^ { 2 } - 2 N _ { a - 1 } Q _ { a } Q _ { a } ^ { \prime } + N _ { a } Q _ { a } ^ { 2 } + c _ { N _ { a } N _ { a - 1 } } \quad ( 5 2 b ) \\ \epsilon \dots \quad \circ \tau \quad \tau \quad \tau \colon \dots \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \end{matrix}$$

for a = 2 3 , , . . . , L . This can be verified by straightforward computations (the interested reader can find details in Appendix C.3).

## 5 Concluding remarks

We give some outlook and conjectures.

- 1. As explained in more detail below, results in [7, 9, 23, 24, 25, 26] suggest to us that, for arbitrary n ∈ Z ≥ 1 , g ∈ C , κ ∈ C \ { 0 } and λ ∈ C n , there exist functions P (1 /g ) λ κ,n ; ( z ; p ) such that Ψ n ( z ; p P ) (1 /g ) λ κ,n ; ( z ; p ) is a solution of the non-stationary eCS equation in (18) with eigenvalues as on the RHS in (22), and

$$\sum _ { j = 1 } ^ { n } z _ { j } \frac { \partial } { \partial z _ { j } } P _ { \lambda ; \kappa, n } ^ { ( 1 / g ) } ( z ; p ) = | \lambda | P _ { \lambda ; \kappa, n } ^ { ( 1 / g ) } ( z ; p ), \quad | \lambda | \equiv \sum _ { j = 1 } ^ { n } \lambda _ { j }. \\ \intertext { a t i m m o n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d I t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d e t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d i t i n d \nolimits }$$

These functions are analytic functions of ( z, p ) and ( λ, κ ) in suitable regions ⊂ C n +1 and, for partitions λ , they are symmetric functions in z that reduce to the Jack polynomials

P (1 /g ) λ,n ( z ) in the limit p → 0. We believe that there is a beautiful mathematical theory of this two-parameter deformation of the Jack polynomials which, at this point, is only partly developed. The result in the present paper is a simple special case or this.

We now discuss in more detail the reasons for these expectations.

- ( ) i It seems that the affine analogue of the Jack polynomials defined and studied in [23, 24] are special cases of these functions P (1 /g ) λ κ,n ; ( z ; p ).
- ( ii ) We constructed series representations of such functions P (1 /g ) λ κ,n ; ( z ; p ) for the simplest non-trivial case n = 2 in [9]. 11

/negationslash

- ( iii ) A construction of singular series solutions (in the sense of [25]) of the non-stationary eCS equation for κ = 0 and arbitrary n ∈ Z ≥ 2 by one of us in [26], which is straightforward to extend this construction to κ = 0. In fact, this construction is easier for non-zero κ [9]: first, if κ = 0, one need not compute the eigenvalue E τ ( ) in (18) and, second, if /Ifractur ( κ ) = 0, it is easy to prove absolute convergence of these singular series solutions in a suitable domain ⊂ C n +1 using the method in [26]. Moreover, using the kernel functions in (38), one can transform these singular solutions to such described in the previous paragraph (this is explained in [25]), but the latter is restricted to κ = kg with k ∈ Z and, for k &gt; 0, these solutions are incomplete (for k &lt; 0 they are over-complete).

/negationslash

/negationslash

- ( iv ) Our main reason are Shiraishi's non-stationary eCS functions f eCS ( z ; p λ κ g | ; | ) in [7, Section 7] mentioned already in the introduction. As Shiraishi conjectured, these functions are such that 12

$$\begin{matrix} z _ { 1 } ^ { \lambda _ { 1 } } \cdots z _ { n } ^ { \lambda _ { n } } f ^ { e c s } ( z ; p | \lambda ; - \kappa | g ) \\.......... = ( 1 / a \lambda........................................................................................................................................................................................................ \end{matrix}$$

have exactly the properties of the functions P (1 /g ) λ κ,n ; ( z ; p ) described in the first paragraph of this section.

- 2. It would be interesting to prove that the elliptic generalizations of the Jack polynomials constructed in the present paper are identical with the functions in (53) in special cases. Since both this functions are given by explicit formulas, this non-trivial test of Shiraishi's conjecture in [7, Section 7] seems feasible.
- 3. It is interesting to speculate what the generalizations of the properties in (6)-(10) of the Jack polynomials to the functions P (1 /g ) λ κ,n ; ( z ; p ) described in Remark 1 above are. We believe that the result in the present paper suggest the following (we stress that the statements in the following three equations are conjectures).

First, the natural elliptic generalization of the scalar product in (5) is

/negationslash

$$\cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \\ \langle f, g \rangle _ { n } ^ { \prime } \equiv \frac { 1 } { n! } \oint _ { | z _ { 1 } | = 1 } \frac { d z _ { 1 } } { 2 \pi i z _ { 1 } } \cdots \oint _ { | z _ { n } | = 1 } \frac { d z _ { n } } { 2 \pi i z _ { n } } \left ( \prod _ { j \neq k } ^ { n } \theta ( 1 - z _ { j } / z _ { k } ; p ) \right ) ^ { g } f ( z ; p ) \overline { g ( z ; \bar { p } ) }, \\ \overline { 1 1 } \text{The routine in } [ 0 ] \, \text{or automatically for the non-attention} \, \text{from anation} \, \text{ but this method as anocial case}$$

11 The results in [9] are actually for the non-stationary Heun equation, but this includes as special case the non-stationary Lam´ e equation which, essentially, is the non-stationary eCS equation for n = 2.

12 Note that our z , κ and g correspond to x , -k and β in [7, Section 7], respectively.

and we expect that the natural generalization of the orthogonality relations in (6) is

$$\left \langle P _ { \lambda ; \kappa, n } ^ { ( 1 / g ) }, P _ { \mu ; - \bar { \kappa }, n } ^ { ( 1 / g ) } \right \rangle _ { n } ^ { \prime } = \delta _ { \lambda \mu } \mathcal { N } _ { \lambda ; \kappa, n } ( p, g ) \\ \text{that may denied on $n$ and $\kappa$ (this is suggested} )$$

for some N λ κ,n ; ( p, g ) that may depend on p and κ (this is suggested to us by the quantum mechanical interpretation of the model, using some quantum mechanics folklore). Second, we expect the following generalization of the generating function in (8) holds true,

$$\prod _ { j = 1 } ^ { n } \prod _ { k = 1 } ^ { m } \frac { 1 } { \theta ( z _ { j } / \xi _ { k } ; p ) ^ { g } } & = \sum _ { \lambda } \frac { 1 } { b _ { \lambda } ; n, m ( p, g ) } P _ { \lambda ; ( n - m ) g, n } ^ { ( 1 / g ) } ( z ; p ) P _ { \lambda ; ( m - n ) g, m } ^ { ( 1 / g ) } ( 1 / \xi ; p ) \quad ( m \leq n ) \\, \quad, \quad \cdot \quad, \cdot \quad, \cdot \quad, \cdot \quad, \cdot \quad, \cdot \quad, \cdot \quad, \cdot \quad, \cdot \quad,$$

where the sum is over all integer vectors λ ∈ Z m , and the constants b λ n,m ; ( p, g ) may depend on n m , and p (this, it seems to us, is the unique equation consistent with all known facts). Third, it is natural to expect that the Pieri relation in (10) remains unchanged,

$$( z _ { 1 } \cdots z _ { n } ) ^ { r } P _ { \lambda ; \kappa, n } ^ { ( 1 / g ) } ( z ; p ) = P _ { \lambda + ( r ^ { n } ) ; \kappa, n } ^ { ( 1 / g ) } ( z ; p ) \ \ ( \lambda \in \mathbf P _ { n }, r \in \mathbb { Z } ).$$

It would be interesting to prove these relations since they would allow to prove the results in the present paper by the very same argument used in Appendix B for p = 0.

- 4. As discussed in Remark 3. above, we believe that the functions P (1 /g ) λ κ,,n ; are an orthogonal system with respect to a natural scalar product 〈· , ·〉 ′ n exactly if p is real and κ purely imaginary, ¯ = p p and ¯ = κ -κ (we assume g &gt; 0). This is exactly the case when the equation in (1) has a natural physical interpretation as non-stationary Schr¨dinger equation with o 〈· , ·〉 ′ n the standard L 2 -scalar product playing an important role in the quantum mechanical interpretation of this equation. To us, this is strong reason to think about the equation in (1) as Schr¨dinger equation (rather than heat equation), even though the results in the o present paper are for real κ values.
- 5. All known CMS-like systems allow for deformations in the sense of [27], and the nonstationary eCS equation is no exception: its deformation is given by

$$\left ( \kappa \frac { \mathrm i } { \pi } \frac { \partial } { \partial \tau } + H _ { n, \tilde { n } } ( x, \tilde { x } ) \, \right ) \psi ( x, \tilde { x } ) = E \psi ( x, \tilde { x } ) \\ \cdot \, \sim \, \sim \, \Re \quad \Phi \quad \Phi.$$

with the deformed eCS Hamiltonian

$$H _ { n, \tilde { n } } ( x, \tilde { x } ) \equiv - \sum _ { j = 1 } ^ { n } \frac { \partial ^ { 2 } } { \partial x _ { j } ^ { 2 } } + \sum _ { j < j ^ { \prime } } ^ { n } 2 g ( g - 1 ) V ( x _ { j } - x _ { j ^ { \prime } } ) \\ + g \sum _ { k = 1 } ^ { \tilde { n } } \frac { \partial ^ { 2 } } { \partial \tilde { x } _ { k } ^ { 2 } } + \sum _ { k < k ^ { \prime } } ^ { \tilde { n } } \frac { 2 ( g - 1 ) } { g } V ( \tilde { x } _ { k } - \tilde { x } _ { k ^ { \prime } } ) + \sum _ { j = 1 } ^ { n } \sum _ { k = 1 } ^ { \tilde { n } } 2 ( 1 - g ) V ( x _ { j } - \tilde { x } _ { k } ) \\ \text{ depending to two arbitrary variable numbers } n \text{ and } \tilde { n } \text{ (we suppress the } \tau \text{-dependence).}$$

depending to two arbitrary variable numbers n and ˜ (we suppress the n τ -dependence). This operator is known to define a quantum integrable system (this is proved in [28] for the case ˜ = 1, but the result is expected to be true for arbitrary ˜). n n Moreover, generalized kernel functions K N,N,M,M ˜ ˜ ( x, x, y, y ˜ ˜) for this non-stationary deformed eCS equation are known [29]: they obey relations as in (28) but for deformed eCS operators H N,N ˜ ( x, x ˜) and

H M,M ˜ ( y, y ˜), and κ = ( N -M g ) is replaced by ( N -M g ) -˜ + N ˜ . M We expect that, using these generalized kernel functions, the construction in the present paper can be generalized to find solutions of the deformed non-stationary eCS equation. We leave this for future work.

It is possible to generalize many results about the Jack polynomials in a natural way to the deformed case [30, 31, 32]. We thus believe that there is a natural generalization of the theory described in Remark 1. above to the deformed case as well.

6. We believe that Theorem 3.1 holds true for complex g with real part /Rfractur ( g ) &gt; 0. It would be interesting to prove this.

## A Special functions

For the convenience of the reader we collect the well-known definitions and properties of special functions that we use (we follow Whittaker &amp; Watson [22]).

The Weierstrass ℘ -function with periods (2 ω , 1 2 ω 2 ) is defined as

$$\wp ( x | \omega _ { 1 }, \omega _ { 2 } ) \equiv \frac { 1 } { x ^ { 2 } } + \sum _ { ( m, n ) \in \mathbb { Z } \times \mathbb { Z } \langle ( 0, 0 ) \} } \left ( \frac { 1 } { ( x - \Omega _ { m, n } ) ^ { 2 } } - \frac { 1 } { \Omega _ { m, n } ^ { 2 } } \right ) \\ \mu \quad \ \. \quad \. \quad \. \quad \. \quad \. \quad \.$$

with Ω m,n ≡ 2 mω 1 +2 nω 2 and /Ifractur ( ω /ω 2 1 ) &gt; 0. Moreover,

$$\vartheta _ { 1 } ( x | \tau ) \equiv 2 \sum _ { n = 0 } ^ { \infty } ( - 1 ) ^ { n } \text{e} ^ { \pi i \tau ( n + 1 / 2 ) ^ { 2 } } \sin \left ( ( 2 n + 1 ) x \right ) \\. \ e d { \dots } \ e d { \dots } \ e d { \dots } \ e d { \dots } \ e d { \dots } \ e d { \dots } \ e d { \dots } \ e d { \dots } \ e d { \dots } \ e d { \dots } \ e d { \dots } \ e d { \dots } \ e d { \dots } \ e d { \dots }$$

with τ ≡ ω /ω 2 1 defines the corresponding Jacobi theta function (we set ω 1 = π in the main text).

We quote two well-known results that we use [22]. First,

$$\wp ( z | \omega _ { 1 }, \omega _ { 2 } ) = - \frac { \eta _ { 1 } } { \omega _ { 1 } } + \frac { \pi ^ { 2 } } { 4 \omega _ { 1 } ^ { 2 } } \sum _ { n \in \mathbb { Z } } \frac { 1 } { \sin ^ { 2 } ( \frac { \pi ( z + 2 n \omega _ { 2 } ) } { 2 \omega _ { 1 } } ) }$$

with and second,

$$\frac { \eta _ { 1 } } { \omega _ { 1 } } = \frac { \pi ^ { 2 } } { \omega _ { 1 } ^ { 2 } } \left ( \frac { 1 } { 1 2 } + \sum _ { n = 1 } ^ { \infty } \frac { 1 } { 2 \sin ^ { 2 } ( \frac { n \pi \omega _ { 2 } } { \omega _ { 1 } } ) } \right ),$$

$$\vartheta _ { 1 } ( x | \tau ) = 2 p ^ { 1 / 8 } G \sin ( x ) \prod _ { m = 1 } ^ { \infty } ( 1 - 2 p ^ { m } \cos ( 2 x ) + p ^ { 2 m } ) \\ \prod ^ { \infty } \,. ( 1 - n ^ { m } ) \text{ and } n = e ^ { 2 \pi i \tau } \text{ (note that our $n$ is $a^{2}$ in [22]}).\text{ These res}$$

with G ≡ ∏ ∞ m =1 (1 -p m ) and p = e 2 π τ i (note that our p is q 2 in [22]). These results imply the identities in (16)-(17), and that the function ϑ x ( ) in (36) equals ϑ 1 ( 1 2 x τ | ) /p 1 8 / G .

## B Proof of Proposition 2.2

We give a concise proof of Proposition 2.2. It is similar to the proof in [5, Section 6], but we emphasize technical details allowing to use circles as integration contours even for non-integer g , different from [5] (see Remark 1.1).

The identity in (8) is in the sense of formal power series. However, we need information about convergence of the series on the RHS in (8). To get this we observe that the LHS in (8) can be expanded in a Taylor series in the variables z j and 1 /ξ k using the binomial theorem repeatedly, and this series is absolutely and uniformly convergent in the region

$$| z _ { j } / \xi _ { j } | < \epsilon \ \forall j, k,$$

for arbitrary /epsilon1 in the range 0 &lt; /epsilon1 &lt; 1. Thus, in such a region, we can re-sum this Taylor series so as to get the series on the RHS in (8), and the interchange of summation and integration below is justified.

The second analytic fact we need is that the radius of the integration contours in the scalar product in (5) can be changed from 1 to an arbitrary value ρ &gt; 0 without changing the result, i.e.,

/negationslash

$$\langle f, g \rangle _ { n } ^ { \prime } \equiv \$$

is independent of ρ &gt; 0 (a concise proof is given at the end of this Appendix for the convenience of the reader, and more detailed discussion can be found in, e.g., [32]). Moreover, for g ∈ Λ n with real coefficients, g (1 /z ¯) = g (1 /z ), and the Jack polynomials have this property.

To prove Proposition 2.2 we use the following identity implied by (8) and (10),

$$\begin{array} {$$

(the interpretation of this for r &lt; 0 is explained in Remark 2.1). For given Jack polynomial P (1 /g ) µ,m ( ξ ), we fix ρ &gt; 0 and assume | z j | &lt; ρ for all j to compute

/negationslash

/negationslash

$$P _ { \mu, m } ^ { ( 1 / g ) } ( \xi ), \$$

$$= \frac { m! \mathcal { N } _ { \mu, m } ( g )$$

In this computation we interchanged summation and integration (first equality), used (54) and that the Jack polynomials have real coefficients (second equality), and inserted the orthogonality relations in (6) (third equality).

It is important to note that, if r &gt; µ m , then 〈 P (1 /g ) µ,m , P (1 /g ) ν +( r m ) ,m 〉 ′ m = 0 for all ν ∈ P m , and thus we have to set P (1 /g ) µ +( r n ) -( r m ) ,n ( z ) ≡ 0 in such a case. Note that

$$\lambda \equiv \mu + ( r ^ { n } ) - ( r ^ { m }$$

so this latter case exactly corresponds to λ = ( λ , . . . , λ 1 n ) not satisfying λ 1 ≥ · · · ≥ λ n . Thus the result of this computation is an integral transform mapping a Jack polynomial P (1 /g ) µ,m to a Jack polynomial P (1 /g ) λ,n provided that λ satisfies λ 1 ≥ · · · ≥ λ n , and the integral is 0 otherwise. Note that our derivation shows that the integral is only well-defined if | z j | &lt; ρ for all j .

We now apply this result to the family of integer vectors λ ( a ) in (11), specializing to ( n, m, r ) = ( N a +1 , N a , r a +1 ) and µ = λ ( a ) . Since λ ( a ) +( r N a +1 a +1 ) -( r N a a +1 ) = λ ( a +1) we get

/negationslash for | z j | &lt; ρ ( j = 1 , . . . , n ). The result in Proposition 2.2 can be obtained from this by iteration, starting out with a = 1 and

$$\begin{array} { c } \vee \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \end{array}$$

$$P _ { \lambda ^ { ( 1 ) }, N _ { 1 } } ^ { ( 1 / g ) } ( z ) = ( z _ { 1 } \cdots z _ { N _ { 1 } } ) ^ { r _ { 1 } } \ \ ( N _ { 1 } = s _ { 1 }, \, \lambda ^ { ( 1 ) } = ( r _ { 1 } ^ { s _ { 1 } } ) ),$$

inserting the result into (55) for a = 2, etc., up to a = L -1.

To keep this appendix self-contained we conclude with a proof that the RHS in (54) is independent of ρ &gt; 0 and equal to the scalar product defined in (5).

Proof of (54) . We denote the RHS in (54) as ( f, g ) ρ . Note that ( f, g ) 1 is the RHS of (5).

To prove the result it is enough to show it for monomial symmetric polynomials, f = m λ and g = m µ with arbitrary λ, µ ∈ P n . Use polar coordinates z j = ρ e i x j to show that

$$( m _ { \lambda }, m _ { \mu } ) _ { \rho } = \rho ^ { | \lambda | - | \mu | } ( m _ { \lambda }, m _ { \mu } ) _ { 1 } = \rho ^ { | \lambda | - | \mu | } \langle m _ { \lambda }, m _ { \mu } \rangle _ { n } ^ { \prime } \quad ( \rho > 0 )$$

with | λ | = ∑ j λ j and similarly for µ . Since the Euler operator D n ≡ ∑ n j =1 z j ∂z ∂ j is selfadjoint with respect to the scalar product 〈· , ·〉 ′ n , and D m n λ = | λ m | λ , 〈 m ,m λ µ 〉 ′ n = 0 unless | λ | - | µ | = 0. This implies the result.

## C Details of the proof in Section 4

## C.1 Analyticity

We prove that the integral transform in (46) defines an analytic function in the specified domain.

We note that the integral in (46) has the form

/negationslash

$$F ( z ; p ) \equiv \left ( \prod _ { j = 1 } ^ { M } \oint _ { | \xi _ { j } | = \rho } \frac { d \xi _ { j } } { 2 \pi \mathrm i \xi _ { j } } \right ) \left ( \prod _ { j \neq k } ^ { M } \theta ( \xi _ { j } / \xi _ { k } ; p ) \right ) ^ { g } f ( z, \xi ; p ) \\., \ e \quad.$$

with the function

Moreover, f ( z, ξ ; p ) is analytic in the following domain ⊂ C N + M +1 ,

$$f ( z, \xi ; p ) = \frac { \left ( \prod _ { j = 1 } ^ { N } z _ { j } ^ { r } \right ) \left ( \prod _ { j = 1 } ^ { M } \xi _ { k } ^ { - r ^ { \prime } } \right ) } { \prod _ { j = 1 } ^ { N } \prod _ { k = 1 } ^ { M } \theta \left ( z _ { j } / \xi _ { k } ; p \right ) ^ { g } } \mathcal { P } _ { M } ( \xi ; p ). \\ \xi ; p ) \, \text{is analytic in the following domain} \subset \mathbb { C } ^ { N + M + 1 },$$

$$| p | < | z _ { j } / \xi _ { k } | < 1 \text{ and } \rho _ { 0 } | p | < | \xi _ { k } | < \rho _ { 0 } \text{ } ( j = 1, \dots, N, \, k = 1, \dots, M ), \text{ } | p | < 1 \text{ } ( 5 7 )$$

with ρ 0 &gt; 0 specified in Lemma 4.2 and ρ 0 | p | &lt; ρ &lt; ρ 0 . Indeed, 1 /θ z /ξ ( j k ) g is analytic in the region | p | &lt; | z /ξ j k | &lt; 1 for all j, k , and the analyticity of the other factors is guaranteed by the assumptions on P M ( ξ ; p ) and r, r ′ in Lemma 4.2. We now claim : For every function f ( z, ξ ; p ) that is analytic in the region (57) , the integral in (56) defines a function F z p ( ; ) that is analytic in the region (47) . This implies the analyticity result for the integral transform in (46) stated in Lemma 4.2.

We are left to prove this claim. Since

/negationslash

$$\left ( \prod _ { j \neq k } ^ { M } \theta ( \xi _ { j } / \xi _ { k } ; p ) \right ) ^ { g } = \prod _ { j < k } ^ { M } \prod _ { m \in \mathbb { Z } } \left [ \left ( 1 - p ^ { | m | } \frac { \xi _ { j } } { \xi _ { k } } \right ) \left ( 1 - p ^ { | m | } \frac { \xi _ { k } } { \xi _ { j } } \right ) \right ] ^ { g } \\ \sim \, \nu \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \end{array}$$

and | ξ /ξ j k | = 1 on the integration contours, it is clear that, for fixed ( z, p ) in the region (47), the integrand in (56) is bounded on the integration contours, and thus F z p ( ; ) is well-defined. Moreover, F z p ( ; ) can be differentiated with respect to each variable z j ( j = 1 , . . . , N ) and p separately by interchanging differentiation with integration, which is justified by the Leibniz integral rule. Thus, F z p ( ; ) is analytic in each variable separately as long as the conditions in (47) are fulfilled and, by Hartogs's theorem, this implies that F z p ( ; ) is analytic in the region (47).

## C.2 Boundary terms

We give details on how to prove the identities in (48) and (49).

Recalling the definition in (15) and using the Leibniz integration rule we write the LHS in (48) as an integral of total derivatives,

$$\dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \\ - \mathfrak { i } \left ( \prod _ { j = 1 } ^ { M } \int _ { - \pi - \mathfrak { i } \epsilon } ^ { \pi - \mathfrak { i } \epsilon } \frac { d y _ { j } } { 2 \pi } \right ) \sum _ { k = 1 } ^ { M } \frac { \partial } { \partial y _ { k } } K _ { N M } ( x, y ) \mathrm e ^ { - \mathrm i Q ^ { \prime } | y | } \psi _ { M } ( y ). \\ \text{definitions we find that this is proportional to (we ignore irrelevant $u$-inc} }$$

Inserting definitions we find that this is proportional to (we ignore irrelevant y -independent factors)

$$& \left ( \prod _ { j = 1 } ^ { M } \int _ { - \pi - i \epsilon } ^ { \pi - i \epsilon } \frac { d y _ { j } } { 2 \pi } \right ) \sum _ { k = 1 } ^ { M } \frac { \partial } { \partial y _ { k } } F ( y ), \quad F ( y ) \equiv \frac { \Psi _ { M } ( y ) ^ { 2 } } { \prod _ { j = 1 } ^ { N } \prod _ { k = 1 } ^ { M } \theta ( e ^ { - i y _ { k } } z _ { j } / \rho ) ^ { g } } e ^ { - i \nu ^ { \prime } | y | } \mathcal { P } _ { M } ( \rho e ^ { i y } ) \\ & \text{where $\rho^{i y}$ is short for $(\rho^{i y _ { 1 }},\dots,\rho^{i y _ { M }} ). $ All factors but $\Psi_{M}(y)^{2}$ in the function $F(y)$ are} \\ & \text{obvious} v \, C ^ { \infty } \text{ and $2\pi$-periodic in each variable $u$ for $u\in[-\pi\pi]^{M}$ (recall our analyticiv}$$

where ρ e i y is short for ( ρ e i y 1 , . . . , ρ e i y M ). All factors but Ψ M ( y ) 2 in the function F y ( ) are obviously C ∞ and 2 π -periodic in each variable y j for y ∈ -[ π, π ] M (recall our analyticity assumptions on P M , and that r ′ ∈ Z and | p | &lt; z | j | /ρ &lt; 1). The only potentially problematic factor that could spoil that the total derivatives of F y ( ) integrate to 0 is Ψ M ( y ) 2 for noninteger g .

Using definitions we can write

$$\overset { \smile } { \Psi } _ { M } ( y ) ^ { 2 } = \left ( \prod _ { j < k } ^ { M } \sin ^ { 2 } \left ( \frac { y _ { j } - y _ { k } } { 2 } \right ) \right ) ^ { g } \left ( \prod _ { j < k } ^ { N } \prod _ { m = 1 } ^ { \infty } ( 1 - 2 p ^ { m } \cos ( y _ { j } - y _ { k } ) + p ^ { 2 m } ) \right ) ^ { 2 g }, \\ \text{which makes manifest that. actually, the only potentially troblematic factor is}$$

which makes manifest that, actually, the only potentially problematic factor is

$$\max \max \max _ { J }, \, \max \max _ { J }, \, \max \max _ { J }, \, \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \min \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \max \text{$\mathrm \Sigma$} } \, \quad \quad \quad$$

(recall that | p | &lt; 1). We now show that Ψ (0) M ( y ) 2 is C 1 and 2 π -periodic in each variable y j for y ∈ -[ π, π ] M . This implies the same properties for the function F y ( ), which proves (48).

The 2 π -periodicity of Ψ (0) M ( y ) 2 is obvious. To prove the C 1 -property we compute the partial derivative ∂ ∂y j Ψ (0) M ( y ) 2 and find that it is a finite sum of terms proportional to

/negationslash

This behaves like ( y j -y k ) -1+2 g as y j → y k , and it is obviously is continuous everywhere else. Since we assume g &gt; 1 / 2, this proves that Ψ (0) M ( y ) 2 is C 1 . This concludes the proof of (48).

$$\cot \left ( \frac { y _ { j } - y _ { k } } { 2 } \right ) \Psi _ { M } ^ { ( 0 ) } ( y ) ^ { 2 } \quad ( k \neq j ). \\. \ e. \ \lambda _ { 1 } + 2 g \ e c \ e. \ \rightharpoonup \ e. \ \text{and it is otherwise}$$

The proof of (49) is similar to the proof of (48), and we are therefore sketchy. Inserting definitions we find that the LHS is an integral of total derivatives,

$$\text{emions we in an image} \, \text{in} \, \text{in} \, \text{is an image} \, \text{or too} \, \text{varivatives}, \\ \left ( \prod _ { j = 1 } ^ { M } \int _ { - \pi - i \epsilon } ^ { \pi - i \epsilon } \frac { d y _ { j } } { 2 \pi } \right ) \sum _ { k = 1 } ^ { M } \frac { \partial } { \partial y _ { k } } \left ( - \left \{ \frac { \partial } { \partial y _ { k } } K _ { N M } ( x, y ) \right \} \, \text{e} ^ { - i Q ^ { \prime } | y | } \psi _ { M } ( y ) \\ & + K _ { N M } ( x, y ) \frac { \partial } { \partial y _ { k } } \epsilon ^ { - i Q ^ { \prime } | y | } \psi _ { M } ( y ) \right ).$$

Computing the terms one finds again that the only potentially problematic term comes from the factor Ψ (0) M ( y ) 2 in the integrand. As in the proof of (48) above one checks that this factor, actually, is not problematic, either. The only difference to (48) is that we now have second order derivative terms ∂ 2 ∂y 2 j Ψ (0) M ( y ) 2 which give rise to terms (we ignore less singular terms)

/negationslash which behaves like ( y j -y k ) -2+2 g as y j → y k (continuity otherwise is again obvious). While these terms are singular for 1 / 2 ≤ g &lt; 1 (they are regular for g ≥ 1), the singularities are integrable for g &gt; 1 / 2, and this is enough to conclude that (49) holds true.

$$\frac { 1 } { \sin ^ { 2 } \left ( \frac { y _ { j } - y _ { k } } { 2 } \right ) } \Psi _ { M } ^ { ( 0 ) } ( y ) ^ { 2 } \quad ( k \neq j ), \\ \intertext { for $1/2<q<1$ (they are regular for $q>$} } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{\ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \$$

## C.3 Eigenvalue computations

We give some details of computations that verify (52a) and (52b).

We start with (52a). We insert (23) for the cases n = N a and n = N a -1 to compute

$$d _ { N _ { a } } - d _ { N _ { a - 1 } } = \sum _ { j = N _ { a - 1 } + 1 } ^ { N _ { a } } \lambda _ { j } = k r _ { a } \\ + 1 \, < \, i \, < \, N \, \text{ on } N \, - \, N \, + 1 \, + 1 \, \text{ to see all } 1 \, + 1$$

using λ j = r a for N a -1 +1 ≤ j ≤ N a and N a = N a -1 + k (recall that we assume that λ is as in (20)). From (51) it follows that this is equal to N Q a a -N a -1 Q ′ a , which proves (52a).

To prove (52b) we use (22) to compute, in a similar manner,

$$1 o \text{ prove } ( \partial 2 D ) \text{ we use } ( Z Z ) \text{ to compute, in a similar manner,} \\ E _ { N _ { a } } - E _ { N _ { a - 1 } } = \sum _ { j = N _ { a - 1 } + 1 } ^ { N _ { a } } \left ( r _ { a } + \frac { 1 } { 2 } g ( N _ { a } + 1 - 2 j ) \right ) ^ { 2 } + \sum _ { j = 1 } ^ { N _ { a - 1 } } \left ( g k \left ( \lambda _ { j } + \frac { 1 } { 2 } g ( N _ { a - 1 } + 1 - 2 j ) \right ) \\ + \frac { 1 } { 4 } ( g k ) ^ { 2 } \right ) + g ^ { 2 } \left [ N _ { a } ( N _ { a } - 1 ) - N _ { a - 1 } ( N _ { a - 1 } - 1 ) \right ] \left ( \frac { \eta _ { 1 } } { \pi } - \frac { 1 } { 1 2 } \right ) \\ = k r _ { a } ^ { 2 } - g k r _ { a } N _ { a - 1 } + \frac { 1 } { 1 2 } k g ^ { 2 } \left ( 3 N _ { a - 1 } ^ { 2 } + 3 N _ { a - 1 } k + k ^ { 2 } - 1 \right ) \\ + g k d _ { N _ { a - 1 } } + g ^ { 2 } k ( 2 N _ { a - 1 } + k - 1 ) \left ( \frac { \eta _ { 1 } } { \pi } - \frac { 1 } { 1 2 } \right ),$$

$$\mu _ { 1 2 } & \lambda _ { 1 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \\ \mu _ { 1 2 } & \lambda _ { 1 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 } \lambda _ { 2 }$$

recalling (23) in the last step. Using (23), (39) and (51) one can check that this equals 2( Q a -Q ′ a ) d N a -1 + N a -1 ( Q ′ a ) 2 -2 N a -1 Q Q a ′ a + N Q a 2 a + c N N a a -1 , which proves (52b).

## References

- [1] F. Calogero, Solution of the one-dimensional N-body problems with quadratic and/or inversely quadratic pair potentials , J. Math. Phys. 12 , 419-436 (1971).
- [2] B. Sutherland, Exact results for a quantum many body problem in one-dimension. II. , Phys. Rev. A5 , 1372-1376 (1972).

| [3]   | M.A. Olshanetsky and A.M. Perelomov, Quantum completely integrable systems con- nected with semisimple Lie algebras , Lett. Math. Phys. 2 , 7-13 (1977).                                                                                            |
|-------|-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [4]   | H. Jack, A class of symmetric polynomials with a parameter , Proc. Roy. Soc. Edinburgh Sect. A 69 , 1-18 (1970).                                                                                                                                    |
| [5]   | H. Awata, Y. Matsuo, S. Odake, J. Shiraishi, Excited states of the Calogero-Sutherland model and singular vectors of the W N algebra , Nucl. Phys. B 449 , 347-374 (1995), arXiv:hep-th/9503043.                                                    |
| [6]   | K. Mimachi and Y. Yamada, Singular vectors of the Virasoro algebra in terms of Jack symmetric polynomials , Comm. Math. Phys. 174 , 447-455 (1995),                                                                                                 |
| [7]   | J. Shiraishi, Affine screening operators, affine Laumon spaces, and conjectures con- cerning non-stationary Ruijsenaars functions , arXiv:1903.07495 [math.QA].                                                                                     |
| [8]   | E. Langmann and K. Takemura, Source identity and kernel functions for Inozemtsev- type systems , J. Math. Phys. 53 , 082105 (2012), arXiv:1202.3544 [math-ph].                                                                                      |
| [9]   | F. Atai and E. Langmann, Series Solutions of the Non-Stationary Heun Equation , SIGMA 14 , Paper 011, 32 pp. (2018), arXiv:1609.02525 [math-ph].                                                                                                    |
| [10]  | V.A. Fateev, A.V. Litvinov, A. Neveu, E. Onofri, A differential equation for a four- point correlation function in Liouville field theory and elliptic four-point conformal blocks , J. Phys. A 42 (2009), 304011, 29 pp. arXiv:0902.1331 [hep-th]. |
| [11]  | E. Langmann, Remarkable identities related to the (quantum) elliptic Calogero- Sutherland model , J. Math. Phys. 47 , 022101 (2006), arXiv:math-ph/0406061.                                                                                         |
| [12]  | V.G. Knizhnik, A.B. Zamolodchikov, Current algebra and Wess-Zumino model in two dimensions , Nucl. Phys. B 247 , 83-103 (1984).                                                                                                                     |
| [13]  | D. Bernard, On the Wess-Zumino-Witten models on the torus , Nucl. Phys. B 303 , 77-93 (1988).                                                                                                                                                       |
| [14]  | G. Felder and A. Varchenko, Special functions, conformal blocks, Bethe ansatz and SL(3 , Z ), Phil. Trans. R. Soc. Lond. A 359 , 1365-1373 (2001), arXiv:math/0101136.                                                                              |
| [15]  | P. I. Etingof and A. A. Kirillov Jr., Representations of affine Lie algebras, parabolic differential equations, and Lam´ e functions , Duke Math. J. 74 , 585-614 (1994). arXiv:hep-th/9310083.                                                     |
| [16]  | G. Felder and A. Varchenko, Integral representation of solutions of the elliptic Knizhnik-Zamolodchikov-Bernard equations , Int. Math. Res. Notices 1995, 221-233 (1995), arXiv:hep-th/9502165.                                                     |
| [17]  | F. Falceto and K. Gaw¸dzki, e Unitarity of the Knizhnik-Zamolodchikov-Bernard con- nection and the Bethe ansatz for the elliptic Hitchin systems , Comm. Math. Phys. 183 , 267-290 (1997), arXiv:hep-th/9604094.                                    |

| [18]   | G. Kuroki and T. Takebe, Bosonization and integral representation of solutions of the Knizhnik-Zamolodchikov-Bernard equations , Comm. Math. Phys. 204 , 587-618 (1999), arXiv:math/9809157.                                             |
|--------|------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [19]   | I.G. Macdonald, Symmetric functions and Hall polynomials (Second Edition), (Oxford University Press, New York, 1995).                                                                                                                    |
| [20]   | I.G. Macdonald, Commuting differential operators and zonal spherical functions , Al- gebraic groups Utrecht 1986, 189-200, Lecture Notes in Math. 1271 (Springer, Berlin, 1987).                                                         |
| [21]   | R.P. Stanley, Some combinatorial properties of Jack symmetric functions , Adv. Math. 77 , 76-115 (1989).                                                                                                                                 |
| [22]   | E.T. Whittaker and G.N. Watson, A course of modern analysis , Fourth Edition, Cam- bridge University Press (1940).                                                                                                                       |
| [23]   | P.I. Etingof and A.A. Kirillov, Jr., On the affine analogue of Jack and Macdonald polynomials , Duke Math. J. 78 , 229-256 (1995), arXiv:hep-th/9403168.                                                                                 |
| [24]   | P.I. Etingof, I.B. Frenkel, A.A. Kirillov Jr., Spherical functions on affine Lie groups , Duke Math. J. 80 , 59-90 (1995), arXiv:hep-th/9407047.                                                                                         |
| [25]   | E. Langmann, Singular eigenfunctions of Calogero-Sutherland type systems and how to transform them into regular ones , SIGMA 3 , 031, 18 pp. (2007), arXiv:math-ph/0702089.                                                              |
| [26]   | E. Langmann, Explicit solution of the (quantum) elliptic Calogero-Sutherland model , Ann. Henri Poincar´ e 15 , 755-791 (2014), arXiv:math-ph/0407050.                                                                                   |
| [27]   | O. Chalykh, M. Feigin, A.P. Veselov, New integrable generalizations of Calogero-Moser quantum problem , J. Math. Phys. 39 695-703 (1998).                                                                                                |
| [28]   | L.A. Khodarinova, Quantum integrability of the deformed elliptic Calogero-Moser prob- lem , J. Math. Phys. 46 , 033506, 22 pp. (2005), arXiv:math-ph/0406066.                                                                            |
| [29]   | E. Langmann, Source identity and kernel functions for elliptic Calogero-Sutherland type systems , Lett. Math. Phys. 94 , 63-75 (2010), arXiv:1003.0857 [math-ph].                                                                        |
| [30]   | A.N. Sergeev and A.P. Veselov, Generalised discriminants, deformed Calogero-Moser- Sutherland operators and super-Jack polynomials , Adv. Math. 192 341-375 (2005), arXiv:math-ph/0307036.                                               |
| [31]   | F. Atai and E. Langmann, Deformed Calogero-Sutherland model and fractional quan- tum Hall effect , J. Math. Phys. 58 011902, 27pp. (2017), arXiv:1603.06157 [math-ph].                                                                   |
| [32]   | F. Atai, M. Halln¨s, a E. Langmann, Orthogonality of super-Jack polynomials and a Hilbert space interpretation of deformed Calogero-Moser-Sutherland operators , Bull. Lond. Math. Soc. 51 , 353-370 (2019), arXiv:1802.02016 [math-ph]. |
Farrokh Atai Edwin Langmann
Exact solutions by integrals of the non-stationary elliptic Calogero-Sutherland equation
We use generalized kernel functions to construct explicit solutions by
integrals of the non-stationary Schr\"odinger equation for the Hamiltonian of
the elliptic Calogero-Sutherland model (also known as elliptic
Knizhnik-Zamolodchikov-Bernard equation). Our solutions provide integral
represenations of elliptic generalizations of the Jack polyomials.