The Partition of Unity Finite Element Method for the Schrödinger Equation

Overview

Abstract

A Schrödinger equation for the system’s wavefunctions in a parallelepiped unit cell subject to Bloch-periodic boundary conditions must be solved repeatedly in quantum mechanical computations to derive the materials’ properties. Recent studies have demonstrated how enriched finite element type Galerkin methods can substantially lower the number of degrees of freedom necessary to produce accurate solutions with respect to the standard plane-waves method. In particular, the flat-top partition of unity finite element method enriched with the radial eigenfunctions of the one-dimensional Schrödinger equation offers a very effective way of solving the three-dimensional Schrödinger eigenvalue problem. We investigate the theoretical properties of this approximation method, its well-posedness and stability, we prove its convergence and derive suitable bound for the h- and p-refinement in the L2 and energy norm for both the eigenvalues and the eigenfunctions. Finally, we confirm these theoretical results by applying this method to the eigenvalue problem of the one-electron Schrödinger equation with the harmonic potential, for which the exact solution is known.

Brief Biography

Francesca Gardini currently works at the Department of Mathematics "F. Casorati", University of Pavia. Francesca does research in Applied Mathematics and Fluid Dynamics. Her scientific interests deals with numerical methods for PDEs, finite element methods, mixed methods, mimetic finite differences, virtual element methods, approximation of eigenvalue problems, a posteriori error estimates, adaptive methods, interaction of fluid and structures. Her most recent publication is 'Virtual Element Method for Second Order Elliptic Eigenvalue Problems.