Variational finite-difference representation of the kinetic energy operator

Abstract

A potential disadvantage of real-space-grid electronic structure methods is the lack of a variational principle and the concomitant increase of total energy with grid refinement. We show that the origin of this feature is the systematic underestimation of the kinetic energy by the finite difference representation of the Laplacian operator. We present an alternative representation that provides a rigorous upper bound estimate of the true kinetic energy and we illustrate its properties with a harmonic oscillator potential. For a more realistic application, we study the convergence of the total energy of bulk silicon using a real-space-grid density-functional code and employing both the conventional and the alternative representations of the kinetic energy operator.

3 pages, 3 figures, 1 table. To appear in Phys. Rev. B. Contribution for the 10th anniversary of the eprint server

Cited by in corpus (8)

• O(N) methods in electronic structure calculations
• Real-space electronic-structure calculations with timesaving double-grid technique
• Numerical Methods for the Inverse Problem of Density Functional Theory
• A review on non-relativistic fully numerical electronic structure calculations on atoms and diatomic molecules
• Comparison of variational real-space representations of the kinetic energy operator
• Properties of Quantum Systems via Diagonalization of Transition Amplitudes II: Systematic Improvements of Short-time Propagation
• Properties of Quantum Systems via Diagonalization of Transition Amplitudes I: Discretization Effects
• Engineering analog quantum chemistry Hamiltonians using cold atoms in optical lattices