Partition of unity finite element method for quantum mechanical materials calculations

TL;DR

This paper introduces a partition of unity finite element method that efficiently solves quantum mechanical equations, reducing degrees of freedom significantly compared to traditional planewave methods, and is applicable to both metals and insulators.

Contribution

The paper develops a general PUFE method incorporating atomic physics into finite element analysis, improving efficiency for large-scale quantum materials calculations.

Findings

Achieves required accuracies with an order of magnitude fewer degrees of freedom than PW methods.

Accurately computes the equation of state for LiH, matching reference results.

Demonstrates efficiency for calculations involving localized atomic states and extreme conditions.

Abstract

The current state of the art for large-scale quantum-mechanical simulations is the planewave (PW) pseudopotential method, as implemented in codes such as VASP, ABINIT, and many others. However, since the PW method uses a global Fourier basis, with strictly uniform resolution at all points in space, it suffers from substantial inefficiencies in calculations involving atoms with localized states, such as first-row and transition-metal atoms, and requires significant nonlocal communications, which limit parallel efficiency. Real-space methods such as finite-differences and finite-elements have partially addressed both resolution and parallel-communications issues but have been plagued by one key disadvantage relative to PW: excessive number of degrees of freedom needed to achieve the required accuracies. We present a real-space partition of unity finite element (PUFE) method to solve the Kohn-Sham equations of density functional theory. In the PUFE method, we build the known atomic physics into the solution process using partition-of-unity enrichment techniques in finite element analysis. The method developed herein is completely general, applicable to metals and insulators alike, and particularly efficient for deep, localized potentials, as occur in calculations at extreme conditions of pressure and temperature. Full self-consistent Kohn-Sham calculations are presented for LiH, involving light atoms, and CeAl, involving heavy atoms with large numbers of atomic-orbital enrichments. We find that the new PUFE approach attains the required accuracies with substantially fewer degrees of freedom, typically by an order of magnitude or more, than the PW method. We compute the equation of state of LiH and show that the computed lattice constant and bulk modulus are in excellent agreement with reference PW results, while requiring an order of magnitude fewer degrees of freedom to obtain.