Real-space grid representation of momentum and kinetic energy operators for electronic structure calculations

J Comput Chem. 2018 Jul 30;39(20):1406-1412. doi: 10.1002/jcc.25208. Epub 2018 Mar 8.

Abstract

We show that the central finite difference formula for the first and the second derivative of a function can be derived, in the context of quantum mechanics, as matrix elements of the momentum and kinetic energy operators on discrete coordinate eigenkets |xn〉 defined on a uniform grid. Starting from the discretization of integrals involving canonical commutations, simple closed-form expressions of the matrix elements are obtained. A detailed analysis of the convergence toward the continuum limit with respect to both the grid spacing and the derivative approximation order is presented. It is shown that the convergence from below of the eigenvalues in electronic structure calculations is an intrinsic feature of the finite difference method. © 2018 Wiley Periodicals, Inc.

Keywords: basis sets and related methodology; computational methodology; density functional theory; electronic structure calculations; finite-difference schemes.