arxiv: v2 [math-ph] 10 Dec 2015

Transcription

1 A superintegrable discrete harmonic oscillator based on bivariate Charlier polynomials Vincent X. Genest, 1, Hiroshi Miki,, Luc Vinet, 3,4, and Guofu Yu 4, 1 Department of Mathematics, Massachusetts Institute of Technology, arxiv: v math-ph] 10 Dec Massachusetts Ave, Cambridge, MA 0135 USA Department of Electronics, Faculty of Science and Engineering, Doshisha University, Kyotanabe City, Kyoto Japan 3 Centre de Recherches Mathématiques, Université de Montréal, P. O. Box 618, Centre-ville Station, Montréal, QC, H3C 3J7, Canada 4 Department of Mathematics, Shanghai Jiao Tong University, Shanghai 0040 China A simple discrete model of the two dimensional isotropic harmonic oscillator is presented. It is superintegrable with su() as its symmetry algebra. It is constructed with the help of the algebraic properties of the bivariate Charlier polynomials. This adds to the other discrete superintegrable models of the oscillator based on Krawtchouk and Meixner orthogonal polynomials in several variables. I. INTRODUCTION For various applications, particularly in optical image processing and in signal analysis, it is useful to have discrete models of the harmonic oscillator that preserve most of the properties of the system in the continuum. Several such discrete models of the one-dimensional quantum oscillator have been developed and the reader might consult ref. 1] for a review. Of particular relevance to the present paper is the model based on the one-variable Charlier polynomials that is constructed in ], see also 3]. It is desirable to work out models in Electronic address: vxgenest@mit.edu Electronic address: hmiki@mail.doshisha.ac.jp Electronic address: vinet@crm.umontreal.ca Electronic address: gfyu@sjtu.edu.cn

2 two dimensions as discrete planar oscillators can be employed for instance to describe wave guides and pixellated screens. In general, combining two discrete one-dimensional oscillator models does not lead to a system that has the symmetries of the two-dimensional problem in the continuum 4, 5]. (Interestingly, we shall here have a counterexample to that statement.) The isotropic harmonic oscillator is the paradigm example of a maximally superintegrable system. In d dimensions, such quantum systems with Hamiltonian H, possess d 1 algebraically independent constants of motion including H. These quantities that commute with H necessarily form a non-abelian algebra and in the case of the isotropic harmonic oscillator they generate the Lie algebra su(d). Superintegrable systems are the hallmark of exactly solvable models, they represent a meaningful laboratory for the study of symmetries and their mathematical description and they prove highly useful in applications. It is thus of considerable interest to enlarge their catalogue. For a review of (mostly continuous) superintegrable systems, see 6]. In the recent past the poor set of known discrete superintegrable models has been significantly increased. Indeed, two discrete models of the -dimensional isotropic harmonic oscillator with su() symmetry have been obtained 7, 8]. While the two-dimensional situation is described in details in the papers quoted, these models have natural extensions to arbitrary dimension d. Their construction crucially relies on the properties of orthogonal polynomials in d variables that have recently received algebraic interpretations. In the first case 7], one called upon the multivariate Krawtchouk polynomials now known 9] to arise as matrix elements of representations of SO(d + 1) on oscillator states. In the second case 8], use was made of the Meixner polynomials in d variables that are similarly interpreted 10] with the help of the pseudo-orthogonal group SO(d,1). In that vein, some of us have also given 11] an algebraic description of the Charlier polynomials in d variables based on the Euclidean group in d dimensions. We here exploit the properties of these polynomials obtained in 11] to construct another discrete model of the isotropic harmonic oscillator thereby complementing the two models we have mentioned. The outline of the paper is as follows. In section, we will review the features of the bivariate Charlier polynomials that stem from their interpretation in terms of E() representations. In section 3, the discrete model of the isotropic oscillator will be defined and shown to be superintegrable with an su() invariance algebra. The wavefunctions will be recorded. In section 4, it will be confirmed that we have a model of the oscillator by showing

3 3 that this system is found in the continuum limit. We shall conclude with some remarks. An Appendix collecting the relevant properties of the orthonormalized Charlier polynomials in one variable is also included. II. THE TWO-VARIABLE CHARLIER POLYNOMIALS The bivariate Charlier polynomials have been seen in 11] to arise in matrix elements of representations of the Euclidean group in two dimensions. Let θ,α and β be the rotation angle and the translations that define a generic Euclidean transformation of the plane with coordinates x and y x cosθx+sinθy + α, y sinθx+cosθy + β. (1a) (1b) These will occur as parameters of the bivariate Charlier polynomials C n1,n (x 1,x ) that can be defined by the generating function e z 1(αcosθ βsinθ) e z (αsinθ+βcosθ) = n 1,n =0 1+ z 1 α cosθ+ z α sinθ ] x1 1 z 1 β sinθ+ z ] x β cosθ C n1,n (x 1,x ) zn 1 1 z n () n1!n!. It can be seen from () that C n1,n (x 1,x ) are polynomials of total degree n 1 + n in the variables x 1 and x with values in N (x 1,x = 0,1,,...). They satisfy the orthogonality relation x 1,x =0 w(x 1,x )C n1,n (x 1,x )C n 1,n (x 1,x ) = δ n1,n 1 δ n,n, (3) where w(x 1,x ) is the product of two independent Poisson distribution with parameters α and β : w(x 1,x ) = e (α +β ) αx 1 β x. (4) x 1!x! The polynomials C n1,n (x 1,x ) can be seen from their generating function to have the following explicit expression as an Aomoto-Gel fand hypergeometric series: C n1,n (x 1,x ) = ( 1)n 1+n n1!n! (αcosθ βsinθ)n 1 (αsinθ+βcosθ) n ρ,σ,µ,ν=0 ( n 1 ) ρ+µ ( n ) σ+ν ( x 1 ) ρ+σ ( x ) µ+ν u ρ ρ!σ!µ!ν! 11u σ 1u µ 1u ν, (5)

4 4 where the parameters u ij are given by and where (a) n stands for cosθ u 11 = α cosθ αβsinθ, u sinθ 1 = α sinθ+αβcosθ, (6) sinθ u 1 = β sinθ αβcosθ, u cosθ = β cosθ+αβsinθ, (a) n = a(a+1) (a+n 1), (a) 0 = 1. (7) The algebraic description of the multivariate Charlier polynomials allows for a very natural characterization of these functions. It leads in particular to an immediate identification of ladder operators that shall prove essential in the construction the superintegrable model. Let T ± x i f(x i ) = f(x i ±1) for i = 1, denote the discrete shift operators in the variable x 1 and x. Introduce the following two operators A (1) + = x 1 α cosθt x 1 x β sinθt x (αcosθ βsinθ), A () + = x 1 α sinθt x 1 + x β cosθt x (αsinθ +βcosθ). (8a) (8b) They act as raising operators on the polynomials C n1,n (x 1,x ): A (1) + C n1,n (x 1,x ) = n 1 +1C n1 +1,n (x 1,x ), (9a) A () + C n1,n (x 1,x ) = n +1C n1,n +1(x 1,x ). (9b) Introduce similarly the operators A (1) and A () defined by A (1) = αcosθt + x 1 βsinθt + x (αcosθ βsinθ), A () = αsinθt + x 1 +βcosθt + x (αsinθ+βcosθ). (10a) (10b) They are the lowering operators: A (1) C n1,n (x 1,x ) = n 1 C n1 1,n (x 1,x ), A () C n1,n (x 1,x ) = n C n1,n 1(x 1,x ). (11a) (11b) The ladder operators(8) and (10) can be combined to produce the two commuting difference operators Y 1 and Y that are diagonalized by the bivariate Charlier polynomials. These operators are given by Y i = A (i) + A (i), i = 1,. (1)

5 5 Explicitly they read β Y 1 = x 1 α cosθsinθt x 1 T x + α x β cosθsinθt+ x 1 Tx x1 ] ] ωcosθ α T x 1 +αt x + x 1 +ωsinθ β T x +βt x + +(x 1 cos θ+x sin θ+ω ), β Y = x 1 α cosθsinθt x 1 T x + α +x β cosθsinθt+ x 1 Tx x1 ] ] ζsinθ α T x 1 +αt x + x 1 ζcosθ β T x +βt x + +(x 1 sin θ+x cos θ+ζ ), (13a) (13b) where ω = αcosθ βsinθ, ζ = αsinθ+βcosθ. (14) The eigenvalue equations are Y i C n1,n (x 1,x ) = n i C n1,n (x 1,x ) i = 1,. (15) Note that they fully determine the polynomials C n1,n (x 1,x ). Let us also write down the form of A (1) + A () and A () + A (1) as difference operators A (1) + A () = x 1 β α cos θtx 1 T x + α x β sin θt x + 1 Tx ] ζ x 1 α cosθt x 1 +αωsinθt x ((x 1 x )sinθcosθ+ωζ), A () + A (1) = x 1 β α sin θtx 1 T x + α +x β cos θt x + 1 Tx ω ] x 1 α sinθt x 1 +αζcosθt x + 1 +((x 1 x )sinθcosθ+ωζ), x ζ β sinθt x βωcosθt + x x ω β cosθt x βζsinθt + x ] ] (16a) (16b) which act on C n1,n (x 1,x ) as follows: A (1) + A () C n1,n (x 1,x ) = (n 1 +1)n C n1 +1,n 1(x 1,x ), A () + A (1) C n1,n (x 1,x ) = n 1 (n +1)C n1 1,n +1(x 1,x ). (17a) (17b)

6 6 III. A DISCRETE AND SUPERINTEGRABLE HAMILTONIAN We present in this section a model which is superintegrable with a Hamiltonian that is a difference operator in two variables. We take for the Hamiltonian the sum of the operators Y 1 and Y introduced in the previous section. We hence posit H = Y 1 +Y. (18) A simple calculation shows that H is the following operator: H = x 1 Tx 1 α T x + 1 x Tx β T x + +x 1 +x +α +β. (19) Let us remark that this is a rather simple expression, that the operators T x + 1 T x and Tx 1 T x + have cancelled in the combination Y 1 +Y, that the parameter θ does not appear and that H is simply the sum of the operators (A7b) associated to the difference equations of univariate Charlier polynomials with parameter α and β. As for the Krawtchouk 7] and the Meixner 8] models, it is straightforward to show that H is superintegrable by exhibiting constants of motion built from the ladder operators. Let J X,J Y and J Z be defined as follows J X = 1 (A(1) + A () +A () + A (1) ), J Y = 1 i (A(1) + A () A () + A (1) ), J Z = 1 (A(1) + A (1) A () + A () ). (0a) (0b) (0c) It is immediate to verify using (8) and (10) that these operators commute with the Hamiltonian (18): H,J X ] = H,J Y ] = H,J Z ] = 0. (1) One also has J X,J Y ] = ij Z, J Y,J Z ] = ij X, J Z,J X ] = ij Y. () This confirms that H is superintegrable with su() as its symmetry algebra. In the realization (0), the su() Casimir operator is related to the Hamiltonian (18) through JX +J Y +J Z = H ( ) H +1. (3) This of course is totally in keeping with the Schwinger realization of su() in terms of the creation and annihilation operators of independent oscillators. We have thus constructed

7 7 in this case, a discrete superintegrable system by the addition of two discrete uni-dimensional models of the oscillator (see ]). We shall return to this point in our concluding remarks. By construction, the wavefunctions of the Hamiltonian (18) are expressed in terms of the two-variable Charlier polynomials C n1,n (x 1,x ). These wavefunctions Φ N,n (x 1,x ) are labelled by the two non-negative integers N and n and read Φ N,n (x 1,x ) = C n,n n (x 1,x ), n = 0,...,N, N = 0,1,,... (4) They are eigenfunctions of H and J Z : HΦ N,n (x 1,x ) = NΦ N,n (x 1,x ), ( J Z Φ N,n (x 1,x ) = n N ) Φ N,n (x 1,x ). (5a) (5b) The energy eigenvalues are given by the non-negative integers N = 0,1,,... and have an (N +1)-fold degeneracy. For a fixed energy, that is, for a given N, the states Φ N,n (x 1,x ) support a (N + 1)-dimensional irreducible representation of su(). With J ± = J X ±ij Y, (6) it is seen from (17) that J + Φ N,n (x 1,x ) = (n+1)(n n)φ N,n+1 (x 1,x ), J Φ N,n (x 1,x ) = n(n n+1)φ N,n 1 (x 1,x ). (7a) (7b) It thus follows that the set of all two-variable Charlier polynomials C n1,n (x 1,x ) such that K = n 1 +n with K a fixed integer, form a (K +1)-dimensional irreducible representation space for su(). In view of (3), the wavefunctions Φ N,n (x 1,x ) are not normalized on the grid (x 1,x ) N N with respect to the standard uniform measure of quantum mechanics. Properly normalized wavefunctions Υ N,n (x 1,x ) are obtained by taking Υ N,n (x 1,x ) = w(x 1,x )C n,n n (x 1,x ), (8) where w(x 1,x ) is given by (4). One has then the orthogonality and completeness relations Υ N,n (x 1,x )Υ N,n (x 1,x ) = δ nn δ NN, (9a) x 1,x =0 N=0 n=0 N Υ N,n (x 1,x )Υ N,n (x 1,x ) = δ x 1 x δ 1 x x. (9b)

8 8 The form of the Hamiltonian or of any other operator that acts on the normalized wavefunctionsυ N,n (x 1,x ) is obtainedby performing thegaugetransformationx wxw 1 1 onthe corresponding operator which has been defined in the Charlier or Φ N,n (x 1,x ) basis. For the Hamiltonian H, this will show that the resulting operator is the sum of one-dimensional Charlier Hamiltonians as given in ]. IV. LIMIT TO THE STANDARD OSCILLATOR IN THE CONTINUUM It will now be shown by taking the continuum limit that the system provided in the last section is a discrete model of the standard harmonic oscillator. This will be done first, by observing that the discrete wavefunctions tend to the usual oscillator wavefunctions that involve a product of two Hermite polynomials. Second we shall indicate that the discrete operators also have the proper limits. We observed that the Hamiltonian H does not depend onθ. Uponsetting θ = 0in thegenerating formula (), we see withthehelp of theunivariate generating formula (A4) that C n1,n (x 1,x ) θ=0 = ( 1) n 1+n α n 1 β n Ĉ n1 (x 1 )Ĉn (x ), (30) where Ĉn(x) are the orthonormalized Charlier polynomials in one-variable defined in (A). This separation of variables suggests using the well-known limit from the univariate Charlier polynomials to the usual Hermite polynomials H n (x). Set x 1 = α x 1 +α, x = β x +β, θ = 0. (31) One has ( lim lim α β e αz 1 1+ z 1 α ) ( x1e βz 1+ z ) x = e z 1 + x 1 z 1 e z + x z. (3) β Upon comparing with the generating function of the Hermite polynomials it is immediate to see from (3) that e t +xt = lim lim C n 1,n (x, x ) θ=0 = α β n=0 H n (x) t n, (33) n! 1 n 1 +n n1!n! H n 1 ( x 1 )H n ( x ). (34)

9 9 Moreover, with the parametrization (31), it is easily found using Stirling s approximation that the Poisson distribution converges to the Gaussian distribution: Hence the wavefunctions lim e ααx1 α x 1! tend to the oscillator wavefunctions φ N,n ( x 1, x ) = = e x 1 π, lim β e ββx x! Υ N,n (x 1,x ) = w(x 1,x )C n,n n (x 1,x ) = e x π. (35) 1 x 1 + x e H n ( x 1 )H N n ( x ). (36) πn n!(n n)! This justifies calling (18) a discrete -dimensional oscillator. Now let us remark that it is not necessary to take θ = 0 in the limiting process. If one keeps θ arbitrary and performs the change of variables (31), one obtains in the limit α,β, a product of Hermite polynomials in the rotated coordinates ˆx 1 = cosθ x 1 sinθ x and ˆx = sinθ x 1 +cosθ x. Let us now examine what happens to the ladder operators under these limits. After the gauge transformations Ā(i) ± = w 1 (x 1,x )A (i) ± w 1 (x 1,x ), the raising operators become and the lowering operators read Ā (1) + = cosθ x 1 T x 1 sinθ x T x (αcosθ βsinθ), (37a) Ā () + = sinθ x 1 T x 1 +cosθ x T x (αsinθ+βcosθ) (37b) Ā (1) = cosθ x 1 +1T + x 1 sinθ x +1T + x (αcosθ βsinθ), Ā () = sinθ x 1 +1T x 1 +cosθ x +1T x (αsinθ+βcosθ). (38a) (38b) Performing again the change of variables (31) and taking the limit α,β, it is seen expectedly that the operators (37) and (38) become the rotated creation/annihilation operators where Ā (1) + cosθa + 1 sinθa+, Ā() + sinθa + 1 +cosθa+, Ā (1) cosθa 1 sinθa, Ā () sinθa 1 +cosθa, a i = x i + xi, a + i = x i xi are the usual ladder operators of the harmonic oscillator. This offers another confirmation that a discretization of the -dimensional harmonic oscillator that preserves its superintegrablity and symmetries has been provided.

10 10 V. CONCLUDING REMARKS Using properties of the bivariate Charlier polynomials, we have constructed a discrete model of the -dimensional isotropic harmonic oscillator which is superintegrable with su() as its symmetry algebra. This complements analogous models of the oscillator which are based on the bivariate Krawtchouk and Meixner polynomials. Interestingly the present model connected to the Charlier polynomials has a feature that distinguishes it from the previous ones. Indeed it is observed that while possessing the full su() symmetry of its continuum counterpart, the model provided here consists ultimately in the addition of two one-dimensional Charlier oscillators. One might wonder why this naive construction of a -dimensional system preserves the continuum symmetries in this case while this does not happen when one-dimensional Krawtchouk or Meixner ocillators are combined. The reason lies simply in the fact, noted in the Appendix, that the ladder operators of the univariate Charlier polynomials do not affect their parameter while this is not so for the Krawtchouk or Meixner polynomials. This said it should be mentioned that even though the Hamiltonian (18) does not depend on θ and can be constructed by adding the Hamiltonians of two onedimensional Charlier oscillators, its wavefunctions Φ N,n (x 1,x ) do depend on θ are not a mere product of univariate Charlier polynomials. As a result the symmetry operators (0) also depend on θ. Only in the continuum can the rotation of the symmetry of the oscillator be used to gauge away θ. Note that the -dimensional model that we discussed simply extends to arbitrary dimensions. Anisotropic discrete oscillators can also be constructed. The Hamiltonian H = k 1 Y 1 +k Y will have k 1 n 1 + k n (n 1,n = 0,1,...) as spectrum and will exhibit degeneracies whenever the frequency ratio k /k 1 is an integer. As in the continuum, these models will be superintegrable with a polynomial symmetry algebra. Explicitly one will have β H = (k k 1 )x 1 α cosθsinθt x 1 T x + α +(k k 1 )x β cosθsinθt+ x 1 Tx x1 ] (k 1 ωcosθ+k ζsinθ) α T x 1 +αt x + 1 +(k 1 ωsinθ k ζcosθ) +k 1 (x 1 cos θ+x sin θ+ω )+k (x 1 sin θ+x cos θ+ζ ) x β T x +βt + x ]

11 11 and the presence of the operators T x + 1 Tx and Tx 1 T x + should be noted in distinction to the isotropic situation. One observes that the use of the multivariate Charlier polynomials leads to a model that is more general than the one that results from the anisotropic combination of two one-dimensional Charlier oscillators. It would be of interest in the future to obtain symmetry-preserving discretizations of other exactly solvable and superintegrable models along the lines followed here. Acknowledgements HM has benefited from the hospitality of the Department of Mathematics of the Shanghai Jiao Tong University where LV gratefully holds a Chair visiting professorship. VXG holds a postdoctoral fellowship from the National Sciences and Engineering Research Council (NSERC) of Canada. HM is supported by JSPS KAKENHI Grant Number 15K The work of LV is supported in part by NSERC. The research of GFY is supported by the National Natural Science Foundation of China (grant no ). Appendix A: Properties of the univariate Charlier polynomials We shall record here for reference the properties of the Charlier polynomials in one variable that are needed. They can all be obtained from 1]. Let C n (x,a) denote the standard Charlier polynomials C n (x,a) = F 0 n, x 1 a (A1) with F 0 theusual hypergeometric series. Weshall provideformulasfortheorthonormalized polynomials Ĉ n (x) = ( a)n n! C n (x,a ). (A) Orthogonality relation Generating function x=0 e az ( 1+ z a a x Ĉ n (x)ĉn x! (x) = δ nn. ) x = n=0 Ĉ n (x) zn n!. (A3) (A4)

12 1 Ladder operators A Ĉ n (x) = nĉn 1(x), A = a(t + 1). (A5a) (A5b) A + Ĉ n (x) = n+1ĉn+1(x), A + = a+ x a T. (A6a) (A6b) Note that the ladder operators do not affect the parameter a included in the definition (A) of Ĉn(x). Difference equation YĈn(x) = nĉn+1(x), Y = A + A = a T + +(x +a ) xt. (A7a) (A7b) 1] N. Atakishiyev, G. Pogosyan and K. B. Wolf, Finite models of the oscillator, Physics of Particles and Nuclei 34, (005); ] N. Atakishiyev, E. I. Jafarov, S. M. Nagiiyev and K. B. Wolf, Meixner oscillators, Revista Mexicana de Fisica 44, (1998); 3] E. I. Jafarov and J. Van der Jeugt, The oscillator model for the Lie superintegrable sh( ) and Charlier polynomials, J. Math. Phys. 54, (013); 4] N. M. Atakishiyev and G. S. Pogosyan, L. E. Vincent, and K. B. Wolf, Finite two-dimensional oscillator: I. The Cartesian model, J. Phys. A: Math. Gen. 34, (001); 5] N. M. Atakishiyev and G. S. Pogosyan, L. E. Vincent, and K. B. Wolf, Finite two-dimensional oscillator: II. The radial model, J. Phys. A: Math. Gen. 34, (001); 6] W. Miller, S. Post and P. Winternitz, Classical and quantum superintegrability with applications, J. Phys. A: Math. Theor. 46, (013); 7] H. Miki, S. Post, L. Vinet and A. Zhedanov, A superintegrable finite oscillator in two dimensions with SU() symmetry, J. Phys. A: Math. Theor. 46, 1507 (013);

13 13 8] J. Gaboriaud, V. X. Genest, J. Lemieux and L. Vinet, A superintegrable discrete oscillator and two-variable Meixner polynomials, J. Phys. A: Math. Theor. 48, 4150 (015); 9] V. X. Genest, L. Vinet and A. Zhedanov, The multivariate Krawtchouk polynomials as matrix elements of the rotation group representations on oscillator states, J. Phys. A: Math. Theor. 46, (013); 10] V. X. Genest, H. Miki, L. Vinet and A. Zhedanov, The multivariate Meixner polynomials as matrix elements of SO(d, 1) representations on oscillator states, J. Phys. A: Math. Theor (014); 11] V. X. Genest, H. Miki, L. Vinet and A. Zhedanov, The multivariate Charlier polynomials as matrix elements of the Euclidean group representation on oscillator states, J. Phys. A: Math. Theor. 47, 1504 (014); 1] R. Koekoek, P. A. Lesky and R. F. Swarttow, Hypergeometric orthogonal polynomials and their q-analogues (Springer, 010)