Multiple Meixner polynomials and non-hermitian oscillator Hamiltonians arxiv:1310.0982v2 [math.ca] 8 Nov 2013 F Ndayiragie 1 and W Van Assche 2 1 Département de Mathématiques, Université du Burundi, Campus Mutanga, P.O. BOX 2700, Buumbura, Burundi 2 Department of Mathematics, KU Leuven, Celestinenlaan 200B box 2400, BE-3001 Leuven, Belgium E-mail: ndayiragiefrancois@yahoo.fr E-mail: walter@wis.kuleuven.be Abstract. Multiple Meixner polynomials are polynomials in one variable which satisfy orthogonality relations with respect to r > 1 different negative binomial distributions Pascal distributions. There are two kinds of multiple Meixner polynomials, depending on the selection of the parameters in the negative binomial distribution. We recall their definition and some formulas and give generating functions and explicit expressions for the coefficients in the nearest neighbor recurrence relation. Following a recent construction of Miki, Tsuimoto, Vinet and Zhedanov for multiple Meixner polynomials of the first kind, we construct r > 1 non-hermitian oscillator Hamiltonians in r dimensions which are simultaneously diagonalizable and for which the common eigenstates are expressed in terms of multiple Meixner polynomials of the second kind. PACS numbers: 02.30.Gp, 03.65.Ge AMS classification scheme numbers: 42C05, 33C47, 81R05 Submitted to: J. Phys. A: Math. Theor. 1. Introduction Meixner polynomials M n x;β,c, n = 0,1,2,... are orthogonal polynomials for the negative binomial distribution Pascal distribution, i.e., M n k;β,cm m k;β,c β k c k = 0, n m, k! where β > 0and 0 < c < 1. These polynomialsareorthogonalon the positive integers. In Chihara [2, Chap.VI, 3] they are called Meixner polynomials of the first kind, but nowadays they are more commonly known as Meixner polynomials see, e.g., [4, 6.1] and [5, 9.10] and the Meixner polynomials of the second kind in [2, Chap.VI,
Multiple Meixner polynomials and non-hermitian oscillator Hamiltonians 2 3] are now known as Meixner-Pollaczek polynomials see, e.g., [4, 5.9] and [5, 9.7]. The Meixner polynomials are explicitly given by M n x;β,c = 2 F 1 n, x β 1 1 = c n n xk k β k 1 c so that the leading coefficient is 1 1/c n /β n. The monic polynomials are given by M n x;β,c = β n c n /c 1 n Mn x;β,c. We will investigate multiple Meixner polynomials, which are polynomials of one variable with orthogonality properties with respect to more than one negative binomial distribution. These multiple Meixner polynomials will be used in Section 4 to construct non-hermitian oscillator Hamiltonians in r dimensions which are simultaneously diagonalizable and for which the eigenstates are expressed in terms of multiple Meixner polynomials. This construction was first done in [7] for multiple Charlier polynomials and in [8] for multiple Meixner polynomials of the first kind, but there is a gap in their construction which we will fix in Section 4.1. We extend this construction to multiple Meixner polynomials of the second kind in Section 4.2. The multiple Meixner polynomials hence give a mathematical approach to a physical phenomenon involving a chain of Hamiltonians. There are two kinds of multiple Meixner polynomials, depending on how one chooses the negative binomial distributions [1], [4]. Let n = n 1,...,n r be a multiindex of size n = n 1 + +n r, then the multiple Meixner polynomials of the first kind M 1 n are the monic polynomials of degree n for which M 1 n kklβ kc k k! c k = 0, l = 0,1,...,n 1, = 1...,r, where β > 0 and 0 < c i c < 1 whenever i. The multiple Meixner polynomials of the second kind M 2 n are the monic polynomials of degree n for which M 2 n kklβ k c k k! = 0, l = 0,1,...,n 1, = 1...,r, where 0 < c < 1 and β > 0 for all = 1,...,r, and β i β / Z when i. For r = 1 we retrieve the Meixner polynomials. An explicit formula for the multiple Meixner polynomials of the first kind can be obtained using the Rodrigues formula [1, 4], M 1 n x = β n c n ΓβΓx+1 c 1 Γβ +x n c x Γ n +β +x Γ n +βγx+1 1 where is the backward difference operator, given by fx = fx fx 1. For this backward difference operator we have the property n fx = n c x n 1 k fx k. 2 k
Multiple Meixner polynomials and non-hermitian oscillator Hamiltonians 3 Combining the Rodrigues formula 1 and 2 we have M 1 n x = M 1 n x = n r c k1 1...c kr r 1k1++kr c 1 1...c r 1 Γx+1 Γx k +1 Γβ + n +x k Γβ +x which gives an explicit expression for the multiple Meixner polynomials of the first kind n r x k c n k c 1 nβ +x n k. 3 These multiple Meixner polynomials are related to the multiple Charlier polynomials by the asymptotic relation lim β M1 n x;β,c 1 = a 1 /β,,c r = a r /β = C n x;a 1,...,a r, where the multiple Charlier polynomial C n is the monic polynomial of degree n satisfying C n kk lak k! = 0, l = 0,1,...,n 1, = 1,...,r, where a 1,...,a r > 0 and a i a whenever i see [1], [4]. An explicit expression for the multiple Meixner polynomials of the second kind can be obtained using the Rodrigues formula [1, 4], n M 2 c n x = Γx+1 β n c 1 c x Γβ Γβ +n +x Γβ +x n Γβ +n k=1 Now we can again use property 2 to find M 2 n x = n 1 n r β k nk Γβ 1 Γβ 1 +x Γβ r Γβ r +x k 2 1 cx Γx+1. 4 c n k k Γx+1 c 1 n 1 Γ x+1 n Γβ 1 +n 1 +x Γβ 1 +n 1 Γβ 2 Γβ 2 +x Γβ 2 +n 2 +x k 2 Γβ 2 +n 2 Γβ r +n r +x Γβ r +n r which gives an explicit expression for the multiple Meixner polynomials of the second kind M 2 n x = n r c n k 1 c 1 n x k β +x k i. n k 5
Multiple Meixner polynomials and non-hermitian oscillator Hamiltonians 4 These multiple Meixner polynomials are related to the multiple Charlier polynomials by the asymptotic relation lim β M2 n x;β 1 = a 1 β,...,β r = a r β,c = 1/β = C n x;a 1,...,a r. We will use the explicit expressions 3 and 5 for obtaining the generating functions in Section 2, for computing the recurrence coefficients in Section 3, and for proving the square summability in the appendix. 2. Generating function 2.1. Multiple Meixner polynomials of the first kind Meixner polynomials have the generating function [4, Eq. 6.1.8], [5, Eq. 9.10.11] M n x;β,c β n t n = 1 t x 1 t n! c x β. 6 For multiple Meixner polynomials of the first kind one has a multivariate generating function with r variables. It was found in [8, 9] where it was proved using the Bargmann realization on the Lie algebra Wr made out of r copies of the Heisenberg- Weyl algebra see also Section 4. Theorem 1. Multiple Meixner polynomials of the first kind have the following multivariate generating function = M 1 1 t r n xt n 1!n r! 1+ t 1 ++ t x r 1+ c 1 t 1 ++ c x β r t r. 7 1 c 1 1 c r 1 c 1 1 c r We will use the following lemma, which is basically an alternative way to express the multinomial theorem [9, Eq. 220 on p. 329]. Lemma 1. The generating function for the multinomial coefficients is x n n 1!n r! s 1 s r = 1 s 1 s r x. 8 This series converges absolutely and uniformly for s 1 + + s r < 1 when x / N and contains a finite number of terms if x N. Proof. The use of the binomial theorem and the multinomial theorem leads to x 1 s 1 s r x = 1 n s 1 ++s r n n n x = s 1 n 1,,n s r 1 n r n = n r x n n 1!n r! s 1 s r.
Multiple Meixner polynomials and non-hermitian oscillator Hamiltonians 5 Proof of Theorem 1. Using 3 the multivariate generating function is n 1 n r Changing the order of summation gives n 1= n r= x k!n 1!!n r! c k1 1 c kr r c 1 1 c r 1 β +x n k t 1 t r. x k!n 1!!n r! c k1 1 c kr r c 1 1 c r 1 β +x n k t 1 t r, and by setting l i = n i k i and putting the factors in l i and k i together x k t1!! c 1 1 l 1=0 k1 kr tr c r 1 l r=0 Now we use Lemma 1 to obtain the desired result. 2.2. Multiple Meixner polynomials of the second kind β +x l1 lr l c1 l 1!l r! c 1 1 t cr 1 c r 1 t r. For multiple Meixner polynomials of the second kind one has the following multivariate generating function with r variables. Theorem 2. Multiple Meixner polynomials of the second kind have the following multivariate generating function 1 t r n xt n 1!n r! = 1 1 [ 1 c M 2 1+ c ] 1 c t x 1+ c x β 1 c t. 9 Proof. If we denote the left hand side in 9 by G, then using 5 the multivariate generating function is G = n 1 n r x k c n k!n 1!!n r!c 1 n β 1 +x β 2 +x n2 k 2 β r +x 1 t 1 t r
Multiple Meixner polynomials and non-hermitian oscillator Hamiltonians 6 Changing the order of summation gives G = n 1= n r= x k c n k!n 1!!n r!c 1 n β 1+x β 2 +x n2 k 2 β r +x 1 t 1 t r, and by setting l i = n i k i and putting the factors in l i and k i together G = l 2=0 β 2 +x l2 l 2! x k!! c c 1 t 2 t1 c 1 l2 k1 kr tr c 1 l 1=0 l r=0 β 1 +x l1 l 1! β 1 +x 1 lr l r! For each sum involving the l i we can use the binomial theorem to find G = t2 c 1 1+ c 1 c t 3 l=0 x k!! If we let a = c/1 c and then this simplifies to G = β l t l = 1 t β, t < 1, l! t1 1+ c c 1 1+ c 1 c t r 1 c t 2 k2 tr 1 c 1 1+ c k1 1 c t r l1 c c 1 t 1 c c 1 t r lr. 1+ c kr 1 1 c t tr r c 1 1+ c x β 1 c t. s = t 1+at +1 1+at t, 1 r, x k!! sk1 1 skr r and by using Lemma 1 we then get G = 1+ a c s 1 ++s r a c x k 1+ c x β 1 c t, 1+ c x β 1 c t. kr
Multiple Meixner polynomials and non-hermitian oscillator Hamiltonians 7 The result then follows by observing that s 1 +s 2 ++s r = t 1 1+at 2 1+at r +t 2 1+at 3 1+at r ++t r 1 1+at r +t r = t 1 +t 2 ++t r +a t i t +a 2 t i t t k ++a r 1 t 1 t 2t r i< i<<k = 1 [ 1 1+at 1 ++t r +a 2 t i t +a ] 3 t i t t k ++a r t 1t r a i< = 1 a which implies 1 1+at 1 1+at 2 1+at r, 1+ a c s 1 +s 2 ++s r = 1 1 c [ 1 i<<k 1+ c ] 1 c t. 3. Recurrence relations For multiple orthogonal polynomials there are always nearest neighbor recurrence relations of the form xp n x = P n+ ek x+b n,k P n x+ a n, P n e x, 10 where k {1,...,r} [4, Thm. 23.1.11], [10], and e k = 0,...,0,1,0,...,0 is the kth unit vector in N r. The recurrence relations for multiple Meixner polynomials of the first and second kind were given in [3] without proof. Here we will work out the details of the computations. 3.1. Multiple Meixner polynomials of the first kind Theorem 3. The nearest neighbor recurrence relations for multiple Meixner polynomials of the first kind are xm 1 n x = M1 n+ e k x+ Proof. From 3 and c k β + n + 1 c k + n i 1 c i M 1 n x c n + n 1M1 1 c 2β n e x. 11 x k β +x n k = 1 k x n + 1 k 1 k + 1 k β+ 2β + n k n k 1 x n 1 + 2 2
Multiple Meixner polynomials and non-hermitian oscillator Hamiltonians 8 we find that where δ n = c c 1 n c k 1 k M 1 n x = x n +δ n x n 1 +, n r k 2 2 + k 2 + 2β +2β + n k n k 1, 2 which leads, with some calculations using the binomial theorem, to δ n = n 2 + n 2β 1 n i +β + n 1 2 c i 1. If we compare the coefficients of x n in 10, then b n,k = δ n δ n+ ek, which for the multiple Meixner polynomials of the first kind gives c k n i b n,k = β + n +. 12 1 c c i For the recurrence coefficients a n, we can use [10, Thm. 3.2, Eq. 3.12] Using 12 we find that a n, a n+ ei, = b n e, b n e,i b n, b n,i. a n, = β + n 1 β + n 2 a n e i,. If we use this formula recursively to go from n to n e by decreasing one of the indices except the th index by one, then we arrive at a n, = β + n 1 β +n 1 a n e,, and a n e, is the recurrence coefficient a 2 n for the Meixner polynomials M n x;β,c, c which is equal to n 1 c β +n 2 1. This gives a n, = c n 1 c 2β + n 1. 3.2. Multiple Meixner polynomials of the second kind Theorem 4. The nearest neighbor recurrence relations for multiple Meixner polynomials of the second kind are xm 2 c n x = M2 n+ e k x+ 1 c n k +β k + n M 2 n 1 c x + cn n +β 1 1 c 2 i n +β β i n +β n i β i M 2 n e x. 13
Multiple Meixner polynomials and non-hermitian oscillator Hamiltonians 9 Proof. From 5 and 1 x k β +x k i n k = 1 k x n + 1 k x n 1 k 2 we find that where δ n = n 1 n r { 1 k k 2 2 2 + k 2 +1 2 1 2β 2 k i +n k 1 n k +, M 2 n x = x n +δ n x n 1 +, + k 2 + 1 2 c n k c 1 n 1 } 2β 2 k i +n k 1 n k, which leads, with some calculations using the formula =2 and the binomial theorem, to δ n = 1 c 1 k 1 k i = k 2 2 1 2 1 n n +β 1+ n i + k 2, n 2 2 + n 2 2β 1. If we compare the coefficients of x n in 10, then b n,k = δ n δ n+ ek, which for the multiple Meixner polynomials of the second kind gives b n,k = c 1 c n k +β k + n 1 c. 14 For the recurrence coefficients a n, we can again use [10, Thm. 3.2, Eq. 3.12] a n, a n+ ei, = b n e, b n e,i b n, b n,i. Using 14 we find that a n, = n +β n i β i +1 n +β n i β i a n ei,. Using this formula recursively to go from n to n e by decreasing each of the indices except the th index by one, we find a n, =,i n +β β i n +β n i β i a n e,,
Multiple Meixner polynomials and non-hermitian oscillator Hamiltonians 10 where a n e, is the recurrence coefficients a 2 n for the Meixner polynomials cn M n x;β,c, which is equal to 1 c n 2 +β 1. This gives a n, = cn n +β 1 1 c 2 i n +β β i n +β n i β i. 4. Non-Hermitian oscillator Hamiltonians In [7] Miki et al. showed how a set of r non-hermitian oscillator Hamiltonians can be obtained which are simultaneously diagonalizable, with real spectra and common eigenstates which are expressed in terms of multiple Charlier polynomials. This idea was further extended in [8] to construct r non-hermitian oscillator Hamiltonians in r dimensions which are simultaneously diagonalizable and for which the spectra are real and the common eigenstates are expressed in terms of multiple Meixner polynomials of the first kind. Here we will construct r non-hermitian oscillator Hamiltonians for which the common eigenstates are expressed in terms of multiple Meixner polynomials of the second kind. The construction uses a Hilbert space H with an orthonormal basis { n,n = 1,2,3,...} such that n m = δ n,m, and on this Hilbert space one has two operators b and b + with the property { b n = n n 1, annihilation operator, b + n = n+1 n+1, creation operator. The ground state 0 has the property b 0 = 0 and all the states n are obtained from the ground state by repeated application of the operator b + in the sense that n = b + n 0 / n!. The annihilation and creation operators b and b + satisfy the commutation relations [b,b] = [b +,b + ] = 0, [b,b + ] = 1, and generate the Lie algebra W with generators {b,b +,1}. If we take r copies W 1,...,W r then we can construct the algebra Wr = where the annihilating operators b 1,...,b r are commuting and the creation operators b + 1,...,b+ r are also commuting, and r W i [b i,b + ] = δ i,. We denote by n 1,...,n r the standard basis in the r-dimensional Hilbert space H H, and we have b i n 1,...,n i,...,n r = n i n 1,...,n i 1,...,n r,
Multiple Meixner polynomials and non-hermitian oscillator Hamiltonians 11 Observe that b + i n 1,...,n i,...,n r = n i +1 n 1,...,n i +1,...,n r. where n i N. b + i b i n 1,...,n i,...,n r = n i n 1...,n i,...,n r, 1 i r, 4.1. Multiple Meixner polynomials of the first kind Miki et al. considered in [8] the Hamiltonians b + k H i = b i + b k + c i β + b + k 1 c c b k + β + i k=1 k=1 for = 1,...,r, and they showed that the states where N x 1 H 1,...,H r : x = N 1 x b + k b c k 1 c 2b+ k=1 M 1 n x!n r! n 1,...,n r, are normalizing constants, are common eigenstates for the operators H i x = x x, 1 i r. The normalizing constants are such that x x = 1, hence [N 1 x ]2 = [M 1 1 n x]2, n 1!n r! however Miki et al. did not show that the infinite sum is finite. In fact this sum diverges for every x R, as we show in the appendix. We propose to renormalize the standard basis n 1,...,n r and to take their norm to be where t 1,...,t r R r + satisfies n 1,...,n r n 1,...,n r = t 1 t r, β n c 2 1 c 2t < 1, and at least one of the t i 1 i r satisfies t i > 1 c i 2. For instance, if c 1 ++c r < 1 then one can take t = 1 c 2 /c for all {1,2,...,r}. Then the norm of the eigenstate in the weighted r-dimensional Hilbert space is [ ˆN 1 x ] 1 = [M 1 1/2 n x]2 t 1 n 1!n r!β t r n and this is finite if and only if x N = {0,1,2,...} see Theorem 5 in the appendix. This way there are countably many eigenstates for the eigenvalues x N.
Multiple Meixner polynomials and non-hermitian oscillator Hamiltonians 12 4.2. Multiple Meixner polynomials of the second kind We will now construct Hamiltonians for which the common eigenstates are expressed in terms of multiple Meixner polynomials of the second kind. This is a new construction but it follows similar ideas as in [8]. Let 0 < c < 1 and β 1,...,β r > 0 be such that β i β / Z whenever i. We first introduce r operators B by B = b + b +β b + l b l β l, 1 r, l=1,l which act on the basis n 1,...,n r as B n 1,...,n r = l=1,l n +β n l β l n 1,...,n r. Hence the n 1,...,n r are eigenstates and the corresponding eigenvalues are nonzero because β β l / Z whenever l. Hence B 1,...,B r are invertible. We define the Hamiltonians H i = b i + c 1 c b+ i b i +β i + 1 1 c + c 1 c 2 and the states b + b +β x = N 2 x k=1,k b + b b + b +β β k B 1 b +, 1 i r, 15 M 2 n x!n r! n 1,...,n r. Note that each H i contains an annihilation operator b i and a number operator b + i b i, but the remaining terms are independent of i and involve the total number operator r b+ b and all the creation operators b + 1,...,b+ r. The interpretation is that there are r different kinds of particles and H i only annihilates particles of the ith kind and creates particles of all kinds, independent of i but with a certain interaction which can be tuned using the parameters β 1,...,β r and c. We normalize the standard basis such that t 1 n 1,...,n r n 1,...,n r = t r, β 1 β r where t 1,...,t r R r + satisfies t < 1 c2 c 2, {1,2,...,r} and at least one of the t i satisfies t i > 1 c 2. One can take, for instance, t = 1 c 2 /c for all {1,2,...,r}. Then the normalizing constants satisfy [N x 2 [M 2 1 ]2 n = x]2 t 1 n 1!n r!β 1 β r t r. It is shown in the appendix Theorem 6 that this multiple sum is finite if and only if x N = {0,1,2,...}, hence this gives a countable number of eigenstates in the weighted r-dimensional Hilbert space.
Multiple Meixner polynomials and non-hermitian oscillator Hamiltonians 13 5. Further work We have given some properties of multiple Meixner polynomials and have shown how these polynomials can be used to construct non-hermitian oscillator Hamiltonians in r dimensions with a countable number of eigenvalues and eigenstates. Other families of multiple orthogonal polynomials, such as multiple Hahn polynomials or multiple Krawtchouk polynomials, may also enable the construction of Hamiltonians for which the eigenstates can explicitly be given. Acknowledgments We thank the referees for suggesting a number of improvements. This research was supported by FWO grant G.0934.13, KU Leuven research grant OT/12/073 and the Belgian Interuniversity Attraction Pole P7/18. F.N. is supported by a VLIR/UOS scholarship. This work was initiated in the Ph.D. dissertation of F.N. at KU Leuven and he thanks the Department of Mathematics at KU Leuven for their hospitality. Appendix If p n n N is a sequence of orthonormal polynomials on the real line for a positive measure µ for which the moment problem is determinate, p n xp m xdµx = δ m,n, then it is well known that R [p n x] 2 1 = µ{x}, A.1 so that this infinite sum is finite if and only if x is a mass point of the measure µ, i.e., x belongs to the support of the discrete part µ d of the measure µ [4, Theorem 2.5.6]. The orthogonality measure for Meixner polynomials is the negative binomial distribution or Pascal distribution which is a discrete measure supported on the integers N = {0,1,2,...}. Following [5, 9.10] the Meixner polynomials are given by M n x;β,c = 2 F 1 n, x β ;1 1 c and then the orthogonality for the Meixner polynomials is 1 c β M n k;β,c M m k;β,c β kc k k!, A.2 = n! c n β n δ m,n. A.3 The orthonormal polynomials p n x and the monic polynomials M n x are given by p n x = M n x;β,c so that A.1 becomes 1 c β βn c n n! [M n x] 2 n!β n c 1 2n c n = n c, M n x = β n Mn x;β,c, c 1 { k! c k β k, x = k N = {0,1,2,...}, +, x R\N. A.4
Multiple Meixner polynomials and non-hermitian oscillator Hamiltonians 14 Observe that A.2 implies that Mn k;β,c = M k n;β,c, hence the orthogonality A.3 also gives 1 c β M n km n lc 1 2n n!β n c n = k! c k δ k,l, β k A.5 for which the case k = l corresponds to A.4. In fact more is true: the generating function 6 for the monic Meixner polynomials is M n x tn n! = 1+ t x 1+ c x β 1 c 1 c t. When x / N we see that the singularity closest to the origin is at t = 1 c, hence the radius of convergence is R = 1 c. This means that and hence lim sup n Mn x 2 n!β n lim sup n Mn x 1/n = limsup n n! 1/n = 1 1 c, Mn x 2 1/n 1/n n! n! 2 = β n 1 1 c 2, so that for x / N [M n x] 2 n!β n t n { < if 0 t < 1 c 2, = if t > 1 c 2. A.6 The case t = 1 c 2 /c > 1 c 2 corresponds to the divergent part of A.4. When x N the first singularity in the generating function is at t = 1 c/c, hence and therefore This means that for x N lim sup n lim sup n [M n x] 2 n!β n t n Mn x n! Mn x 2 1/n = n!β n 1/n = c 1 c, c 2 1 c 2. { < if 0 t < 1 c 2 /c 2, = if t > 1 c 2 /c 2. A.7. The case t = 1 c 2 /c < 1 c 2 /c 2 corresponds to the convergent part of A.4. We now prove two similar results for the multiple Meixner polynomials, which are needed in Section 4 Theorem 5. Let A 1 = {t 1,...,t r R r + : c1 1 c 1 2 2 cr t 1 ++ t r < 1}, 1 c r
Multiple Meixner polynomials and non-hermitian oscillator Hamiltonians 15 and B i = {t 1,...,t r R r + : t i > 1 c i 2 }. For multiple Meixner polynomials of the first kind one has for t 1,...,t r A 1 r B i [M 1 n x]2 t 1 n 1!n r!β t r n { <, if x N = {0,1,2,...}, = if x R\N. Proof. Suppose that x R\N and t i > 1 c i 2, then [M 1 n x]2 n 1!n r!β n t 1 t r n i=0 [M 1 0,...,0,n i,0,...,0 x]2 n i!β ni and this single sum divergessincem 1 0,...,0,n i,0,...,0 x isthe monic Meixnerpolynomial M n x;β,c i, see A.6. This shows the divergent part of the theorem. Suppose that x N, say x = k, then the multiple sum 3 only contains a finite number of terms since k k = 0 whenever k > k. This means that M 1 n k = c n c 1 n k k k k β +k n k t ni i, c k. A simple estimation, using k k k! and β+k n k = β n β+ n k k /β k gives M 1 n k c n 1 c nk!β n = k k c n 1 c nk!β n Pn 1,...,n r, β + n k k β k c k where Pn 1,...,n r is a polynomial in the variables n 1,...,n r. This means that k! 2 [M 1 n k]2 t 1 n 1!n r!β t r n [Pn 1,...,n r ] 2 c 2 1 cr 2 β n 1 c 1 2 1 c r 2 n 1!n r! t 1 t r. The latter sum converges for t 1,...,t r A 1 see Lemma 1 and Abel s lemma [6, Prop. 2.3.4 on p. 87]. This proves the convergent part of the theorem. Theorem 6. Let 2 1 c A 2 = {t 1,...,t r R r + : t i <, 1 i r} c and B i = {t 1,...,t r R r + : t i > 1 c 2 }.
Multiple Meixner polynomials and non-hermitian oscillator Hamiltonians 16 For multiple Meixner polynomials of the second kind one has for t 1,...,t r A 2 r B i [M 2 n x]2 t 1 n 1!n r!β 1 β r t r Proof. Suppose that x R\N and t i > 1 c 2, then [M 2 n x]2 n 1!n r!β 1 β r t 1 t r { <, if x N = {0,1,2,...}, =, if x R\N. n i=0 [M 2 0,...,0,n i,0,...,0 x]2 n i!β i ni and this single sum divergessincem 2 0,...,0,n i,0,...,0 x isthe monic Meixnerpolynomial M n x;β i,c, see A.6. This shows the divergent part of the theorem. Suppose that x N, say x = k, then the multiple sum 5 only contains a finite number of terms since k k = 0 whenever k > k. This means that M 2 n k = c n c 1 n k k k k c k t ni i, 1 β +k k i. n k A simple estimation, using k k k! and β + k l n k = β n β + n k k l/β k l, gives M 2 n k c n 1 c n k! β n k k c n = 1 c n k! β n Qn 1,...,n r, c k β +n k ki β k 1 ki where Qn 1,...,n r is a polynomial in the variables n 1,...,n r. This means that [M 2 n k]2 t 1 n 1!n r!β 1 β r t r k! 2 [Qn 1,...,n r ] 2 c 2 n β 1 β r 1 c 2 n n 1!n r! t 1 t r. The latter sum converges for t 1,...,t r A 2 since the multiple sum without the polynomial Q factors into a product of r binomial series and the multiple sum with the polynomial Q still convergesby Abel s lemma [6, Prop. 2.3.4on p. 87]. This proves the convergent part of the theorem. Similar results are true for Charlier polynomials and multiple Charlier polynomials. These were considered in [7]. The Charlier polynomials C n a are orthogonal polynomials for the Poisson distribution e a C a n kca m kak k! = an n!δ m,n,
Multiple Meixner polynomials and non-hermitian oscillator Hamiltonians 17 where a > 0. They are given by C a n x = an 2F 0 n, x 1 a from which we find that C n a k/ a n = C a k n/ ak. Hence the orthogonality also gives C a k nca k m 1 a k k! = ea a n n!δ m,n. The orthonormal Charlier polynomials are p n x = C n a x/ a n n!, so that the last sum for n = m becomes [p k n] 2 = e a a n n!, which indeed corresponds to A.1. In fact more is true: the generating function for Charlier polynomials is [2, VI.1] C a n xtn n! = 1+tx e at. When x / N we see that this function has a singularity at t = 1, hence the radius of convergence is R = 1, which implies that lim sup n C a 1/n n x = 1. n!, But then lim sup n C a n x 2 n! 1/n = limsup n C a n x 2 n! 2 1/n n! 1/n =, so that for x / N [C a n x] 2 t n = n! for every t > 0. When x N the generating function has no singularities and hence the radius of convergence is. One has C n a k = an k =0 n k a so that giving C a n k an k! [C a k =0 n a k!a+1 n, n k] 2 t n k! 2 a+1 2n t n = e a+12 t n! n!
Multiple Meixner polynomials and non-hermitian oscillator Hamiltonians 18 so that this series converges for every t > 0. Hence one has for every t > 0 { [C n a x] 2 t n <, if x N = {0,1,2,...},. n! =, if x R\N. Multiple Charlier polynomials are defined by n 1 C n x = 1 n n r x k1++ a k1 1 a kr r, where a i > 0 for 1 i r and a i a whenever i. The square summability for these multiple Charlier polynomials, needed for the normalization of the eigenfunctions in [7], is given by Theorem 7. For multiple Charlier polynomials one has for all t 1,...,t r R r + { [C n x] 2 <, if x N = {0,1,2,...}, n 1!n r! t 1 t r =, if x R\N. The proof is similar and in fact a bit easier to the proofs before. The choice t 1 = = t r = 1 shows that the eigenstates in [7] are indeed those corresponding to the eigenvalues x = k N. References [1] J. Arvesú, J. Coussement, W. Van Assche, Some discrete multiple orthogonal polynomials, J. Comput. Appl. Math. 153 2003, 19 45. [2] T.S. Chihara, An Introduction to Orthogonal Polynomials, Mathematics and its Applications 13, Gordon and Breach, New York, 1978. [3] M. Haneczok, W. Van Assche, Interlacing properties of zeros of multiple orthogonal polynomials, J. Math. Anal. Appl. 389 2012, 429 438. [4] M.E.H. Ismail, Classical and Quantum Orthogonal Polynomials in One Variable, Encyclopedia of Mathematics and its Applications 98, Cambridge University Press, 2005 paperback edition, 2009. [5] R. Koekoek, P.A. Lesky, R.F. Swarttouw, Hypergeometric Orthogonal Polynomials and Their q-analogues, Springer Monographs in Mathematics, Springer-Verlag, Berlin, 2010. [6] S.G. Krantz, Function Theory of Several Complex Variables, John Wiley & Sons, New York, 1982. [7] H. Miki, L. Vinet, A. Zhedanov, Non-Hermitian oscillator Hamiltonians and multiple Charlier polynomials, Phys. Lett. A 376,. 2 2011, 65 69. [8] H. Miki, S. Tsuimoto, L. Vinet and A. Zhedanov, An algebraic model for the multiple Meixner polynomials of the first kind, J. Phys. A: Math. Theor. 45,. 32 2012, 325205. [9] H.M. Srivastava, P.W. Karlsson, Multiple Gaussian Hypergeometric Series, Ellis Horwood Ltd, Chichester, 1985. [10] W. Van Assche, Nearest neighbor recurrence relations for multiple orthogonal polynomials, J. Approx. Theory 163 2011, 1427 1448.