Discrete Quantum Harmonic Oscillator

In this paper, we propose a discrete model for the quantum harmonic oscillator. The eigenfunctions and eigenvalues for the corresponding Schrödinger equation are obtained through the factorization method. It is shown that this problem is also connected with the equation for Meixner polynomials.


Introduction
The aim of this paper is to propose a discrete version of the stationary Schrödinger equation for the quantum harmonic oscillator and to solve the associated spectral problem, i.e., to find the eigenfunctions and eigenvalues for the corresponding Hamiltonian. There are in the literature several ways for introducing such a discretization (see, for example, [1][2][3]). The most obvious one consists of simply substituting in place of the quadratic potential V(x) = x 2 of the continuous case the function V(n) = n 2 in the standard central difference expression for the corresponding Hamiltonian. In general, it is not simple to find explicit expressions for the corresponding eigenfunctions and eigenvalues. Moreover, even if we would find them, it is not always true that the spectrum in the discrete case will satisfy the fundamental property of the continuous case, namely that the associated energy levels form a discrete set of equidistant eigenvalues. The discrete version of the infinite square well and the relationship with the tight-binding model was studied by Boykin and Klimeck [4]. In [5], some properties of the mapping of the discrete Schrödinger equation into a two-term wave evolution equation have been considered. In [6,7], the authors showed that supersymmetric quantum mechanics is a simple, powerful tool for generating potentials with known spectra starting from a given initial solvable one. In [8], the factorization method has been used to deal with second-order difference equation on time scales. The construction of creation and annihilation operators for the harmonic oscillator and hydrogen atom on the lattice starting from Rodrigues formula was done by Lorente [2]. The exact discretization of the Schrödinger equation was proposed in [9] based on Fourier transforms.
In a recent paper, we introduced a discrete version of Darboux transformations and Crum formula for generating new potentials starting from an initial solvable one [10], and included a different proposal for the discrete harmonic oscillator potential given by V(n) = n 2 + n + 1, which has an explicit solution for zero energy. Motivated by this, here we use such an alternative choice of discretization and we solve the corresponding discrete stationary Schrödinger equation. We show that the associated spectrum consists of an infinite set of equidistant energy levels, as it happens in the continuous case.
To do that, we employ the factorization method for generating a sequence of Hamiltonians whose eigenfunctions are interrelated through first-order difference operators. We also show that the solution to the eigenvalue problem is connected with the difference equation for the Meixner polynomials.

Preliminaries
In this section, we recall some notation from difference calculus and difference equations, and well-known facts about classical discrete orthogonal polynomials.
We denote by (Z; R) and (Z; C) the sets of real-and complex-valued sequences, respectively. We use the standard notation for shift operators as well as forward and backward difference operators acting on (Z; R) and (Z; C) by Thus, the forward second difference is given by A linear second-order difference equation can be written in general as follows where {a}, {b} and {c} are given sequences. It is well-known how to solve the above equation when the coefficients are constant (see, e.g., [11]).
We want to propose here a discrete version of the Schrödinger equation for the harmonic oscillator There are different approaches for discretizing the one-dimensional time-independent Schrödinger equation (see [4,5,9]). Very often, a discretization appears from the standard central difference formula −∆∇ + n 2 ψ(n) = λψ(n).
In our considerations below, we modify the above equation as follows −∆ 2 + (n 2 + n + 1) ψ(n) = λψ(n + 1), which can be written explicitly as i.e., we consider one step forward shifted eigenfunctions. In addition, there is a relationship of Equation (1) with classical orthogonal polynomials of a discrete variable. Indeed, it is well-known that the equation of hypergeometric type describes such discrete classical orthogonal polynomials as Charlier, Meixner, Kravchuk, and Hahn polynomials, among others (see [12,13]). The sequences σ(n) and τ(n) are second-and first-degree polynomials, respectively, and If σ(n)ρ(n)n i n=a,b = 0, for all i ∈ N ∪ {0}, and the function ρ satisfies the analogue of Pearson equation then the polynomial solutions of Equation (2) are orthogonal with respect to the weight function ρ In particular, the Meixner polynomials are important in our research [14,15]. The Meixner polynomials M γ,µ m (n) are obtained from the following sequences and constants In this case, the scalar product is defined on the interval [a, b] = [0, ∞). These polynomials can be also defined through a generating function in the way The explicit formulas for the Meixner polynomials are the following They can be also expressed in terms of the hypergeometric series m F n a 1 , . . . , a m , b 1 , . . . , b n x := where the Pochhammer symbol is defined by and (a) 0 = 1.

Sequence of Discrete Quantum Harmonic Oscillators
Let k (Z; C), k ∈ Z, be the sets of complex-valued sequences {ψ k } defined on Z and square-summable ψ k |ψ k k < +∞, with the scalar products where ρ k are weight functions. The main object of our considerations is the sequence of second-order operators H k : k (Z, C) −→ k (Z, C) given by H k = −T + + n(n − 1 − k)T − + 2 − k and the problem of eigenvalues for these operators Here, m labels the energy levels and associated eigenfunctions. In this sequence, we recognize the problem of eigenvalues for the discrete quantum harmonic oscillator (Equation (1)) for k = 0.
The family of operators H k can be factorized as a product of first-order operators where the annihilation and creation operators A k : k (Z, C) −→ k−1 (Z, C) and A * k : k−1 (Z, C) −→ k (Z, C) are given by In turn, the constants α k are given by This type of factorizations was presented in detail for the general situation in a discrete case in [16][17][18] and in [19][20][21] for the τ-and q-cases too. It is based on classical methods taken from the work of some founders of quantum mechanics such as Schrödinger [22] (see also [23][24][25]).
In general, the creation operator A * k is the adjoint of the annihilation operator A k relative to the scalar product (6), where the weight functions ρ k satisfy the Pearson type equations (see [17]). In our case, we have and the recursion relation The boundary conditions read as follows If ρ k were a positive weight function, then the operator H k would be self-adjoint and its eigenfunctions corresponding to different eigenvalues would be orthogonal From above expression, after a direct calculation, we find We see that the function ρ k is not positive, i.e., ·|· k is not a scalar product, but it can be interpreted as a pseudo-scalar product. In this case, an explicit calculation using Equations (8) and (9) also yields Let us note that there exists an important class of solutions of Equation (7) such that λ 0 k = α k . If we consider the first-order difference equation then any nonzero solution to this equation is automatically a solution of Equation (7) with eigenvalue It is easy to find now a solution to Equation (14) (the so-called ground state) where C k is a constant. If k ≥ 0, then the series ψ 0 k |ψ 0 k k is convergent. The operators T ± operate on the level of ground states for different k as follows: . . . ψ 0 Moreover, the function ψ 0 k (n) can be used to construct the excited state solutions of Equation (7) through the formula or, equivalently, ψ m k (n) = (nT − − I) m ψ 0 k (n + m). These eigenfunctions correspond to the eigenvalues λ m k = α k−m = 2 − 2k + 2m. For example, a straightforward calculation leads to Moreover, the conditions in Equations (10)-(12) mean that we restrict the sequences {ψ m k } to the spaces k ({k + 1, k + 2, . . . }, C).
These expressions suggest to look for solutions of the eigenvalue problem in Equation (7) for a fixed k in the form ψ m k (n) = P m k (n)ψ 0 k (n).
Substituting this expression into the difference equation (Equation (7)), after some calculations, we get the equation for the sequence P m k (n): We recognize this as the equation for the Meixner polynomials (see [26]). Comparing Equations (16) and (2), we identify Finally, from the comparison of Equations (17) To summarize the previous treatment, the eigenfunctions of the operators H k are given by families of polynomials, which are indexed by the parameter k, with ground states such that From Equation (13) and restricting to the subspace k ({k + 1, k + 2, . . . }, C), we only get that because the function ψ 0 k is not defined at the point k. This means that we have the following relationship

Conclusions
In this work, we propose a discrete version of the harmonic oscillator, which has a spectrum analogous to the corresponding problem in the continuous case. Our construction is based on a factorization method applied to second-order difference equations. In this case, the eigenfunctions turn out to be associated with the Meixner polynomials.
Author Contributions: All authors contributed equally to this research and to writing the paper.