Submitted Successfully!
To reward your contribution, here is a gift for you: A free trial for our video production service.
Thank you for your contribution! You can also upload a video entry or images related to this topic.
Version Summary Created by Modification Content Size Created at Operation
1 We introduce exact finite differences for the derivative, we developed the calculus of them. To the future, we can find exact finite differences for any function and they can be used in Quantum mechanics, for the determination of time eigenstates + 2609 word(s) 2609 2019-08-07 21:31:41 |
2 format correction -85 word(s) 2524 2019-11-01 06:51:22 | |
3 Roll back entry to make some changes. + 85 word(s) 2609 2020-10-28 07:33:50 |

Video Upload Options

Do you have a full video?


Are you sure to Delete?
If you have any further questions, please contact Encyclopedia Editorial Office.
Torres-Vega, G. A discrete quantum momentum operator. Encyclopedia. Available online: (accessed on 15 June 2024).
Torres-Vega G. A discrete quantum momentum operator. Encyclopedia. Available at: Accessed June 15, 2024.
Torres-Vega, Gabino. "A discrete quantum momentum operator" Encyclopedia, (accessed June 15, 2024).
Torres-Vega, G. (2019, August 07). A discrete quantum momentum operator. In Encyclopedia.
Torres-Vega, Gabino. "A discrete quantum momentum operator." Encyclopedia. Web. 07 August, 2019.
A discrete quantum momentum operator

We introduce finite-differences derivatives intended to be exact when applied to the real exponential function. We want to recover the known results of continuous calculus with our finite differences derivatives but in a discrete form. The purpose of this work is to have a discrete momentum operator suitable for use as an operator in discrete quantum mechanics theory.

discrete derivative discrete symmetric operator, discrete quantum mechanics


The subject of finite differences is an old, but  useful method with wide application in science and engineering. The story starts with Newton and Leibnitz themselves with the very  definition of the derivative of a function; the limit when the differences become zero, as is well known. The formal definition of the derivative \( f'(x_0) \) of some function \( f(x) \) of a single variable \( x \), at \( x_0 \), is obtained from the limit \( x\to x_0 \) of any of the following ratios of differences

      \( \displaystyle f'(x_0) :=\lim_{x\to x^-_0} \frac{f(x_0)-f(x)}{x_0-x}, \quad f'(x_0) :=\lim_{x\to x^+_0} \frac{f(x)-f(x_0)}{x-x_0}. \)

These ratios are known as backward and forward differences, respectively.

The limit of vanishing finite differences involves variables which are continuous. However, many times it is not possible to perform such limit  at all. A situation like that appears when doing numerical simulations in a computer due to the limit in how small a number can be represented in the computer. Another situation in which it is not possible to take the infinitesimal difference limit is when the independent variable is a discrete variable by nature, as is the case found in Quantum Mechanics theory regarding the spectrum of quantum operators.

Finite differences are necessary in numerical computations to approximate several quantities like derivatives and integrals of functions, as well as, differential and integral equations. Their use as a numerical tool is a well developed subject, see for instance the series of books Numerical Recipes.[1] An interesting development is the use of a complex finite differences to evaluate the real derivative of a function.[2]

On the other hand, the reader can learn about the calculus of finite differences  from the classical works of Kopal[3], Boole[4]. Jordan[5], Richardson[6], for instance. 

Another branch of the finite differences tree, is known as the exact finite differences technique. That scheme was developed by researchers like Potts,[7][8] Ronald E. Mickens[9] with the purpose of obtaining exact finite differences representations of continuous differential equations and of their solutions.

Another use of finite differences was developed by Armando Martínez-Pérez and Gabino Torres-Vega[10] with the intention of obtaining discrete operators for use in Quantum Mechanics theory. It is common that a quantum operator has a discrete spectrum and a derivative with respect to the spectrum is necessary some times. The aim is to obtain discrete operators which comply with discrete versions of the properties that a quantum operator must have.

The explicit form of the exact finite difference derivative depends on the functions that we want to consider. In this article, we discuss the use of an exact finite differences derivative when its eigenfunction is the real exponential function. This will provide with a momentum operator to be used in discrete Quantum Mechanics.


The partition

Let \( \{q_{i}\}_{i=-N}^{N} \) be a \( N+1 \) points partition of the real interval \( [-a,a] \), that is

          \( \displaystyle -a=q_{-N} < q_{-N+1} < \dots < q_{N-1} < q_{N}=a. \)

The separations between mesh points are

          \( \displaystyle \Delta_j = q_{j+1}-q_{j}, \quad j = -N, \dots, N-1. \)

These separations can or cannot be equal and are finite. The discrete variable $q_j$ is the independent variable with respect to we will calculate the derivative of a function.

The exact finite differences derivative

For a given \( v\in\mathbb{R} \), backward and forward finite differences  derivatives, of a vector \( g=(g_{-N}, g_{-N+1},\dots,g_{N-1},g_N)^T \), at \( q_j \), defined on the partition, are

         \( \displaystyle (D_bg)_j=\frac{g_{j}-g_{j-1}}{\eta_{j}}, \)

          \( \displaystyle (D_fg)_j=\frac{g_{j+1}-g_j}{\chi_{j}}, \)
where the denominators are defined as

          \( \displaystyle \eta_{j}=\frac{2}{v}\,e^{-v\Delta_{j-1}/2} \sinh\left(\frac{v}{2}\Delta_{j-1}\right) =\Delta_{j-1}+\frac{v}{2}\Delta^2_{j-1}+{\cal O}(\Delta^3_j), \)

          \( \displaystyle \chi_{j}=\frac{2}{v}\,e^{v\Delta_{j}/2} \sinh\left(\frac{v}{2}\Delta_{j}\right) =\Delta_j+\frac{v}{2}\Delta^2_j+{\cal O}(\Delta^3_j). \)

The discrete variable \( q_j \) is the independent variable with respect to which we will calculate the derivative of a function. The functions in the denominators of these expressions, the functions \( \eta_{j} \) and \( \chi_{j} \), makes sure that the real exponential function \( e^{vq} \) be, exactly, an eigenfunction of the discrete derivative operation with real eigenvalue \( v \),

          \( \displaystyle (D_be^{vq})_j=(D_fe^{vq})_j=v\,e^{vq_j}, \)

which is the same property as the continuous variable derivative has.

When \( |\eta_{j}-\Delta_j| \) and \( |\chi_{j}-\Delta_j| \) are less than some small \( \epsilon \), our method and the usual finite-differences derivative give similar results for any vector defined on the mesh. Note that the denominators \( \eta_{j} \) and \( \chi_{j} \) become \( \Delta_j \) when \( v=0 \).


Properties of the finite-differences derivative

Some properties of the exact finite-differences derivative are

The connection between forward and backward discrete derivatives. The equality \( \displaystyle \eta_{j+1}=e^{-v\Delta_j} \chi_{j} \) implies that

          \( (D_fg)_j=e^{-v\Delta_j}\,(D_bg)_{j+1}. \)

The summation of the derivative. The finite differences versions of \( \displaystyle \int_a^xdy\,g'(y)=g(x)-g(a) \) are

\( \displaystyle \sum_{j=-N}^n\Delta_j(D_fg)_j =-\frac{\Delta_{-N}}{\chi_{-N}}g_{-N} +\sum_{j=-N+1}^n\left( \frac{\Delta_{j-1}}{\chi_{j-1}} -\frac{\Delta_j}{\chi_{j}}\right)g_j +\frac{\Delta_n}{\chi_{n}}g_{n+1}, \)

and\( \displaystyle \sum_{j=-N+1}^n\Delta_{j-1}(D_bg)_j =-\frac{\Delta_{-N}}{\eta_{-N+1}}g_{-N} +\sum_{j=-N+1}^{n-1}\left( \frac{\Delta_{j-1}}{\eta_{j}} -\frac{\Delta_j}{\eta_{{j+1}}}\right)g_j +\frac{\Delta_{n-1}}{\eta_{n}}g_n, \)

where \( \)\( n \)<\( N \). The summation term at the right hand side of these equalities involve asymmetry terms that vanish when the separation between mesh points is the same.

The eigenfunction of the summation. The finite differences versions of \( \displaystyle \int_a^xdx\,v\,e^{vx}=e^{vx}-e^{va} \) are\( \displaystyle \sum_{j=-N}^n\Delta_jv\,e^{vq_j} =\sum_{j=-N}^n\Delta_j(D_f\,e^{vq})_j =-\frac{\Delta_{-N}}{\chi_{-N}}e^{vq_{-N}} +\sum_{-N+1}^n\left( \frac{\Delta_{j-1}}{\chi_{j-1}} -\frac{\Delta_j}{\chi_{j}}\right)e^{vq_j} +\frac{\Delta_n}{\chi_{n}}e^{vq_{n+1}}, \)

and\( \displaystyle \sum_{j=-N+1}^n\Delta_{j-1}v\,e^{vq_j} =\sum_{j=-N+1}^n\Delta_{j-1}(D_b\,e^{vq})_j =-\frac{\Delta_{-N}}{\eta_{-N+1}}e^{vq_{-N}} +\sum_{j=-N+1}^{n-1}\left( \frac{\Delta_{j-1}}{\eta_{j}} -\frac{\Delta_j}{\eta_{j+1}}\right)e^{vq_j} +\frac{\Delta_{n-1}}{\eta_{n}}e^{vq_n}, \)where \( -N \)<\( n \)<\( N \).

The derivative of a constant function \( c \)

          \( (D_f\,c)_j=0, \quad\text{and}\quad (D_b\,c)_j=0. \)
The derivatives of \( q \).

\( \displaystyle (D_f\,q)_j =\frac{q_{j+1}-q_j}{\chi_{j}} =\frac{\Delta_j}{\chi_{j}}, \quad\text{and}\quad (D_b\,q)_j =\frac{q_j-q_{j-1}}{\eta_{j}} =\frac{\Delta_{j-1}}{\eta_{j}}. \)

These derivatives will approach to one in the limit of small \( \Delta_j \). In particular, both \( (D_f\,q)_j \) and \( (D_b\,q)_j \) are equal to one when \( v=0 \).

The chain rule. For the forward scheme this rule turns to be

          \( \displaystyle (D_fg(h(q)))_j=({\cal D}_fg(h))_j\,(D_fh)_j, \)

          \( \displaystyle ({\cal D}_fg(h))_j =\frac{g(h(q_{j+1}))-g(h(q_j))}{h(q_{j+1})-h(q_j)}. \)
For the backward finite differences we have

        \( \displaystyle (D_bg(h(q)))_j=({\cal D}_bg(h))_j\,(D_bh)_j, \)

          \( \displaystyle ({\cal D}_bg(h))_j =\frac{g(h(q_j))-g(h(q_{j-1}))}{h(q_{j})-h(q_{j-1})}. \)

The derivative of a product of vectors. There are four equalities

         \( (D_f\,gh)_j=g_{j+1}(D_fh)_j+(D_fg)_jh_j., \)
         \( \displaystyle (D_f\,gh)_j=(D_fg)_jh_{j+1}+g_j(D_fh)_j, \)
         \( \displaystyle (D_b\,gh)_j=g_j(D_bh)_j+(D_bg)_jh_{j-1}, \)
          \( \displaystyle (D_b\,gh)_j=(D_bg)_jh_j+g_{j-1}(D_bh)_j. \)

The derivative of a vector of inverses.

         \( \displaystyle \left(D_f\,\frac{1}{h}\right)_j=-\frac{(D_fh)_j}{h_jh_{j+1}}, \quad\text{and}\quad \left(D_b\,\frac{1}{h}\right)_j=-\frac{(D_bh)_j}{h_jh_{j-1}}, \)

provided \( h_j,h_{j\pm 1} \neq 0 \).

The derivative of a vector of ratios. There are four versions of this property

          \( \displaystyle \left(D_f\,\frac{g}{h}\right)_j=\frac{(D_f\,g)_j}{h_j} -\frac{(D_fh)_j}{h_jh_{j+1}}g_{j+1}, \)

          \( \displaystyle \left(D_f\,\frac{g}{h}\right)_j=\frac{(D_f\,g)_j}{h_{j+1}} -\frac{(D_fh)_j}{h_jh_{j+1}}g_j, \)

          \( \displaystyle \left(D_b\,\frac{g}{h}\right)_j=\frac{(D_b\,g)_j}{h_{j-1}} -\frac{(D_bh)_j}{h_jh_{j-1}}g_j, \)

        \( \displaystyle \left(D_b\,\frac{g}{h}\right)_j=\frac{(D_b\,g)_j}{h_j} -\frac{(D_bh)_j}{h_jh_{j-1}}g_{j-1}. \)

Summation by parts.

          \( \displaystyle \sum_{j=-N}^{N-1}\Delta_jg_{j+1}(D_fh)_j +\sum_{j=-N}^{N-1}\Delta_j(D_fg)_jh_j=\cal{I}_1, \)

where        \( \displaystyle {\cal I}_1 =\sum_{j=-N}^{N-1}\Delta_j(D_f\,gh)_j =-\frac{\Delta_{-N}}{\chi_{-N}}g_{-N}h_{-N} +\sum_{-N+1}^{N-1}\left( \frac{\Delta_{j-1}}{\chi_{j-1}} -\frac{\Delta_j}{\chi_{j}}\right)h_jg_j +\frac{\Delta_{N-1}}{\chi_{{N-1}}}h_Ng_N. \)The terms in the summation vanish when \( \Delta_j\to0 \), i.e. when the separation between the mesh points is the same. But we can chose the vectors \( \mathbf{h} \) and \( \mathbf{g} \) is such a way that the last sum vanishes. Also

          \( \displaystyle \sum_{j=-N+1}^n\Delta_{j-1}g_j(D_bh)_j +\sum_{j=-N+1}^n\Delta_{j-1}(D_bg)_jh_{j-1} =\cal{I}_2, \)
          \( \displaystyle \cal{I}_2 =-\frac{\Delta_{-N}}{\chi_{b,-N+1}}g_{-N}h_{-N} +\sum_{j=-N+1}^{N-1}\left( \frac{\Delta_{j-1}}{\chi_{b,j}} -\frac{\Delta_j}{\chi_{b,{j+1}}}\right)g_jh_j +\frac{\Delta_{N-1}}{\chi_{b,N}}g_Nh_N. \)

The commutator between \( q \) and \( D \). The discrete version of the relationship \( \displaystyle \frac{d}{dq}qh(q)-q\frac{d}{dq}h(q)=h(q) \) becomes

          \( \displaystyle (D_f\,qh)_j-q_{j+1}(D_f\,h)_j=\frac{\Delta_j}{\chi_{j}}q_j, \)
          \( \displaystyle (D_b\,qh)_j-q_j(D_b\,h)_j=\frac{\Delta_j}{\eta_{j}}q_{j-1}. \)

Translation of the exponential function. The discrete derivative is the generator of translations of the exponential function
          \( \displaystyle (e^{s\,D_{b,f}}e^{vq})_j=e^{v(q_j+s)}. \)

We will need of the bounded Fourier transform of a given function \( g(v) \) which is defined as

          \( \displaystyle \tilde g_j=(Fg)_j=\frac{1}{\sqrt{{\cal V}}} \int_{-{\cal V}/2}^{{\cal V}/2}e^{ivq_j}g(v)dv, \)

on the mesh.

We also need of the discrete Fourier transform of a vector \( h_j \) which is defined as
          \( \displaystyle \tilde h(v)=({\cal F}h)(v) =\sum_{j=-N}^{N-1}\sqrt{\frac{\Delta_j}{\cal V}}\,e^{-ivq_j}h_j. \)
These transformations preserve the norm of vectors and functions. In the equidistant case in which \( q_{j+1}-q_j=2\pi/(N-1) \) these transforms can be identified with the usual Fourier transform on \( [-{\cal V},{\cal V}] \)
and with the Fourier series, respectively.

The adjoint of \( D \). Consider the equality          \( \displaystyle v({\cal F}h)_j=\sum_{j=-N}^{N-1}\sqrt{\frac{\Delta_j}{\cal V}}\, e^{-ivq_{j+1}}(-iD_fh)_j +i\sum_{j=-N}^N\sqrt{\frac{\Delta_j}{\cal V}}\, (D_fe^{-ivq}h)_j. \)
Then, there is the relationship
          \( \displaystyle p\leftrightarrow e^{-iv\Delta_j}(-iD_fh)_j,\quad j=-N,\dots,N-1, \)
provided that the asymmetry term
          \( \displaystyle \sum_{j=-N}^{N-1}\sqrt{\frac{\Delta_j}{\cal V}}\, (D_fe^{-ivq}h)_j, \)

The conjugate of \( q_j \). The relationship 
          \( \displaystyle q_j\tilde g_j =i\frac{1}{\sqrt{{\cal V}}} \int_{-{\cal V}/2}^{{\cal V}/2} e^{ivq_j}\frac{d\,g(v)}{dv}dv -i\frac{1}{\sqrt{{\cal V}}} e^{ivq_j}g(v)\Big|_{v=-{\cal V}/2}^{{\cal V}/2} \)
indicates that
          \( \displaystyle q_j\tilde g_j =i\frac{1}{\sqrt{{\cal V}}} \int_{-{\cal V}/2}^{{\cal V}/2} e^{ivq_j}\frac{d\,g(v)}{dv}dv -i\frac{1}{\sqrt{{\cal V}}} e^{ivq_j}g(v)\Big|_{v=-{\cal V}/2}^{{\cal V}/2} \)
provided boundary term
          \( \displaystyle \frac{i}{\sqrt{{\cal V}}} e^{ivq_j}g(v)\Big|_{v=-{\cal V}/2}^{{\cal V}/2} \)


Eigenvetors of the coordinate operator

The normalized eigenvector of the coordinate operator with eigenvalue \( q_n \) in the \( v \) representation is
          \( \displaystyle e_{q_n}(v)=\frac{e^{-ivq_n}}{\sqrt{\cal V}}. \)

These functions are orthonormal in a discrete sense, i.e. 
          \( \displaystyle \int_{-{\cal V}/2}^{{\cal V}/2}dv\,e^*_{q_m}(v) e_{q_n}(v)=\text{sinc} \left[\frac{\cal V}{2}(q_m-q_n)\right] \)

Now, the conjugate to the coordinate eigenvector \( e_{q_n}(v) \) is

\( \displaystyle \)\( \displaystyle e_{q_n}(m) =\frac{1}{\sqrt{{\cal V}}}\int_{-{\cal V}/2}^{{\cal V}/2} e^{ivq_m}\frac{e^{-ivq_n}}{\sqrt{\cal V}}dv =\sqrt{\frac{{\cal V}}{2a}}\, \text{sinc}\left[\frac{\cal V}{2}(q_m-q_n)\right], \)

and the orthonormality between these vectors reads

\( \displaystyle \sum_{m=-N}^N\Delta_me^*_{q_j}(m)e_{q_n}(m)=\sum_{m=-N}^N \text{sinc}\left[\frac{\cal V}{2}(q_m-q_j)\right]\, \text{sinc}\left[\frac{\cal V}{2}(q_m-q_n)\right], \)
a result which becomes the Kronecker delta \( \delta_{jn} \) when \( {\cal V}\to\infty \).

Eigenvectors of the derivative operator

Now, the normalized eigenvector of the discrete derivative with eigenvalue \( v \) in the coordinate representation is
          \( \displaystyle e_{v}(n)=\frac{e^{ivq_n}}{\sqrt{\Delta_n(2N+1)}}, \)

and the orthogonality for these states reads
(31)          \( \displaystyle \sum_{n=-N}^N\Delta_ne^*_{v'}(n)e_{v}(n) =\frac{1}{2N+1}\sum_{n=-N}^Ne^{i(v-v')q_n} \mathop{\longrightarrow}_{N\to\infty} \delta(v-v'). \)

The conjugate function to \( e_{v}(n)=e^{ivq_n}/\sqrt{2N+1} \) is
          \( \displaystyle e_{v}(v') =\sum_{j=-N}^N\frac{e^{i(v'-v)q_j}} {\sqrt{{\cal V}(2N+1)}} , \)

This function approximates  a delta function with small noise in it.

Concluding remarks

This is a step more in the theory of discrete operators. It shows that it is possible to have discrete operators very similar to the usual operator of continuous variable theory. We will explore more things about this operator,  Things like its inverse, its use in obtaining self-adjoint extensions, for instance.

We have discussed a local approach to the finite differences first derivative. Another point of view is obtained by collecting the finite differences at each point of the mesh in a single matrix, the subject of another a future work.





  1. Press, W. H., Teukolsky, S. A., Vetterling, W. T. and Flannery, B. P.. Numerical recipes the Art of Numerical Computing, third edition; Cambridge University Press: Cambridge, UK, 207; pp. 1.
  2. William Squire; George Trapp; Using Complex Variables to Estimate Derivatives of Real Functions. SIAM Review 1997, 40, 110-112, 10.1137/s003614459631241x.
  3. Kopal, Z.. Numerical Analysis: With Emphasis on the Application of Numerical Techniques to Problems of Infinitesimal Calculus in Single Variable; John Wiley \& Sons Inc.: New York, 1961; pp. 1.
  4. Boole, George. A treatise on the calculus of finite differences; Cambridge University Press: New York, 1860, 2009; pp. 1.
  5. Jordan, Charles. Calculus of finite differences, 2nd ed; Chelsea Publishing Co.: New York, USA, 1950; pp. 1.
  6. Richardson, C. H.. An introduction to the calculus of finite differences; Van Nostrand Co., Inc.: New York, USA, 1954; pp. 1.
  7. Potts, R. B.; Differential and Difference Equations. Am Math Monthly 1982, 89, 402-407.
  8. Potts, R. B.; Ordinary and partial differences equations. J Austral. Math. Soc. Ser. B, 1986, 27, 488.
  9. Mickens, R.E.. Nonstandard finite difference models of differential equations; World Scientific: Singapore, 1994; pp. 1.
  10. Martínez-Pérez, A. and Torres-Vega, G.; Exact Finite Differences. The Derivative on Non Uniformly Spaced Partitions. Symmetry 2017, 9, 217.
  11. Potts, R. B.; Differential and Difference Equations. Am Math Monthly 1982, 89, 402-407.
  12. Potts, R. B.; Ordinary and partial differences equations. J Austral. Math. Soc. Ser. B, 1986, 27, 488.
  13. Mickens, R.E.. Nonstandard finite difference models of differential equations; World Scientific: Singapore, 1994; pp. 1.
  14. Martínez-Pérez, A. and Torres-Vega, G.; Exact Finite Differences. The Derivative on Non Uniformly Spaced Partitions. Symmetry 2017, 9, 217.
Contributor MDPI registered users' name will be linked to their SciProfiles pages. To register with us, please refer to :
View Times: 1.9K
Revisions: 3 times (View History)
Update Date: 24 Aug 2021
Video Production Service