Quasi-Exactly Solvable N-Body Spin Hamiltonians with Short-Range Interaction Potentials

We review some recent results on quasi-exactly solvable spin models presenting near-neighbors interactions. These systems can be understood as cyclic generalizations of the usual Calogero-Sutherland models. A nontrivial modification of the exchange operator formalism is used to obtain several infini...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Datum:2006
Hauptverfasser: Enciso, A., Finkel, F., González-López, A., Rodríguez, M.A.
Format: Artikel
Sprache:English
Veröffentlicht: Інститут математики НАН України 2006
Schriftenreihe:Symmetry, Integrability and Geometry: Methods and Applications
Online Zugang:http://dspace.nbuv.gov.ua/handle/123456789/146108
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Назва журналу:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Zitieren:Quasi-Exactly Solvable N-Body Spin Hamiltonians with Short-Range Interaction Potentials / A. Enciso, F. Finkel, A. González-López, M.A. Rodríguez // Symmetry, Integrability and Geometry: Methods and Applications. — 2006. — Т. 2. — Бібліогр.: 42 назв. — англ.

Institution

Digital Library of Periodicals of National Academy of Sciences of Ukraine
id irk-123456789-146108
record_format dspace
spelling irk-123456789-1461082019-02-08T01:23:39Z Quasi-Exactly Solvable N-Body Spin Hamiltonians with Short-Range Interaction Potentials Enciso, A. Finkel, F. González-López, A. Rodríguez, M.A. We review some recent results on quasi-exactly solvable spin models presenting near-neighbors interactions. These systems can be understood as cyclic generalizations of the usual Calogero-Sutherland models. A nontrivial modification of the exchange operator formalism is used to obtain several infinite families of eigenfunctions of these models in closed form. 2006 Article Quasi-Exactly Solvable N-Body Spin Hamiltonians with Short-Range Interaction Potentials / A. Enciso, F. Finkel, A. González-López, M.A. Rodríguez // Symmetry, Integrability and Geometry: Methods and Applications. — 2006. — Т. 2. — Бібліогр.: 42 назв. — англ. 1815-0659 2000 Mathematics Subject Classification: 81Q05; 35Q40 http://dspace.nbuv.gov.ua/handle/123456789/146108 en Symmetry, Integrability and Geometry: Methods and Applications Інститут математики НАН України
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection DSpace DC
language English
description We review some recent results on quasi-exactly solvable spin models presenting near-neighbors interactions. These systems can be understood as cyclic generalizations of the usual Calogero-Sutherland models. A nontrivial modification of the exchange operator formalism is used to obtain several infinite families of eigenfunctions of these models in closed form.
format Article
author Enciso, A.
Finkel, F.
González-López, A.
Rodríguez, M.A.
spellingShingle Enciso, A.
Finkel, F.
González-López, A.
Rodríguez, M.A.
Quasi-Exactly Solvable N-Body Spin Hamiltonians with Short-Range Interaction Potentials
Symmetry, Integrability and Geometry: Methods and Applications
author_facet Enciso, A.
Finkel, F.
González-López, A.
Rodríguez, M.A.
author_sort Enciso, A.
title Quasi-Exactly Solvable N-Body Spin Hamiltonians with Short-Range Interaction Potentials
title_short Quasi-Exactly Solvable N-Body Spin Hamiltonians with Short-Range Interaction Potentials
title_full Quasi-Exactly Solvable N-Body Spin Hamiltonians with Short-Range Interaction Potentials
title_fullStr Quasi-Exactly Solvable N-Body Spin Hamiltonians with Short-Range Interaction Potentials
title_full_unstemmed Quasi-Exactly Solvable N-Body Spin Hamiltonians with Short-Range Interaction Potentials
title_sort quasi-exactly solvable n-body spin hamiltonians with short-range interaction potentials
publisher Інститут математики НАН України
publishDate 2006
url http://dspace.nbuv.gov.ua/handle/123456789/146108
citation_txt Quasi-Exactly Solvable N-Body Spin Hamiltonians with Short-Range Interaction Potentials / A. Enciso, F. Finkel, A. González-López, M.A. Rodríguez // Symmetry, Integrability and Geometry: Methods and Applications. — 2006. — Т. 2. — Бібліогр.: 42 назв. — англ.
series Symmetry, Integrability and Geometry: Methods and Applications
work_keys_str_mv AT encisoa quasiexactlysolvablenbodyspinhamiltonianswithshortrangeinteractionpotentials
AT finkelf quasiexactlysolvablenbodyspinhamiltonianswithshortrangeinteractionpotentials
AT gonzalezlopeza quasiexactlysolvablenbodyspinhamiltonianswithshortrangeinteractionpotentials
AT rodriguezma quasiexactlysolvablenbodyspinhamiltonianswithshortrangeinteractionpotentials
first_indexed 2025-07-10T23:11:02Z
last_indexed 2025-07-10T23:11:02Z
_version_ 1837303410782633984
fulltext Symmetry, Integrability and Geometry: Methods and Applications Vol. 2 (2006), Paper 073, 11 pages Quasi-Exactly Solvable N -Body Spin Hamiltonians with Short-Range Interaction Potentials? A. ENCISO, F. FINKEL, A. GONZÁLEZ-LÓPEZ ∗ and M.A. RODRÍGUEZ Depto. F́ısica Teórica II, Universidad Complutense, 28040 Madrid, Spain ∗ E-mail: artemio@fis.ucm.es Received September 15, 2006, in final form October 23, 2006; Published online November 03, 2006 Original article is available at http://www.emis.de/journals/SIGMA/2006/Paper073/ Abstract. We review some recent results on quasi-exactly solvable spin models presenting near-neighbors interactions. These systems can be understood as cyclic generalizations of the usual Calogero–Sutherland models. A nontrivial modification of the exchange operator formalism is used to obtain several infinite families of eigenfunctions of these models in closed form. Key words: Calogero–Sutherland models; exchange operators; quasi-exact solvability 2000 Mathematics Subject Classification: 81Q05; 35Q40 1 Introduction In the early 1970s, F. Calogero [6] and B. Sutherland [37] introduced the quantum integrable systems that nowadays bear their names. Apart from their intrinsic mathematical interest [25, 21, 3], Calogero–Sutherland (CS) models have played a central role in Physics due to their relevant applications to as diverse topics as soliton theory [23, 32], quantum field and string theory [18, 9], quantum Hall effect [2], fractional statistics [28] and random matrix theory [36]. The first satisfactory explanation of the integrability of these models was given by Olshanetsky and Perelomov [26], who connected these models with the radial Laplacian of symmetric spaces associated to the root system AN . This unified view enabled them to introduce important generalizations, including different root systems and elliptic potentials. During the last decade CS models have been extended to the case of particles with internal degrees of freedom, which we shall henceforth call spin. There are two different approaches to spin CS models, namely the supersymmetric [5] and Dunkl operator [10, 29] formalisms. These methods have allowed to solve totally or partially several rational, trigonometric and elliptic spin models, both in their AN and BCN versions [8, 42, 16, 17]. The interest in spin CS models also increased as a consequence of their direct connection with the Haldane–Shastry (HS) spin chain [19, 34], which was laid bare by Polychronakos through the so called “freezing trick” [30]. This technique was also used in the construction of solvable spin chains associated to different potentials and root systems [31, 12]. Auberson, Jain and Khare [22, 1] introduced partially solvable versions of the CS models in which each particle only interacts with its nearest and next-to-nearest neighbors. Similar scalar models were also studied by Ezung, Gurappa, Khare and Panigrahi [15]. There are two reasons that make these kind of systems very promising from a physical point of view. First, some of them are related to the short-range Dyson model in random matrix theory [4]. Second, ?This paper is a contribution to the Proceedings of the Workshop on Geometric Aspects of Integ- rable Systems (July 17–19, 2006, University of Coimbra, Portugal). The full collection is available at http://www.emis.de/journals/SIGMA/Coimbra2006.html file:artemio@fis.ucm.es http://www.emis.de/journals/SIGMA/2006/Paper073/ http://www.emis.de/journals/SIGMA/Coimbra2006.html 2 A. Enciso, F. Finkel, A. González-López and M.A. Rodŕıguez the HS chains associated to these models would occupy an interesting intermediate position between the Heisenberg chain (short-range, position-independent interactions) and the usual HS chains (long-range, position-dependent interactions). A first step towards the construction of these chains was Deguchi and Ghosh’s definition of spin 1/2 versions of the Jain and Khare Hamiltonians [7] using the supersymmetric formalism. Unfortunately, all these authors solely managed to construct a few exact solutions, and all of them with trivial spin dependence, and the procedures developed to obtain exact solutions are by no means systematic. In [13, 14] we introduced three new families of spin near-neighbors models and used a nontriv- ial modification of the Dunkl operator method to obtain a wide range of fully explicit solutions. In this article we aim to review the main ideas underlying our constructions and explain their connection with the usual CS models. In Section 2 we define the Hamiltonians that we shall deal with and state our main result, which is a description of the algebraic states of the models. In Section 3 we review the philosophy underlying the calculation of exact solutions of Schrödinger’s equation by algebraic methods. In Section 4 we sketch the main logical steps that the construc- tion of the invariant flags rests upon. We finish the paper by showing in Section 5 how the actual computation of the algebraic eigenfunctions can be carried out. For the sake of brevity, complete proofs are not given in this review. 2 Main result Let Σ be the Hilbert space of internal degrees of freedom of N particles of spin M ∈ 1 2N. Let us fix a basis B = { |s1, . . . , sN 〉 : si ∈ {−M,−M + 1, . . . ,M} } of Σ and define the spin exchange operators Sij as Sij |s1, . . . , si, . . . , sj , . . . , sN 〉 = |s1, . . . , sj , . . . , si, . . . , sN 〉. These operators can be expressed in terms of the (normalized) fundamental generators of SU(2M + 1) as Sij = (2M + 1)−1 + ∑ a J a i J a j , the index sum ranging from 1 to 4M(M + 1). The Hamiltonians of the models we shall be concerned with are given by (cf. [14]) Hε = − ∑ i ∂2 xi + Vε, ε = 0, 1, 2, (1) where V0 = ω2r2 + ∑ i 2a2 (xi − xi−1)(xi − xi+1) + ∑ i 2a (xi − xi+1)2 (a− Si,i+1), (2a) V1 = ω2r2 + ∑ i b(b− 1) x2 i + ∑ i 8a2x2 i (x2 i − x2 i−1)(x 2 i − x2 i+1) + 4a ∑ i x2 i + x2 i+1 (x2 i − x2 i+1) 2 (a− Si,i+1), (2b) V2 = 2a2 ∑ i cot(xi − xi−1) cot(xi − xi+1) + 2a ∑ i csc2(xi − xi+1)(a− Si,i+1), (2c) with r2 = ∑ i x 2 i and a, b > 1/2. Here and in what follows, all sums and products run from 1 to N unless otherwise stated, with the identifications x0 ≡ xN and xN+1 ≡ x1. There is a hyperbolic potential analogous to (2c) that is recovered substituting xi by −xi and V2 by −V2. The scalar QES N -Body Spin Hamiltonians with Short-Range Interaction Potentials 3 reductions Hsc ε of these models are obtained from the above Hamiltonians by the substitution Si,i+1 → 1. A few remarks on the configuration spaces of these models are now in order. In all three models the potential diverges as (xi−xi+1)−2 on the hyperplanes xi = xi+1, so that the particles i and i+1 cannot overtake one another. Since we are interested in models with nearest and next- to-nearest neighbors interactions, we shall henceforth assume that x1 < · · · < xN . For the second potential (2b) we shall take in addition x1 > 0, due to the double pole at xi = 0. A first observation concerning the eigenstates of the models (1) is that if ψ is an eigenfunction of Hsc ε with energy E, then the factorized state Ψ = ψ|s〉 is an eigenfunction of Hε with the same energy for any spin state |s〉 symmetric under permutation of particles. A second observation is that H2 commutes with the total momentum P = −i ∑ i ∂xi . Hence the movement of the center of mass decouples, and only the eigenfunctions of H2 with zero total momentum need to be considered. The next theorem summarizes the main results presented in [13, 14], which yield a fully explicit description of several families of algebraic eigenfunctions. It is not difficult to realize, however, that these eigenfunctions do not exhaust the whole spectrum of the models. We shall denote by Λ the projection operator on states totally symmetric under the simultaneous permutation of both the spatial and spin coordinates. We define the spin vectors |si〉, |sij〉 as Λ(x1|s〉) = ∑ i xi|si〉, Λ(x1x2|s〉) = ∑ i<j xixj |sij〉. We shall consider the subspace Σ′ ⊂ Σ of spin vectors |s〉 such that ∑ i |si〉 is symmetric. A thorough characterization of this space is given in [14]. Theorem 1. Let l, m be nonnegative integers and let |s〉 be an arbitrary spin vector, and denote by x the center of mass coordinate 1 N ∑ i xi. Then the following statements hold: 1. Let α = N(a + 1 2) − 3 2 , β ≡ β(m) = 1 − m − N(a + 1 2), t = 2r2 Nx2 − 1, and µlm = e− ω 2 r2 xmL−β l (ωr2) ∏ i |xi − xi+1|a. The Hamiltonian H0 possesses the following families of spin eigenfunctions with eigenvalue Elm = E0 + 2ω(2l +m), where E0 = Nω(2a + 1) is the ground state energy: Ψ(0) lm = µlmP (α,β) [m 2 ] (t)Λ|s〉, m ≥ 0, Ψ(1) lm = µl,m−1 [ P (α+1,β) [m−1 2 ] (t)Λ(x1|s〉)− xP (α+1,β) [m−1 2 ] (t)Λ|s〉 ] , m ≥ 1, Ψ(2) lm = µl,m−2 [ P (α+2,β) [m 2 ]−1 (t) ( Λ(x2 1|s〉)− 2xΛ(x1|s〉) ) + x2 ( P (α+2,β) [m 2 ]−1 (t)− 2(α+ 1) 2[m−1 2 ] + 1 P (α+1,β) [m 2 ]−1 (t) ) Λ|s〉 ] , m ≥ 2, Ψ(3) lm = µl,m−3 [ 2 3N P (α+3,β) [m−3 2 ] (t) ∑ i x3 i + x3ϕm(t) ] Λ|s〉, m ≥ 3, Ψ(4) lm = µl,m−4 [ 3 2([m−3 2 ] + 1 2) x2P (α+3,β) [m 2 ]−2 (t)Λ(x2 1|s〉) + (3 2 x3φm(t)− 1 N P (α+4,β) [m 2 ]−2 (t) ∑ i x3 i ) Λ(x1|s〉) + ( 1 N xP (α+4,β) [m 2 ]−2 (t) ∑ i x3 i + 3 2 x4χm(t) ) Λ|s〉 ] , m ≥ 4. 4 A. Enciso, F. Finkel, A. González-López and M.A. Rodŕıguez When |s〉 ∈ Σ′ there is an additional family of eigenfunctions given by Ψ̃(2) lm = µl,m−2 [ P (α+2,β) [m 2 ]−1 (t) ( Λ(x1x2|s〉)− 2xΛ(x1|s〉) ) + x2 ( P (α+2,β) [m 2 ]−1 (t) + 2(α+ 1)( 2[m−1 2 ] + 1 ) (N − 1) P (α+1,β) [m 2 ]−1 (t) ) Λ|s〉 ] , m ≥ 2. The functions ϕm, φm and χm are polynomials given explicitly by ϕm = m+ 2α+ 2 m− 1 P (α+2,β−2) m 2 − P (α+3,β−1) m 2 −1 − 4α+ 7 m− 1 P (α+2,β−1) m 2 −1 + 1 3 P (α+3,β) m 2 −2 , φm = P (α+4,β−1) m 2 −1 − 2P (α+3,β−1) m 2 −1 − m+ 2α+ 3 (m− 1)(m− 3) P (α+2,β−1) m 2 −1 − 1 3 P (α+4,β) m 2 −2 + m+ 2α− 1 m− 3 P (α+3,β) m 2 −2 , χm = 3m+ 2α (m− 1)(m− 3) P (α+2,β−1) m 2 −1 + 2m− 7 m− 3 P (α+3,β−1) m 2 −1 − P (α+4,β−1) m 2 −1 − m+ 2α+ 2 (m− 1)(m− 3) P (α+2,β) m 2 −2 − m+ 2α m− 3 P (α+3,β) m 2 −2 + 1 3 P (α+4,β) m 2 −2 , for even m, and ϕm = 2P (α+2,β−1) m−1 2 − P (α+3,β−1) m−1 2 + 1 3 P (α+3,β) m−3 2 + m+ 2α+ 2 m(m− 2) P (α+1,β) m−3 2 − m+ 2α+ 2 m− 2 P (α+2,β) m−3 2 , φm = P (α+4,β−1) m−3 2 − 2m− 5 m− 2 P (α+3,β) m−3 2 − 1 3 P (α+4,β) m−5 2 + m+ 2α− 1 m− 2 P (α+3,β) m−5 2 , χm = 2m− 3 m(m− 2) P (α+2,β−1) m−3 2 + 2(m− 3) m− 2 P (α+3,β−1) m−3 2 − P (α+4,β−1) m−3 2 − m+ 2α+ 1 m(m− 2) P (α+2,β) m−5 2 − m+ 2α m− 2 P (α+3,β) m−5 2 + 1 3 P (α+4,β) m−3 2 , for odd m. 2. The Hamiltonian H1 possesses the following families of spin eigenfunctions with eigenvalue Ek = E0 + 4kω, where E0 = Nω(4a+ 2b+ 1) is the ground state energy: Ψ(0) k = µLα−1 k (ωr2)Λ|s〉, k ≥ 0, Ψ(1) k = µLα+1 k−1 (ωr2) [ NΛ(x2 1|s〉)− r2Λ|s〉 ] , k ≥ 1, Ψ(2) k = µLα+3 k−2 (ωr2) [ N(α+ 1) ∑ i x4 i − βr4 ] Λ|s〉, k ≥ 2, with α = N(2a+ b+ 1 2), β = N(4a+ b+ 3 2) and µ = e− ω 2 r2 ∏ i |x2 i − x2 i+1| a xb i . 3. The Hamiltonian H2 possesses the following spin eigenfunctions with zero momentum Ψ0 = µΦ(0), Ψ1,2 = µ ∑ i { cos sin }( 2(xi − x) ) |si〉, Ψ3 = µ [ 2a 2a+ 1 Φ(0) + ∑ i6=j cos ( 2(xi − xj) ) |sj〉 ] , Ψ4 = µ ∑ i6=j sin ( 2(xi − xj) ) |sj〉, where µ = ∏ i sin a |xi − xi+1|. Their energies are respectively given by E0, E1,2 = E0 + 4 ( 2a+ 1− 1 N ) , E3,4 = E0 + 8(2a+ 1), where E0 = 2Na2 is the ground state energy. QES N -Body Spin Hamiltonians with Short-Range Interaction Potentials 5 3 Calogero models and Dunkl operators Let us consider a self-adjoint operator H acting on a given Hilbert space H. The idea underlying the construction of algebraic eigenfunctions of H (see, e.g., [40]) is that the explicit knowledge of a finite-dimensional subspace H1 ⊂ H which is invariant under the operator H allow ones to compute dimH1 eigenfunctions and eigenvalues of H by algebraic methods, i.e., diagonalizing the matrix H|H1 . In this case the operator is said to be quasi-exactly solvable (QES). In fact, the Hamiltonians of the models (2a) and (2b) possess an infinite flag H1 ⊂ H2 ⊂ · · · of known finite-dimensional invariant subspaces, which yield an arbitrary large number of eigenvalues and eigenfunctions. Although such models are sometimes termed “exactly solvable” [39], we will not use this terminology since the algebraic eigenfunctions of the models (2a) and (2b) are not an orthonormal basis of the Hilbert space. A particularly convenient method of carrying out this program [38, 35] is through a Lie algebra g of first-order differential operators Ja = ∑ i ξ ai(z)∂zi + ηa(z) (a = 1, . . . , r) acting on a finite-dimensional module M⊂ C∞(M), with M a domain in RN . Let us assume that there exists a second-order differential operator H̃ = r∑ a,b=1 cabJ aJb + r∑ a=1 caJ a + c0 that is equivalent, up to gauge transformation H̃ 7→ µH̃µ−1 (µ ∈ C∞(M,R+)) and global change of variables z ∈ M 7→ x ∈ RN , to a Schrödinger operator H = −∆g + V , ∆g standing for the Laplace–Beltrami operator in (RN , g). If µM|z 7→x ⊂ L2(RN , √ gdx), then one can obtain dimM eigenfunctions of H by algebraic methods. The Schrödinger operators amenable to this treatment are termed Lie-algebraic. It should be observed that the operators Ja are not symmetries of H: g is an algebra of hidden symmetries of H. Dunkl operators [20] were originally introduced to study spherical harmonics associated to measures invariant under a Coxeter group [10]. Actually, let v,z ∈ RN and define the reflection σvz = z − 2|v|−2(z · v)v. Let W be a Coxeter group, which can be assumed to be the Weyl group of a (possibly nonreduced) root system R, and let va = ∑ i v i aei (a = 1, . . . , r) be a basis of positive roots. Define an action of W on C[z] as Kaf = f ◦ σva . Dunkl operators were originally defined as Ji = ∂zi + r∑ a=1 g2 av i a z · va (1−Ka), (3) where the real parameters g2 a are chosen so that they are constant on each orbit of W. The operators Ji can be understood as deformations of the partial derivatives ∂zi that commute with the deformed Laplacian, i.e., ∑ i J 2 i Jj = Jj ∑ i J 2 i . It can be verified that {Ji,Ka} span a degenerate Hecke algebra [24]. Currently, the definition of Dunkl operators has been generalized to mean a set {Ji}N i=1 ⊂ End C[z] of first order differential operators which leave invariant finite- dimensional polynomial subspaces and such that {Ji,Ka} span a Hecke algebra. In their original form, Dunkl operators are directly connected with the rational Calogero model of type R [11]. In the simplest case (R = AN−1, W = SN ) and writing x instead of z, the Dunkl operators (3) read Ji = ∂xi + g2 ∑ j 6=i 1 xi − xj (1−Kij), where (Kijf)(x1, . . . , xi, . . . , xj , . . . , xN ) = f(x1, . . . , xj , . . . , xi, . . . , xN ) (4) denotes the reflection operator associated to the root ei − ej . Actually, let us consider the 6 A. Enciso, F. Finkel, A. González-López and M.A. Rodŕıguez Calogero ground state function µ(x) = e−ωr2 ∏ i<j |xi − xj |a and the auxiliary operator J0 = ∑ i xi∂xi . (5) A straightforward calculation shows that one can use the ground state function µ to gauge transform the deformed Laplacian so that µHCµ −1 = − ∑ i ∂2 xi + a ∑ i6=j 1 (xi − xj)2 (a−Kij) + ωr2, (6) with g2 = a(a− 1) and HC = ∑ i J2 i + 2ωJ0 + E0, yields the Calogero model of AN−1 type when acting on symmetric functions ψ ∈ ΛL2(RN ). Since Ji, J0 preserve the space of polynomials Pn = {f ∈ C[x] : deg f ≤ n} (7) for any n = 0, 1, . . . , the Dunkl operators provide a very convenient fashion of exploring the solvability properties of Calogero–Sutherland models. It should be remarked, however, that for an arbitrary Coxeter group W the Dunkl operators (3) do not form a Lie algebra; nevertheless, this technique captures most of the relevant features of the Lie-algebraic method. As we shall now outline, Dunkl operators can also be used to introduce internal degrees of freedom in the picture without breaking the solvability properties of the models. Given a scalar differential-difference operator D linear in Kij , let us denote by D∗ the differential operator acting on C∞ ⊗ Σ obtained from D by the replacement Kij → Sij . It is clear that the actions of D and D∗ coincide on the (bosonic) Hilbert space H = Λ(L2(RN )⊗ Σ). (8) Therefore the scalar operator (6) coincides with that of the spin Calogero model HC ≡ µH∗ Cµ −1 = − ∑ i ∂2 xi + a ∑ i6=j 1 (xi − xj)2 (a− Sij) when acting on symmetric states Ψ ∈ H. 4 Invariant subspaces The proof of the main theorem rests on the construction of appropriate invariant flags for the Hamiltonians (1) using a modification of the Dunkl operator formalism. This modification turns out to be rather nontrivial, ultimately due to the fact the cyclic group is not of Coxeter type. As the first step, we consider the second-order differential-difference operators Tε given by Tε = ∑ i zε i∂ 2 i + 2a ∑ i 1 zi − zi+1 (zε i∂i − zε i+1∂i+1)− 2a ∑ i ϑε(zi, zi+1) (zi − zi+1)2 (1−Ki,i+1), (9) where ∂i = ∂zi , zN+1 ≡ z1, and ϑ0(x, y) = 1, ϑ1(x, y) = 1 2 (x+ y), ϑ2(x, y) = xy. QES N -Body Spin Hamiltonians with Short-Range Interaction Potentials 7 Each Hamiltonian Hε is related to a linear combination Hε = cTε + c−J − + c0J 0 + E0 (10) of its corresponding operator Tε and the auxiliary first-order differential operators J− = ∑ i ∂i and J0 = ∑ i zi∂i via a change of variables, a gauge transformation, and the star mapping defined in the previous section, that is, Hε = µ ·H∗ ε ∣∣ zi=ζ(xi) · µ−1. (11) The constants c, c−, c0, E0, the gauge factor µ, and the change of variables ζ for each model are listed in Table 1. Hence the construction of the models (1) is analogous to that of the usual Calogero–Sutherland models, with the operators (9) being a cyclic analog of the sum of the squares of the Dunkl operators. Table 1. Parameters, gauge factor and change of variable in equations (10) and (11). ε = 0 ε = 1 ε = 2 c −1 −4 4 c− 0 −2(2b+ 1) 0 c0 2ω 4ω 4(1− 2a) E0 Nω(2a+ 1) Nω(4a+ 2b+ 1) 2Na2 µ(x) e− ω 2 r2 ∏ i |xi − xi+1|a e− ω 2 r2 ∏ i |x2 i − x2 i+1| a xb i ∏ i sina |xi − xi+1| ζ(x) x x2 e±2ix It is obvious that the operators (10) preserve the polynomial space (7) for any n ∈ N, so one may be tempted to believe that the usual arguments for spin CS models should yield an invariant flag for the Hamiltonians (1). Nevertheless, this is not the case. In fact, the standard construction [16, 17] is based on the fact that the actions of HC and µHCµ −1 on symmetric states coincide and these operators commute with the SN symmetrizer Λ. Unfortunately, Hε (or Hε) do not commute with Λ and cyclic symmetry does not suffice to exchange Hε and µHεµ −1, so this procedure does not grant the existence of any nontrivial invariant subspaces, not even of direct products M⊗ ΛΣ (M⊂ C∞(RN )). We shall now review the actual construction of the invariant spaces, which is considerably more involved. We shall not provide complete proofs, but merely a sketch of the main logical steps the construction rests upon. Classical results on the theory of invariants [41] are responsible for the success of studying CS models through symmetric polynomials [27, 33]. We shall extend this approach to deal with the Hamiltonians (1). Let us first introduce two bases {σk} and {τk} of the space of symmetric polynomials in z: σk = ∑ i zk i , τk = ∑ i1<···<ik zi1 · · · zik ; k = 1, . . . , N. The operators Tε consist of three summands which are of second, first and zeroth order in the derivatives. Let us denote each summand by Lε, 2aXε and −2aAε respectively. It is not difficult to realize that {Lε} span a Lie algebra isomorphic to sl(2), as in the CS case. It turns out that the first-order differential operators Xε leave invariant a flag of symmetric polynomials, as stated in the following easy lemma. 8 A. Enciso, F. Finkel, A. González-López and M.A. Rodŕıguez Lemma 1. For each n = 0, 1, . . . , the operator Xε leaves invariant the linear space X n ε , where X n 0 = C[σ1, σ2, σ3] ∩ Pn, X n 1 = C[σ1, σ2, τN ] ∩ Pn, X n 2 = C[σ1, τN−1, τN ] ∩ Pn. Remark 1. It should be noted that these flags cannot be trivially enlarged, since, e.g., 1 4 X0σ4 = 2σ2 + ∑ i zizi+1, 1 3 X1σ3 = 2σ2 + ∑ i zizi+1, X1τN−1 = τN ∑ i (zizi+1)−1, 1 2 X2σ2 = 2σ2 + ∑ i zizi+1, X2τN−2 = NτN−2 − τN ∑ i (zizi+1)−1 are not symmetric polynomials. In the next proposition we characterize subspaces of the flags described in Lemma 1 that are preserved by the whole operator Tε. If f ∈ C[σ1, σ2, σ3, τN−1, τN ], we adopt the convenient notation fk = { ∂σk f, k = 1, 2, 3, ∂τk f, k = N − 1, N. Proposition 1. For each n = 0, 1, . . . , the operator Tε leaves invariant the linear space Sn ε , where Sn 0 = {f ∈ X n 0 | f33 = 0}, Sn 1 = {f ∈ X n 1 | f22 = fNN = 0}, Sn 2 = {f ∈ X n 2 | f11 = fN−1,N−1 = 0}. Proposition 1 implies that each operator Tε preserves product symmetric subspaces Sn ε ⊗ΛΣ spanned by factorized states. The main result on invariant subspaces shows that in fact the latter operator leaves invariant a richer flag of nontrivial finite-dimensional subspaces of Λ(Pn ⊗ Σ). Theorem 2. Let T n 0 = 〈 f(σ1, σ2, σ3)Λ|s〉, g(σ1, σ2, σ3)Λ(z1|s〉), h(σ1, σ2)Λ(z2 1 |s〉), h̃(σ1, σ2)Λ(z1z2|s′〉) | f33 = g33 = 0 〉 , T n 1 = 〈 f(σ1, σ2, τN )Λ|s〉, g(σ1, τN )Λ(z1|s〉) | f22 = fNN = gNN = 0 〉 , T n 2 = 〈 f(σ1, τN−1, τN )Λ|s〉, g(τN−1, τN )Λ(z1|s〉), τNq(σ1, τN )Λ(z−1 1 |s〉) | f11 = fN−1,N−1 = gN−1,N−1 = q11 = 0 〉 , where |s〉 ∈ Σ, |s′〉 ∈ Σ′, deg f ≤ n, deg g ≤ n−1, deg h ≤ n−2, deg h̃ ≤ n−2, deg q ≤ n−N+1, and deg is the total degree in z. Then T n ε is invariant under Tε for all n = 0, 1, . . . . From this theorem one easily obtains the following corollary, which is crucial for the compu- tation of the algebraic eigenfunctions of the models (1). Corollary 1. For each ε = 0, 1, 2, the gauge Hamiltonian Hε leaves invariant the space Hn ε defined by Hn 0 = T n 0 , Hn 1 = T n 1 ∣∣ fN=gN=0 , Hn 2 = T n 2 . (12) QES N -Body Spin Hamiltonians with Short-Range Interaction Potentials 9 5 Spectrum and eigenfunctions As happens with the usual CS models, the algebraic eigenvalues of the Hamiltonians (1) can be easily obtained by choosing a basis of the invariant spaces µHn ε in which the action of Hε is triangular. In all three cases, the algebraic eigenvalue E0 is the ground state energy, since the corresponding eigenfunctions do not vanish in the configuration space Cε. Now we shall outline how the algebraic eigenfunctions in Theorem 1 were calculated. The easiest case is ε = 2, since states related by a multiplicative factor τk N only differ by the movement of the center of mass, which is conserved. Hence the invariant subspaces Hn 2 solely allow one to compute five eigenfunctions of zero total momentum, which are the ones listed in the main theorem. In the case ε = 1 the eigenvalue equation reads (H1 − E0 − 4ωn)Φ = 0, where Φ is a symmetric vector-valued polynomial in Hn 1 of degree n. Setting t = ωσ1 and Φ = [p(t) + σ2q(t)]Λ|s〉+ g(t)Λ(z1|s〉), this equation can be written as Lα+1 k−1g = Lα+3 k−2q = 0, Lα−1 k p = − α Nω g − 2β Nω2 tq, (13) with the Laguerre operator Lλ ν defined as (Lλ νf)(t) = tf ′′(t) + (λ+ 1− t)f ′(t) + νf(t). With some effort one can obtain all the solutions of these equations, as shown in the following proposition. Proposition 2. The polynomial solutions of the system of ODE’s (13) are spanned by Φ(0) n = Lα−1 n (t)Λ|s〉, n ≥ 0, Φ(1) n = Lα+1 n−1(t) [ NωΛ(z1|s〉)− tΛ(z1|s〉) ] , n ≥ 1, Φ(2) n = Lα+3 n−2(t) [ Nω2(α+ 1)σ2 − βt2 ] Λ|s〉, n ≥ 2. These solutions are characterized by the conditions q = g = 0, q = 0 and g = 0 respectively and correspond to the algebraic eigenfunctions of H1 presented in Theorem 1. The case ε = 0 is similar, but the computations become more involved due to the rich structure of the invariant flag. Writing Φ = (p+ σ3q)Λ|s〉+ (u+ σ3v)Λ(z1|s〉) + hΛ(z2 1 |s〉) + h̃Λ(x1x2|s〉), where deg Φ = k and h̃ = 0 if |s〉 6∈ Σ′, the equation (H0 − E0 − 2ωk)Φ = 0 reduces to the systems of PDE’s[ L0 − 2ω(k − 2) ] h̃− 8h̃2 = 0, (14a)[ L0 − 2ω(k − 2) ] h− 8h2 = 6v, (14b)[ L0 − 2ω(k − 1) ] u− 4u2 = 4h1 + 4h̃1 + 6σ2v1 + 6(2a+ 1)σ1v, (14c)[ L0 − 2ω(k − 4) ] v − 16v2 = 0, (14d)( L0 − 2ωk ) p = 2u1 + 2(2a+ 1)h− 4a N − 1 h̃+ 6σ2q1 + 6(2a+ 1)σ1q, (14e)[ L0 − 2ω(k − 3) ] q − 12q2 = 2v1, (14f) 10 A. Enciso, F. Finkel, A. González-López and M.A. Rodŕıguez with L0 = − ( N∂2 σ1 + 4σ1∂σ1∂σ2 + 4σ2∂ 2 σ2 + 2(2a+ 1)N∂σ2 ) + 2ω(σ1∂σ1 + 2σ2∂σ2). One can check by inspection that the system (14) possesses six families of polynomial solutions, as collected in Table 2. These correspond to the six families of algebraic eigenfunctions listed in Theorem 1. Table 2. The six types of polynomial solutions of the system (14) and their corresponding eigenfunctions. Conditions Corresponding eigenfunction q = u = v = h = h̃ = 0, p 6= 0 Ψ(0) lm u = v = h = h̃ = 0, q 6= 0 Ψ(3) lm q = v = h = h̃ = 0, u 6= 0 Ψ(1) lm q = v = h̃ = 0, h 6= 0 Ψ(2) lm q = v = h = 0, h̃ 6= 0 Ψ̃(2) lm h̃ = 0, v 6= 0 Ψ(4) lm Acknowledgements This work was partially supported by the DGI under grant no. FIS2005-00752. A.E. acknow- ledges the financial support of the Spanish Ministry of Education through an FPU scholarship. [1] Auberson G., Jain S.R., Khare A., A class of N -body problems with nearest- and next-to-nearest-neighbour interactions, J. Phys. A: Math. Gen., 2001, V.34, 695–724, cond-mat/0004012. [2] Azuma H., Iso S., Explicit relation of the quantum Hall effect and the Calogero–Sutherland model, Phys. Lett. B, 1994, V.331, 107–113, hep-th/9312001. [3] Baker T.H., Forrester P.J., The Calogero–Sutherland model and generalized classical polynomials, Comm. Math. Phys., 1997, V.188, 175–216, solv-int/9608004. [4] Bogomolny E.B., Gerland U., Schmit C., Models of intermediate statistics, Phys. Rev. E, 1999, V.59, R1315– R1318. [5] Brink L., Turbiner A., Wyllard N., Hidden algebras of the (super) Calogero and Sutherland models, J. Math. Phys., 1998, V.39 1285–1315, hep-th/9705219. [6] Calogero F., Solution of the one-dimensional N -body problem with quadratic and/or inversely quadratic pair potentials, J. Math. Phys., 1971, V.12, 419–436. [7] Deguchi T., Ghosh P.K., Spin chains from super-models, J. Phys. Soc. Japan, 2001, V.70, 3225–3237, hep-th/0012058. [8] Desrosiers P., Lapointe L., Mathieu P., Supersymmetric Calogero–Moser–Sutherland models and Jack su- perpolynomials, Nucl. Phys. B, 2001, V.606, 547–582, hep-th/0103178. [9] D’Hoker E., Phong D.H., Calogero–Moser systems in SU(N) Seiberg–Witten theory, Nucl. Phys. B, 1998, V.513, 405–444, hep-th/9709053. [10] Dunkl C.F., Differential-difference operators associated to reflection groups, Trans. Amer. Math. Soc., 1989, V.311, 167–183. [11] Dunkl C.F., Orthogonal polynomials of types A and B and related Calogero models, Comm. Math. Phys., 1998, V.197, 451–487, q-alg/9710015. [12] Enciso A., Finkel F., González-López A., Rodŕıguez M.A., Haldane–Shastry spin chains of BCN type, Nucl. Phys. B, 2005, V.707, 553–576, hep-th/0406054. [13] Enciso A., Finkel F., González-López A., Rodŕıguez M.A., Solvable scalar and spin models with near- neighbors interactions. Phys. Lett. B, 2005, V.605, 214–222, hep-th/0407274. [14] Enciso A., Finkel F., González-López A., Rodŕıguez M.A., Exchange operator formalism for N -body spin models with near-neighbors interactions, nlin.SI/0604073. http://arxiv.org/abs/cond-mat/0004012 http://arxiv.org/abs/hep-th/9312001 http://arxiv.org/abs/solv-int/9608004 http://arxiv.org/abs/hep-th/9705219 http://arxiv.org/abs/hep-th/0012058 http://arxiv.org/abs/hep-th/0103178 http://arxiv.org/abs/hep-th/9709053 http://arxiv.org/abs/q-alg/9710015 http://arxiv.org/abs/hep-th/0406054 http://arxiv.org/abs/hep-th/0407274 http://arxiv.org/abs/nlin.SI/0604073 QES N -Body Spin Hamiltonians with Short-Range Interaction Potentials 11 [15] Ezung M., Gurappa N., Khare A., Panigrahi P.K., Quantum many-body systems with nearest and next-to- nearest neighbor long-range interactions, Phys. Rev. B, 2005, V.71, 125121(8), cond-mat/0007005. [16] Finkel F., Gómez-Ullate D., González-López A., Rodŕıguez M.A., Zhdanov R., AN -type Dunkl operators and new spin Calogero–Sutherland models, Comm. Math. Phys., 2001, V.221, 477–497, hep-th/0102039. [17] Finkel F., Gómez-Ullate D., González-López A., Rodŕıguez M.A., Zhdanov R., New spin Calogero–Suther- land models related to BN -type Dunkl operators, Nucl. Phys. B, 2001, V.613, 472–496, hep-th/0103190. [18] Gorsky A., Nekrasov N., Hamiltonian systems of Calogero type, and two dimensional Yang–Mills theory, Nucl. Phys. B, 1994, V.414, 213–238, hep-th/9304047. [19] Haldane F.D.M., O(3) nonlinear σ model and the topological distinction between integer- and half-integer- spin antiferromagnets in two dimensions, Phys. Rev. Lett., 1988, V.60, 635–638. [20] Heckman G.J., Dunkl operators, Astérisque, 1997, V.245, 223–246. [21] Hitchin N., Stable bundles and integrable systems, Duke Math. J., 1987, V.54, 91–114. [22] Jain S.R., Khare A., An exactly solvable many-body problem in one dimension and the short-range Dyson model, Phys. Lett. A, 1999, V.262, 35–39, cond-mat/9904121. [23] Kasman, A., Bispectral KP solutions and linearization of Calogero–Moser particle systems, Comm. Math. Phys., 1995, V.172, 427–448, hep-th/9412124. [24] Kirillov A.A., Jr., Lectures on affine Hecke algebras and Macdonald’s conjectures, Bull. Amer. Math. Soc. (N.S.), 1997, V.34, 251–292, math.QA/9501219. [25] Moser J., Three integrable Hamiltonian systems connected to isospectral deformations, Adv. Math., 1975, V.16, 197–220. [26] Olshanetsky M.A., Perelomov A.M., Quantum integrable systems related to Lie algebras, Phys. Rep., 1983, V.94, 313–403. [27] Perelomov A.M., Remark on the completeness of the coherent state system, Theor. Math. Phys., 1971, V.6, 213–224. [28] Polychronakos A.P., Nonrelativistic bosonization and fractional statistics, Nucl. Phys. B, 1989, V.324, 597– 622. [29] Polychronakos A.P., Exchange operator formalism for integrable systems of particles, Phys. Rev. Lett., 1992, V.69, 703–706, hep-th/9202057. [30] Polychronakos A.P., Lattice integrable systems of Haldane–Shastry type, Phys. Rev. Lett., 1993, V.70, 2329–2332, hep-th/9210109. [31] Polychronakos A.P., Exact spectrum of an SU(N) spin chain with inverse-square exchange, Nucl. Phys. B, 1994, V.419, 553–566. [32] Polychronakos A.P., Waves and solitons in the continuum limit of the Calogero–Sutherland model, 1995, Phys. Rev. Lett., 1995, V.74, 5153–5156, hep-th/9411054. [33] Rühl W., Turbiner A.V., Exact solvability of the Calogero and Sutherland models, Mod. Phys. Lett. A, 1995, V.10, 2213–2221, hep-th/9506105. [34] Shastry B.S., Exact solution of an S = 1/2 Heisemberg antiferromagnetic chain with long-ranged interac- tions, Phys. Rev. Lett., 1988, V.60, 639–642. [35] Shifman M.A., Turbiner A.V., Quantal problems with partial algebraization of the spectrum, Comm. Math. Phys., 1989, V.126, 347–365. [36] Shukla P., Non-hermitian random matrices and the Calogero–Sutherland model, Phys. Rev. Lett., 2001, V.87, 194102(5). [37] Sutherland B., Exact results for a quantum many-body problem in one dimension, I, II, Phys. Rev. A, 1971, V.4, 2019–2021, 1972, V.5, 1372–1376. [38] Turbiner A.V., Quasi-exactly solvable problems and sl(2) algebra, Comm. Math. Phys., 1988, V.118, 467– 474. [39] Turbiner A.V., Lie algebras and polynomials in one variable, J. Phys. A: Math. Gen., 1992, V.25, L1087– L1093. [40] Turbiner A.V., Lie algebras and linear operators with invariant subspaces, in Lie Algebras, Cohomologies and New Findings in Quantum Mechanics, Editors N. Kamran and P.J. Olver, Contemporary Mathematics, Vol. 160, Providence, AMS, 1994, 263–310, funct-an/9301001. [41] Weyl H., The classical groups, Princeton, Princeton University Press, 1997. [42] Yamamoto T., Multicomponent Calogero model of BN -type confined in a harmonic potential, Phys. Lett. A, 1995, V.208, 293–302, cond-mat/9508012. http://arxiv.org/abs/cond-mat/0007005 http://arxiv.org/abs/hep-th/0102039 http://arxiv.org/abs/hep-th/0103190 http://arxiv.org/abs/hep-th/9304047 http://arxiv.org/abs/cond-mat/9904121 http://arxiv.org/abs/hep-th/9412124 http://arxiv.org/abs/math.QA/9501219 http://arxiv.org/abs/hep-th/9202057 http://arxiv.org/abs/hep-th/9210109 http://arxiv.org/abs/hep-th/9411054 http://arxiv.org/abs/hep-th/9506105 http://arxiv.org/abs/funct-an/93010012 http://arxiv.org/abs/cond-mat/9508012 1 Introduction 2 Main result 3 Calogero models and Dunkl operators 4 Invariant subspaces 5 Spectrum and eigenfunctions