  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.
